 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 867 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
惜夢(mèng)尋草新蟲 (初入文壇)
|
[求助]
matlab程序調(diào)試
|
|
|
如題，程序如下，不知道為什么會(huì)報(bào)錯(cuò)，希望能幫助調(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) |
回復(fù)此樓» 猜你喜歡
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有53人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
-
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有1人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- matlab程序出錯(cuò)調(diào)試
已經(jīng)有7人回復(fù)
- 搞地學(xué)的朋友，請(qǐng)問怎么用matlab做Mann-Kendall趨勢(shì)檢驗(yàn)啊？
已經(jīng)有10人回復(fù)
- matlab調(diào)試問題
已經(jīng)有3人回復(fù)
- 求幫忙啊~ matlab一多項(xiàng)式結(jié)果太長(zhǎng)，怎么寫到txt里？
已經(jīng)有14人回復(fù)
- matlab調(diào)試問題
已經(jīng)有3人回復(fù)
- while語句死循環(huán)了，求教
已經(jīng)有6人回復(fù)
- 求助matlab編程
已經(jīng)有13人回復(fù)
- 【求助】VS2005數(shù)值計(jì)算程序調(diào)試【急求�。。�！加金幣】
已經(jīng)有9人回復(fù)
惜夢(mèng)尋草
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 875.5
- 帖子: 25
- 在線: 5.7小時(shí)
- 蟲號(hào): 2118505
- 注冊(cè): 2012-11-10
- 性別: MM
- 專業(yè): 光學(xué)
4樓2013-01-08 19:49:12
【答案】應(yīng)助回帖
★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+1, 謝謝 2013-01-08 13:58:02
惜夢(mèng)尋草: 金幣+2 2013-01-09 18:57:26
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+1, 謝謝 2013-01-08 13:58:02
惜夢(mèng)尋草: 金幣+2 2013-01-09 18:57:26
|
In an assignment A(I) = B, the number of elements in B and I must be the same. Error in Untitled (line 30) w(m)=abs(np(m)./lamdap-ns(m)./lamdass(m)-ni(m)./lamdaii(m)-1./A); >> 估計(jì)問題出在那個(gè) 1./A 上  |
2樓2013-01-08 10:10:51
3樓2013-01-08 13:46:49
lgycjpcqu
金蟲 (正式寫手)
- 應(yīng)助: 35 (小學(xué)生)
- 金幣: 566.8
- 散金: 411
- 帖子: 371
- 在線: 176.9小時(shí)
- 蟲號(hào): 1327328
- 注冊(cè): 2011-06-20
- 性別: GG
- 專業(yè): 電化學(xué)
【答案】應(yīng)助回帖
★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+1, 謝謝 2013-01-09 09:22:05
惜夢(mèng)尋草: 金幣+2 2013-01-09 18:57:35
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+1, 謝謝 2013-01-09 09:22:05
惜夢(mèng)尋草: 金幣+2 2013-01-09 18:57:35
|
A向量中共有401個(gè)元素，m的取值范圍由于判定語句 if lamdai(n)>lamdas(n); 的存在必然會(huì)小于401（如果等于此判定就沒有意義） 所以錯(cuò)誤提示說，w和A中的元素要求相等。 如果語句 w(m)=abs(np(m)./lamdap-ns(m)./lamdass(m)-ni(m)./lamdaii(m)-1./A); 的意思說當(dāng)上面的判定成立時(shí)，用A向量對(duì)應(yīng)位置的元素計(jì)算w的數(shù)值的話，改成1/A（m）即可。 我也是初學(xué)者，呵呵希望對(duì)你有幫助！ |
5樓2013-01-09 08:58:18
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +9 | 景隨風(fēng) 2026-03-29 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 343求調(diào)劑085601 +3 | 要努力學(xué)習(xí)x 2026-03-29 | 3/150 |
|
|
[考研] 070300化學(xué)354求調(diào)劑 +6 | 101次希望 2026-03-28 | 6/300 |
|
|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績(jī)338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么！ 2026-03-26 | 7/350 |
|
|
|
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
|
[考研]
 085602 289分求調(diào)劑
+8
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|