数学建模社区-数学中国

标题: Maple的for循环求助 [打印本页]

作者: No-one    时间: 2020-3-27 14:25
标题: Maple的for循环求助
本帖最后由 No-one 于 2020-3-27 14:26 编辑 . u4 T. c& g( _( d+ l$ G

7 l- F' P6 p1 o3 ]2 C想通过Maple的for循环来把数据依次写入到一个零矩阵里,但结果是产生了一个table,矩阵里的0没有被替换,请问大佬们该如何修改。谢谢啦- S- o* d; |0 P# l! s
#######################################: k( b$ E) N4 a
Qlist := [seq(i, i = 1 .. 2)];
: M. O: R# p, }# h7 B5 U8 Yklist := [seq(.1*i, i = 1 .. 2)];1 A- Z$ J0 f  n/ }% h9 y  e
WL_Matrix := Matrix(nops(Qlist), nops(klist), 0);
& P, `* \/ W3 O3 [  H; Dn := nops(Qlist);
7 h& H- s* P8 m+ M9 B$ E4 g0 G. I6 k7 S& r
for i to n do8 e9 `4 ~. t7 L8 Z9 b- K( V# J6 F5 R

/ W) p7 S) p0 K for j to n do
4 C* }  t6 b  T  i# T' A/ {
2 r; [& m: z6 F- q# Z2 Heq5 := [solve(subs([Q = Qlist, k = klist[j]], W^2*W0*k+Q^2) = 0, W)];
0 @: U% g4 x1 P6 i# y$ L
% e" _; ^- @2 G- w3 x- nWL_Marix[i, j] := eq5[1]; print(WL_Marix); " ]1 H6 V7 }# A% w6 L' d

' ^: _# @* m/ S. Keq5 := assign['eq5'] % t( j# d9 D6 R5 J% i7 ^7 r
' p) L0 K1 |8 Y
end do / ~6 h+ p4 ~$ _5 J& d, U& d

7 B6 V# Q" v! e2 R5 ~end do;
: ^! c4 x  G/ |3 A0 M#################################4 o8 x( H0 k% ?. B% I1 ?, D7 Z

8 K/ F# C- P) c- w- x1 S! X5 a! j  ~

" C& _, w0 a, g' ]7 e- h6 J, V

# H: U" E, p  c, z" E* P
$ k! K, f2 ]+ {4 N) a! ]2 O( \- _: u5 q- u$ T- m4 o
2 |) q" W+ P  C. E0 _. x* \5 i! r

Capture1.PNG (9.48 KB, 下载次数: 454)

Maple18运行结果

Maple18运行结果






欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5