| 查看: 1074 | 回复: 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论文润色/翻译怎么收费?
已经有262人回复
-
澳门大学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人回复
2楼2013-01-08 10:10:51
3楼2013-01-08 13:46:49
4楼2013-01-08 19:49:12
lgycjpcqu
金虫 (正式写手)
- 应助: 35 (小学生)
- 金币: 566.8
- 散金: 411
- 帖子: 371
- 在线: 176.9小时
- 虫号: 1327328
- 注册: 2011-06-20
- 性别: GG
- 专业: 电化学
【答案】应助回帖
★ ★ ★
感谢参与，应助指数 +1
csgt0: 金币+1, 谢谢 2013-01-09 09:22:05
惜梦寻草: 金币+2 2013-01-09 18:57:35
感谢参与，应助指数 +1
csgt0: 金币+1, 谢谢 2013-01-09 09:22:05
惜梦寻草: 金币+2 2013-01-09 18:57:35
|
A向量中共有401个元素，m的取值范围由于判定语句 if lamdai(n)>lamdas(n); 的存在必然会小于401（如果等于此判定就没有意义） 所以错误提示说，w和A中的元素要求相等。 如果语句 w(m)=abs(np(m)./lamdap-ns(m)./lamdass(m)-ni(m)./lamdaii(m)-1./A); 的意思说当上面的判定成立时，用A向量对应位置的元素计算w的数值的话，改成1/A（m）即可。 我也是初学者，呵呵希望对你有帮助！ |
5楼2013-01-09 08:58:18
6楼2013-01-09 18:57:04
★
感谢参与，应助指数 +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
lgycjpcqu
金虫 (正式写手)
- 应助: 35 (小学生)
- 金币: 566.8
- 散金: 411
- 帖子: 371
- 在线: 176.9小时
- 虫号: 1327328
- 注册: 2011-06-20
- 性别: GG
- 专业: 电化学
8楼2013-01-09 22:20:08
