惜夢(mèng)尋草

應(yīng)助: 0 (幼兒園)
金幣: 875.5
帖子: 25
在線: 5.7小時(shí)
蟲(chóng)號(hào): 2118505
注冊(cè): 2012-11-10
性別: MM
專(zhuān)業(yè): 光學(xué)

[求助] 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年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有182人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有30人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有62人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

matlab程序出錯(cuò)調(diào)試已經(jīng)有7人回復(fù)
搞地學(xué)的朋友，請(qǐng)問(wèn)怎么用matlab做Mann-Kendall趨勢(shì)檢驗(yàn)��？已經(jīng)有10人回復(fù)
matlab調(diào)試問(wèn)題已經(jīng)有3人回復(fù)
求幫忙啊~ matlab一多項(xiàng)式結(jié)果太長(zhǎng)，怎么寫(xiě)到txt里？已經(jīng)有14人回復(fù)
matlab調(diào)試問(wèn)題已經(jīng)有3人回復(fù)
while語(yǔ)句死循環(huán)了，求教已經(jīng)有6人回復(fù)
求助matlab編程已經(jīng)有13人回復(fù)
【求助】VS2005數(shù)值計(jì)算程序調(diào)試【急求�。。�！加金幣】已經(jīng)有9人回復(fù)

1樓 2013-01-07 22:50:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shutaham

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 16 (小學(xué)生)
金幣: 788.5
散金: 50
紅花: 2
帖子: 634
在線: 413.7小時(shí)
蟲(chóng)號(hào): 707269
注冊(cè): 2009-02-24
專(zhuān)業(yè): 光學(xué)

【答案】應(yīng)助回帖

★ ★ ★
感謝參與，應(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ì)問(wèn)題出在那個(gè) 1./A 上

贊一下

回復(fù)此樓

2樓2013-01-08 10:10:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jiahongen

新蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 180.9
散金: 9
紅花: 1
帖子: 48
在線: 26.7小時(shí)
蟲(chóng)號(hào): 633313
注冊(cè): 2008-10-22
專(zhuān)業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

呵呵，正在學(xué)習(xí)中

回復(fù)此樓

3樓2013-01-08 13:46:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

惜夢(mèng)尋草

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 875.5
帖子: 25
在線: 5.7小時(shí)
蟲(chóng)號(hào): 2118505
注冊(cè): 2012-11-10
性別: MM
專(zhuān)業(yè): 光學(xué)

引用回帖:

2樓: Originally posted by shutaham at 2013-01-08 10:10:51
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);
>&g ...

謝謝！那應(yīng)該怎么改呢？A應(yīng)該不是一個(gè)值

贊一下

回復(fù)此樓

4樓2013-01-08 19:49:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lgycjpcqu

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 35 (小學(xué)生)
金幣: 566.8
散金: 411
帖子: 371
在線: 176.9小時(shí)
蟲(chóng)號(hào): 1327328
注冊(cè): 2011-06-20
性別: GG
專(zhuān)業(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

A向量中共有401個(gè)元素，m的取值范圍由于判定語(yǔ)句
if lamdai(n)>lamdas(n);
的存在必然會(huì)小于401（如果等于此判定就沒(méi)有意義）
所以錯(cuò)誤提示說(shuō)，w和A中的元素要求相等。
如果語(yǔ)句
w(m)=abs(np(m)./lamdap-ns(m)./lamdass(m)-ni(m)./lamdaii(m)-1./A);
的意思說(shuō)當(dāng)上面的判定成立時(shí)，用A向量對(duì)應(yīng)位置的元素計(jì)算w的數(shù)值的話(huà)，改成1/A（m）即可。
我也是初學(xué)者，呵呵希望對(duì)你有幫助！

贊一下

回復(fù)此樓

5樓2013-01-09 08:58:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

惜夢(mèng)尋草

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 875.5
帖子: 25
在線: 5.7小時(shí)
蟲(chóng)號(hào): 2118505
注冊(cè): 2012-11-10
性別: MM
專(zhuān)業(yè): 光學(xué)

引用回帖:

5樓: Originally posted by lgycjpcqu at 2013-01-09 08:58:18
A向量中共有401個(gè)元素，m的取值范圍由于判定語(yǔ)句
if lamdai(n)>lamdas(n);
的存在必然會(huì)小于401（如果等于此判定就沒(méi)有意義）
所以錯(cuò)誤提示說(shuō)，w和A中的元素要求相等。
如果語(yǔ)句
w(m)=abs(np(m)./lamdap- ...

修改了以后運(yùn)行一下還是有問(wèn)題
??? Undefined function or variable 'A1'.
不知道怎么改了

贊一下

回復(fù)此樓

6樓2013-01-09 18:57:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

connors

禁蟲(chóng) (小有名氣)

★
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+1, 謝謝 2013-01-10 10:35:48

本帖內(nèi)容被屏蔽

7樓2013-01-09 20:15:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lgycjpcqu

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 35 (小學(xué)生)
金幣: 566.8
散金: 411
帖子: 371
在線: 176.9小時(shí)
蟲(chóng)號(hào): 1327328
注冊(cè): 2011-06-20
性別: GG
專(zhuān)業(yè): 電化學(xué)

【答案】應(yīng)助回帖

★ ★
csgt0: 金幣+1, 謝謝 2013-01-10 10:35:54
惜夢(mèng)尋草: 金幣+1, ★有幫助, 謝謝啦！ 2013-01-11 08:16:16

引用回帖:

6樓: Originally posted by 惜夢(mèng)尋草 at 2013-01-09 18:57:04
修改了以后運(yùn)行一下還是有問(wèn)題
??? Undefined function or variable 'A1'.
不知道怎么改了

...

我看了一下，修改后
if min(w)<10^-5
的判斷不能成立，所以A1沒(méi)有數(shù)值，
按照你的程序j計(jì)算后好像
min(w)=3.054*10^(-5)
Undefined function or variable 'A1'.
就是A1沒(méi)有進(jìn)行定義或是你就算中沒(méi)有賦值，好好檢查一下計(jì)算A1的計(jì)算程序即可發(fā)現(xiàn)問(wèn)題。

贊一下

回復(fù)此樓

8樓2013-01-09 22:20:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 288求調(diào)劑 +14	于海海海海 2026-03-19	14/700	2026-03-19 19:40 by maocaozhuxi
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	10/500	2026-03-19 16:05 by 余麟余
[考研] 一志愿武漢理工材料工程專(zhuān)碩調(diào)劑 +4	Doleres 2026-03-19	4/200	2026-03-19 15:27 by 丁丁*
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專(zhuān)利，省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開(kāi)心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 0856專(zhuān)碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專(zhuān)家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086