作者easy1612 (up)
看板MATLAB
标题[问题] 回圈取值
时间Sun Jan 1 22:10:49 2012
请问各位大大
小弟得程式码如下
1.clc;
2.n=2.168;
3.X=[];
4.b=[];
5.for L=0:11
6.M21=[1 L ; 0 1];
7.M22=[1 0; 0 1];
8.M23=[1 11-L; 0 1];
9.M24=[1 0 ; 0 1/n];
10.M25=[1 8/n ;0 1 ];
11.M26=[1 0 ;0 n];
12.M27=[1 5; 0 1];
13.M28=[1 0; -2/150 1];
14.M2t=M21*M22*M21*M23*M24*M25*M26*M27*M28*M27*M24*M25*M26*M23;
15.A=M2t(1,1);
16.B=M2t(1,2);
17.C=M2t(2,1);
18.D=M2t(2,2);
19.det(M2t);
20.q1=(A-D-((D-A)^2+4*B*C)^0.5)/(2*C);
21.q11=1/q1;
22.R=1/real(q11);
23.W=(-(1.064*(10^-18))/(pi*(imag(q11))))^0.5;
24.b=[L,W]
25.X=[X,b'];
26.end
27.X
理想为回圈内取出L和W
问题出在24行
出来的b值只有L，W都是0
确定W是有值
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.115.42.207
※ 编辑: easy1612 来自: 140.115.42.207 (01/01 22:13)
1F：推 sunev:把23行的分号去掉，看看情况，再把21的也去掉。 01/01 23:04
2F：→ easy1612:修正过後结果与修正前相同 01/02 00:28
3F：→ Gwaewluin:结果当然是相同的，去掉分号只是要你检查途中的计算结果 01/02 07:32