| 查看: 265 | 回复: 0 | ||
[求助]
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 ] |
回复此楼» 猜你喜欢
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
filecode=后面第一个是大写字母
已经有4人回复
-
时间戳又变了8-15
已经有20人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有58人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
找到一些相关的精华帖子，希望有用哦~
matlab 求助程序
已经有5人回复
- matlab 语言错误
已经有4人回复
- 求助一个Matlab语句错误，Error==>plus Matrix dimensions must agree
已经有3人回复
- matlab的m文件关联问题
已经有17人回复
- matlab-出错提示Undefined method for input arguments type 'double'.
已经有5人回复
- 求助MATLAB做三维图问题
已经有4人回复
- matlab语句问题
已经有5人回复
- 求助 matlab隐函数画图问题
已经有4人回复
- matlab禁忌算法求解vrp问题
已经有10人回复
- matlab程序出错调试
已经有7人回复
- 求助一个matlab关于分类讨论的问题
已经有9人回复
- matlab编程中遇到的问题
已经有4人回复
- Matlab中save问题
已经有5人回复
- matlab符号函数求导及梯度的问题
已经有5人回复
- 跪求matlab的for循环语句的编写
已经有20人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- 请教关于matlab矩阵转换问题
已经有6人回复
- 急求一个matlab编码，各位如果知道一定要帮我一把啊
已经有10人回复
- 申请TUD语言问题
已经有25人回复
- 关于Matlab作图的问题
已经有14人回复
- matlab下bvp4c求解非线方程问题
已经有5人回复
- matlab非线性参数拟合问题
已经有7人回复
- MATLAB 画图小问题
已经有11人回复
- 100BB急求matlab编程中的问题解决，如何修改这个程序呢？
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
