| 查看: 895 | 回复: 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) |
回复此楼» 猜你喜欢
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有72人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
考博自荐
已经有3人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
5