数学建模社区-数学中国
标题:
Maple的for循环求助
[打印本页]
作者:
No-one
时间:
2020-3-27 14:25
标题:
Maple的for循环求助
本帖最后由 No-one 于 2020-3-27 14:26 编辑
# T Z2 T3 D; q6 v% v8 J8 Q
% _3 C% O0 h, E5 y/ A% q6 ^
想通过Maple的for循环来把数据依次写入到一个零矩阵里，但结果是产生了一个table,矩阵里的0没有被替换，请问大佬们该如何修改。谢谢啦
7 P3 A- N7 b8 D8 j1 s
#######################################
; r3 S2 t* ]; U( V7 l
Qlist := [seq(i, i = 1 .. 2)];
) ]/ H( d" v3 o' s- c7 i" o# N
klist := [seq(.1*i, i = 1 .. 2)];
% m3 X3 y# F* ~) t7 Z9 T3 p% j* e
WL_Matrix := Matrix(nops(Qlist), nops(klist), 0);
* n( T' \3 E' V6 ~* c( P! o) e
n := nops(Qlist);
2 j0 |4 D7 h0 e! P3 u5 {8 |8 |
9 n) h3 Q5 b& N1 l1 A4 c0 O9 u6 M
for i to n do
6 l3 g1 `6 m4 W9 ~: t* L, O
c) k3 G, U/ v: ?
for j to n do
# k8 b" n6 ], Y( o4 G
3 ~" U1 J! X/ r0 T3 N" X3 i
eq5 := [solve(subs([Q = Qlist
, k = klist[j]
], W^2*W0*k+Q^2) = 0, W)];
% Z0 d2 l5 n' ^5 N) T
. Y7 D0 b: {" {7 v& @- W; N* K
WL_Marix[i, j] := eq5[1]; print(WL_Marix);
7 C/ R* n& @4 L+ }/ ?
1 U3 @/ }! Y6 ~% b1 x( q
eq5 := assign['eq5']
& O. M5 [6 T" b
0 @( ]. F8 _3 a' @9 j
end do
' e" B2 ]9 N9 E/ ^
2 D N5 Q: M! Z; w
end do;
% @. D7 C. y& T( Y
#################################
) s+ O' d( Z) g: X
6 ]$ H$ M; u$ i& P9 K
; T4 ^3 d1 ?3 X; L8 u$ W
) L+ O$ Y/ L. P3 D" }! N" P
. ~: a: b4 S7 G: \0 F7 p/ q# R
4 }2 K) C* t8 i. B. {, Q4 {
1 u/ k( V4 }/ Q) G+ o0 c
* R8 O1 @9 i" t8 f; ]
# R) H3 i% P9 a3 d$ p' j8 u
Capture1.PNG
(9.48 KB, 下载次数: 410)
2020-3-27 14:24 上传
点击文件名下载附件
Maple18运行结果
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
