| 5 | 1/1 | 返回列表 |
| 查看: 869 | 回復: 7 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[求助]
matlab程序調(diào)試
|
|||
|
如題,程序如下,不知道為什么會報錯,希望能幫助調(diào)試一下,謝謝! 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) |
|
本帖內(nèi)容被屏蔽 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 283分材料與化工求調(diào)劑 +11 | 羅KAKA 2026-04-02 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 一志愿哈工大 材料與化工 +30 | 洛神哥哥 2026-03-31 | 32/1600 |
|
|
[考研] 309求調(diào)劑 +8 | 呆菇不是戴夫 2026-04-02 | 8/400 |
|
|
[考研] 324求調(diào)劑 +6 | 想上學求調(diào) 2026-04-02 | 6/300 |
|
|
[考研] 286分調(diào)劑 +20 | Faune 2026-03-30 | 22/1100 |
|
|
[考研] 262求調(diào)劑 +5 | 勵志一定發(fā)文章 2026-04-02 | 6/300 |
|
|
[考研] 311求調(diào)劑 +9 | 勇敢的小吳 2026-04-02 | 9/450 |
|
|
[考研]
|
shzhou_ 2026-04-01 | 6/300 |
|
|
[考研] 各位老師好,我的一志愿為北京科技大學085601材料專碩 +13 | Koxui 2026-03-28 | 13/650 |
|
|
[考研] 求調(diào)劑,一志愿南京師范大學計算機專碩,初試373,六級通過, +3 | 計算機追夢人 2026-04-01 | 3/150 |
|
|
[考研] 291求調(diào)劑 +20 | Y-cap 2026-03-29 | 25/1250 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +7 | 烤麥芽 2026-03-27 | 10/500 |
|
|
[考研] 材料專碩調(diào)劑 +17 | 椰椰。 2026-03-29 | 17/850 |
|
|
[考研] 0710生物學求調(diào)劑! +6 | 敘述文 2026-03-31 | 6/300 |
|
|
[考研] 0856材料化工調(diào)劑 總分330 +18 | zhubinhao 2026-03-27 | 18/900 |
|
|
[考研] 085601 329分調(diào)劑 +6 | yzsa12 2026-03-31 | 6/300 |
|
|
[考研] 本科211安全工程,初試290分,求調(diào)劑 +3 | 2719846834 2026-03-28 | 3/150 |
|
|
[考研] 調(diào)劑求院校招收 +7 | 鶴鯨鴿 2026-03-28 | 7/350 |
|
|
[考研] 330一志愿中國海洋大學 化學工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|