| 查看: 1079 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
matlab程序调试
|
|||
|
如题，程序如下，不知道为什么会报错，希望能帮助调试一下，谢谢！ clear clc a1=5.756; a2=0.0983; a3=0.2020; a4=189.32; a5=12.52; a6=1.32e-2; b1=2.860e-6; b2=4.700e-8; b3=6.113e-8; b4=1.516e-4; T=25; f=(T-24.5).*(T+570.82); A=28.0:0.01:32.0; lamdap=1.064; lamdas=1.4:0.01:5.4; lamdai=1./(1./lamdap-1./lamdas); c=1; c2=1; for i=1:length(A) m=1; for n=1:length(lamdas) if lamdai(n)>lamdas(n); lamdass(m)=lamdas(n); lamdaii(m)=lamdai(n); np(m)=(a1+b1*f+(a2+b2*f)/(lamdap^2-(a3+b3*f)^2)+(a4+b4*f)/(lamdap^2-a5^2)-a6*lamdap^2)^(1/2); ns(m)=(a1+b1*f+(a2+b2*f)/(lamdass(m)^2-(a3+b3*f)^2)+(a4+b4*f)/(lamdass(m)^2-a5^2)-a6*lamdass(m)^2)^(1/2); ni(m)=(a1+b1*f+(a2+b2*f)/(lamdaii(m)^2-(a3+b3*f)^2)+(a4+b4*f)/(lamdaii(m)^2-a5^2)-a6*lamdaii(m)^2)^(1/2); w(m)=abs(np(m)./lamdap-ns(m)./lamdass(m)-ni(m)./lamdaii(m)-1./A); m=m+1; end end if min(w)<10^-5 [x1,y1]=min(w); A1(c)=A(i); lamdas1(c)=lamdass(y1); lamdai1(c)=lamdaii(y1); lamdap2=lamdas1(c); lamdas2=(lamdap2+0.001):0.0001:6; lamdai2=1./(1./lamdap2-1./lamdas2); c=c+1; m2=1; for n2=1:length(lamdas2) if lamdai2(n2)>lamdas2(n2) lamdass2(m2)=lamdas2(n2); lamdaii2(m2)=lamdai2(n2); ns2(m2)=(a1+b1*f+(a2+b2*f)/(lamdass2(m2)^2-(a3+b3*f)^2)+(a4+b4*f)/(lamdass2(m2)^2-a5^2)-a6*lamdass2(m2)^2)^(1/2); ni2(m2)=(a1+b1*f+(a2+b2*f)/(lamdaii2(m2)^2-(a3+b3*f)^2)+(a4+b4*f)/(lamdaii2(m2)^2-a5^2)-a6*lamdaii2(m2)^2)^(1/2); np2(m2)=(a1+b1*f+(a2+b2*f)/(lamdap2^2-(a3+b3*f)^2)+(a4+b4*f)/(lamdap2^2-a5^2)-a6*lamdap2^2)^(1/2); w2(m2)=abs(np2(m2)/lamdap2-ns2(m2)/lamdass2(m2)-ni2(m2)/lamdaii2(m2)-1./A); m2=m2+1; end end if min(w2)<10^-6.07 [x2,y2]=min(w2); A2(c2)=A(i); lamdasss2(c2)=lamdass2(y2); lamdaiii2(c2)=lamdaii2(y2); c2=c2+1; end end end plot(A1,lamdas1,A1,lamdai1,A2,lamdasss2,A2,lamdaiii2) |
回复此楼» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有208人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有112人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有69人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab程序出错调试
已经有7人回复
- 搞地学的朋友，请问怎么用matlab做Mann-Kendall趋势检验啊？
已经有10人回复
- matlab调试问题
已经有3人回复
- 求帮忙啊~ matlab一多项式结果太长，怎么写到txt里？
已经有14人回复
- matlab调试问题
已经有3人回复
- while语句死循环了，求教
已经有6人回复
- 求助matlab编程
已经有13人回复
- 【求助】VS2005数值计算程序调试【急求！！！！加金币】
已经有9人回复
★
感谢参与，应助指数 +1
csgt0: 金币+1, 谢谢 2013-01-10 10:35:48
感谢参与，应助指数 +1
csgt0: 金币+1, 谢谢 2013-01-10 10:35:48
|
本帖内容被屏蔽 |
7楼2013-01-09 20:15:58
2楼2013-01-08 10:10:51
3楼2013-01-08 13:46:49
4楼2013-01-08 19:49:12
