数学建模社区-数学中国
标题:
Maple的for循环求助
[打印本页]
作者:
No-one
时间:
2020-3-27 14:25
标题:
Maple的for循环求助
本帖最后由 No-one 于 2020-3-27 14:26 编辑
$ o1 D3 U' o# U) e5 l
9 e+ I2 C% r$ ~5 T3 r7 e
想通过Maple的for循环来把数据依次写入到一个零矩阵里，但结果是产生了一个table,矩阵里的0没有被替换，请问大佬们该如何修改。谢谢啦
, `: S9 E. v( R1 b3 v6 b
#######################################
/ V' k! A+ l( p) t
Qlist := [seq(i, i = 1 .. 2)];
% A R+ L& i' |7 Y5 w
klist := [seq(.1*i, i = 1 .. 2)];
* _1 Y+ t! k. G7 c9 j5 V
WL_Matrix := Matrix(nops(Qlist), nops(klist), 0);
& @. @3 m. k+ ^: V' B
n := nops(Qlist);
( u% g- f3 k: N/ Y
! F# [% }4 T9 r' ]) B
for i to n do
! Q5 N1 R" G( ^6 Q9 N" G
! X+ _ s. X F, e: [+ K& Q
for j to n do
8 N3 v5 S) {/ {5 Y; V5 l
/ ^5 L, S0 T& J5 q" P3 \
eq5 := [solve(subs([Q = Qlist
, k = klist[j]
], W^2*W0*k+Q^2) = 0, W)];
+ n8 Z; G, ~; `) ^6 v3 G
: b- v* R! J! [) K
WL_Marix[i, j] := eq5[1]; print(WL_Marix);
a; X8 p$ ]6 q5 m/ f& k4 v
N0 p; b6 U& x( e T& _% p
eq5 := assign['eq5']
) i. f6 ?/ | i$ _( B0 L! c
# X' \) [& }5 u1 D6 p
end do
. A- f: X' i J% [+ D: V; [
" t2 Y' }1 \; | A2 }: u
end do;
2 K6 N. j" ?% w% y% m1 R' p
#################################
( M( S1 m( `4 M( ^4 q- n' ~3 s
, [* |+ f- ] \; q& m, [2 S
$ L; C0 `0 ~ ~9 w7 B# A
& l4 c5 l/ G4 F" D+ O
3 g7 J2 E/ [8 z0 N+ \, y
$ G1 e! h; \% p, Z, x9 N
, f; F2 m# Z) n" M* o( f
" H8 a; @! F; T( k! Y; Z! ?" |5 E
g! @- X6 e9 i4 C$ R- _$ ^
Capture1.PNG
(9.48 KB, 下载次数: 467)
2020-3-27 14:24 上传
点击文件名下载附件
Maple18运行结果
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
