作业帮 > 综合 > 作业

用Mathematica 怎样写语句 从文件中读入一维数组,怎样输出数组到另一个文件中?

来源:学生作业帮 编辑:拍题作业网作业帮 分类:综合作业 时间:2024/05/01 17:01:44
用Mathematica 怎样写语句 从文件中读入一维数组,怎样输出数组到另一个文件中?
是把数组存到一个txt文件吗?然后怎样写语句读入这个数组进入程序运算,之后有怎样把计算结果(一维数组)输出到另外一个txt文件?假设我把这两个txt文件存到与mathematica notebook文件的同一个文件夹下.
1.txt 中存的是d的以为数组 代表d分别取5个值,通过以下公式运算得出相应的5个v的实数解,并将v的计算结果依次存到2.txt文件中.
Solve[d*1.309*10^(-3) (0.0186 + v^2/2) ==
2.36086*10^(-8)* v^1.4 + 10^(-7),v]
d={0.001,0.00123,0.00138,0.0015,0.0016}
(**生成一维数组,存到C:/Users/zhou/Desktop/11.txt **)
Data=Table[w,{w,0,1,0.01}]
Export["C:/Users/zhou/Desktop/11.txt", Data, "Table"]
(**从文件中读入一维数组 **)
snum = OpenRead["C:/Users/zhou/Desktop/11.txt"]
snum1 = ReadList[snum, Number]
Close[snum]
(**计算 **)
snum2 = snum1/3
(**把计算结果(一维数组)输出到C:/Users/zhou/Desktop/22.txt **)
Export["C:/Users/zhou/Desktop/22.txt", snum2, "Table"]
再问: aa = OpenRead["C:/Users/DONNA/Desktop/1.txt"] bb = ReadList[aa, Number] Close[aa] d = bb Solve[v*d == 0.2, v] cc = v Export["C:/Users/DONNA/Desktop/2.txt", cc, "Table"] Close[cc] 请帮忙看看错在哪里,为什么不能输出想要的计算结果?谢谢!文件的读入语句是对的。是Solve[v*d == 0.2, v] 和cc = v这俩个语句的问题,这只是一个简单的方程,但是如果计算复杂的方程(隐函数),这种情况计算语句怎么写?
再答: 请提供 1.txt 如果1.txt内的数据为一维数组,Solve[v*d == 0.2, v]的解集为{}。 你想解的是什么样的方程? 给出一个你想计算复杂的方程(隐函数)
再问: 1.txt 中存的是d的以为数组 代表d分别取5个值, 通过以下公式运算得出相应的5个v的实数解,并将v的计算结果依次存到2.txt文件中。谢谢! Solve[d*1.309*10^(-3) (0.0186 + v^2/2) == 2.36086*10^(-8)* v^1.4 + 10^(-7), v]
再答: 你要解的方程是 Solve[d*1.309*10^(-3)*(0.0186 + v^2/2) == 2.36086*10^(-8)* v^1.4 + 10^(-7), v] ? d分别取5个值是多少? Solve[d*1.309*10^(-3)*(0.0186 + v^2/2) == 2.36086*10^(-8)* v^1.4 + 10^(-7), v, Reals] 的通解为: {{v -> ConditionalExpression[ Root[-3.125*10^28 + 3.80428*10^31 d - 1.85249*10^34 d^2 + 4.51032*10^36 d^3 - 5.49073*10^38 d^4 + 2.6737*10^40 d^5 + (1.02266*10^33 d - 9.95961*10^35 d^2 + 3.63736*10^38 d^3 - 5.90401*10^40 d^4 + 3.59369*10^42 d^5) #1^2 + (-1.33866*10^37 d^2 + 9.77785*10^39 d^3 - 2.38065*10^42 d^4 + 1.93209*10^44 d^5) #1^4 + (8.76151*10^40 d^3 - 4.2664*10^43 d^4 + 5.19379*10^45 d^5) #1^6 - 2.29193*10^25 #1^7 + (-2.8672*10^44 d^4 + 6.9809*10^46 d^5) #1^8 + 3.75317*10^47 d^5 #1^10 &, 1], 0.00410721 < d < 0.00410722]}, {v -> ConditionalExpression[ Root[-3.125*10^28 + 3.80428*10^31 d - 1.85249*10^34 d^2 + 4.51032*10^36 d^3 - 5.49073*10^38 d^4 + 2.6737*10^40 d^5 + (1.02266*10^33 d - 9.95961*10^35 d^2 + 3.63736*10^38 d^3 - 5.90401*10^40 d^4 + 3.59369*10^42 d^5) #1^2 + (-1.33866*10^37 d^2 + 9.77785*10^39 d^3 - 2.38065*10^42 d^4 + 1.93209*10^44 d^5) #1^4 + (8.76151*10^40 d^3 - 4.2664*10^43 d^4 + 5.19379*10^45 d^5) #1^6 - 2.29193*10^25 #1^7 + (-2.8672*10^44 d^4 + 6.9809*10^46 d^5) #1^8 + 3.75317*10^47 d^5 #1^10 &, 2], 0 < d < 0.00410721 || 0.00410721 < d < 0.00410722]}}
再问: d={0.001, 0.00123, 0.00138, 0.0015, 0.0016} 求出v的 5 个相应实数解后,输出到 2.txt 文件里。
再答: sol = Solve[d1*1.309*10^(-3)*(0.0186 + v^2/2) == 2.36086*10^(-8)* v^1.4 + 10^(-7), v, Reals][[2]]; d={0.001, 0.00123, 0.00138, 0.0015, 0.0016} ; rr = Simplify[v /. sol /. Table[{d1 -> d[[i]]}, {i, 5}]] Export["C:/Users/zhou/Desktop/22.txt", rr, "Table"]
再问: 非常感谢zzz700的热心帮忙!!已经提高了悬赏:) 还有一个问题请教: sol = Solve[d1*1.309*10^(-3)*(0.0186 + v^2/2) == 2.36086*10^(-8)* v^1.4 + 10^(-7), v, Reals][[2]]; 这句中的最后的[[2]]是代表什么意思?
再答: sol = Solve[d1*1.309*10^(-3)*(0.0186 + v^2/2) == 2.36086*10^(-8)* v^1.4 + 10^(-7), v, Reals] 的解有两个。 你的d={0.001, 0.00123, 0.00138, 0.0015, 0.0016} 0 < d < 0.00410721 最后的[[2]]是代表取第二个解。