24小时热门版块排行榜

返回列表

f2012

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 237.7
散金: 8
帖子: 92
在线: 31.2小时
虫号: 1682710
注册: 2012-03-11
专业: 钢铁冶金

[求助] matlab 语言错误

>> T=1700;
p=[45584.0871,39280.9406,10069.7385,34484.3661]; %c1=b1/mn  c2=a1/si
q=[-28.50974,-6.8812,-2.1133,-17.838];

k=zeros(1,4);
for i=1:1:4;
k(i)=exp(p(i)/T+q(i));
end

N1=zeros(1,26);
N2=zeros(1,26);
N3=zeros(1,26);
N4=zeros(1,26);
N5=zeros(1,26);
N6=zeros(1,26);
S=zeros(1,3);
U=zeros(1,3);
g=zeros(1,2);
h=zeros(1,2);
x=zeros(1,2);

a=[0.05 0.057 0.1 0.111 0.15 0.18 0.2 0.25 0.257 0.3 0.334 0.35 0.4 0.46 0.5
0.56 0.6 0.625 0.65 0.666 0.7 0.75 0.8 0.85 0.9 0.95];
b=[0.95 0.943 0.9 0.889 0.85 0.82 0.8 0.75  0.743 0.7 0.666 0.65 0.6 0.54 0.5
0.44 0.4 0.375 0.35 0.334 0.3 0.25 0.2 0.15 0.1 0.05];
E=0.0000001;

for i=1:1:26
if i==1;
g(1)=0.02*a(i);
h(1)=b(i);
elseif i>1;
g(1)=N1(i-1)*a(i)/a(i-1);
h(1)=N2(i-1)*b(i)/b(i-1);
end

while j<100000;
g(2)=g(1)+E;
h(2)=h(1)+E;
x(1)=g(1);
x(2)=h(1);
W=x(1)+x(2)+k(1)*(x(1)^3)*x(2)+k(2)*(x(1)^5)*(x(2)^3)+k(3)*x(1)* x(2)+k(4)*
(x(1)^11)* (x(2)^19)-1;
V=b(i)*x(1)-a(i)*x(2)+(3*b(i)-a(i))*k(1)*(x(1)^3)* x(2)+ (5*b(i)- 3*a(i))*k
(2)*(x(1)^5)* (x(2)^3)+(b(i)-a(i))*k(3)*x(1)*x(2)+(11*b(i)-19*a(i))*k(4)*(x
(1)^11)* (x(2)^19);
S(1)=W;
U(1)=V;

x(1)=g(2);
x(2)=h(1);
W=x(1)+x(2)+k(1)*(x(1)^3)*x(2)+k(2)*(x(1)^5)*(x(2)^3)+k(3)*x(1)* x(2)+k(4)*
(x(1)^11)* (x(2)^19)-1;
V=b(i)*x(1)-a(i)*x(2)+(3*b(i)-a(i))*k(1)*(x(1)^3)* x(2)+ (5*b(i)- 3*a(i))*k
(2)*(x(1)^5)* (x(2)^3)+(b(i)-a(i))*k(3)*x(1)*x(2)+(11*b(i)-19*a(i))*k(4)*(x
(1)^11)* (x(2)^19);
S(2)=W;
U(2)=V;

x(1)=g(1);
x(2)=h(2);
W=x(1)+x(2)+k(1)*(x(1)^3)*x(2)+k(2)*(x(1)^5)*(x(2)^3)+k(3)*x(1)* x(2)+k(4)*
(x(1)^11)* (x(2)^19)-1;
V=b(i)*x(1)-a(i)*x(2)+(3*b(i)-a(i))*k(1)*(x(1)^3)* x(2)+ (5*b(i)- 3*a(i))*k
(2)*(x(1)^5)* (x(2)^3)+(b(i)-a(i))*k(3)*x(1)*x(2)+(11*b(i)-19*a(i))*k(4)*(x
(1)^11)* (x(2)^19);
S(3)=W;
U(3)=V;

S(2)=S(2)-S(1);
U(2)=U(2)-U(1);
S(3)=S(3)-S(1);
U(3)=U(3)-U(1);
D=(S(2)*U(3)-S(3)*U(2))/(E^2);

x(1)=g(1)+(-S(1)*U(3)+S(3)*U(1))/E/D;
x(2)=h(1)+(-S(2)*U(1)+S(1)*U(2))/E/D;
W=x(1)+x(2)+k(1)*(x(1)^3)*x(2)+k(2)*(x(1)^5)*(x(2)^3)+k(3)*x(1)* x(2)+k(4)*
(x(1)^11)* (x(2)^19)-1;
V=b(i)*x(1)-a(i)*x(2)+(3*b(i)-a(i))*k(1)*(x(1)^3)* x(2)+ (5*b(i)- 3*a(i))*k
(2)*(x(1)^5)* (x(2)^3)+(b(i)-a(i))*k(3)*x(1)*x(2)+(11*b(i)-19*a(i))*k(4)*(x
(1)^11)* (x(2)^19);

if abs(W)>0.0000005;
g(1)=x(1);
h(1)=x(2);
elseif abs(V)>0.0000005
g(1)=x(1);
h(1)=x(2);
else
break;
end
j=j+1;

end
N1(i)=x(1);
N2(i)=x(2);
N3(i)=k(1)*(N1(i)^3)*N2(i);
N4(i)=k(2)*(N1(i)^5)*(N2(i)^3);
N5(i)=k(3)*N1(i)*N2(i);
N6(i)=k(4)*(N1(i)^11)*(N2(i)^19);

end
activety=[N1;N2;N3;N4;N5;N6;]

>> 运行结果都是0

[ Last edited by f2012 on 2013-10-31 at 23:22 ]

回复此楼