数学建模社区-数学中国
标题:
Maple的for循环求助
[打印本页]
作者:
No-one
时间:
2020-3-27 14:25
标题:
Maple的for循环求助
本帖最后由 No-one 于 2020-3-27 14:26 编辑
& g2 ?9 S4 D! v9 M7 o% p4 X
& x7 c- I; I) Y+ v6 ~/ ^
想通过Maple的for循环来把数据依次写入到一个零矩阵里，但结果是产生了一个table,矩阵里的0没有被替换，请问大佬们该如何修改。谢谢啦
3 ]5 a( \+ g( J
#######################################
3 k Q0 S6 ^! J/ H
Qlist := [seq(i, i = 1 .. 2)];
+ {5 B+ N4 D1 W9 H; \7 J
klist := [seq(.1*i, i = 1 .. 2)];
& B# Y: y; B' p, c O
WL_Matrix := Matrix(nops(Qlist), nops(klist), 0);
8 P8 s. }% s: Z, D7 T+ G m
n := nops(Qlist);
# g( S1 r# j" Z- J
, k# x9 Q8 {0 R8 e" f
for i to n do
3 f7 N2 n5 z; J1 L7 T( m
/ _, }, Z- R7 i8 Z) f3 g
for j to n do
; b: g* R3 t2 t# e) s
9 Y9 H ~9 W* c+ b* @
eq5 := [solve(subs([Q = Qlist
, k = klist[j]
], W^2*W0*k+Q^2) = 0, W)];
1 h0 |/ p- X) S; i \6 y- k
. }' T& y' E" t/ q( e
WL_Marix[i, j] := eq5[1]; print(WL_Marix);
" F j& D' r; ~" R1 D
* l. t5 d% V) e% X4 Q' x2 o, G) K
eq5 := assign['eq5']
# O& b+ R5 Y- t5 f
9 v" m/ v* @ h/ r$ H$ B1 B% M% e
end do
% V" P# C, ^9 `9 @* w- y
5 m V# H; i! U- T' p( k. p; u
end do;
' @/ @) Q) X$ @, u L% `. ~! {
#################################
- o! I- \3 `$ T( F. n6 g
3 V9 K) P4 ~+ @ G9 R7 c1 f; |
9 m8 [ v8 Y' W* ^% H
1 I& A/ Y, t5 A9 ~
# o! A$ ~8 w" c, f+ Q
7 D: r' \$ }+ A5 K$ Y
& U+ n: `" F' l7 R+ \4 ]2 C
; r& R* F1 b! z
1 Z2 e( |; E. t" _, \0 p. A
Capture1.PNG
(9.48 KB, 下载次数: 421)
2020-3-27 14:24 上传
点击文件名下载附件
Maple18运行结果
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
