惜夢(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年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有251人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(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)有60人回復(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ù)值的話，改成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ā)表

[考研] 271材料工程求調(diào)劑 +6	.6lL 2026-03-18	6/300	2026-03-19 15:41 by haoshis
[考研] 274求調(diào)劑 +6	S.H1 2026-03-18	6/300	2026-03-19 09:34 by 花店相見(jiàn)
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 265求調(diào)劑 +3	梁梁校校 2026-03-17	3/150	2026-03-18 09:12 by zhukairuo
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專(zhuān)家開(kāi)始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥(niǎo)于飛Chao
[考研] 296求調(diào)劑 +5	大口吃飯身體健 2026-03-13	5/250	2026-03-17 21:05 by 不惑可樂(lè)
[考研] 268求調(diào)劑 +7	好運(yùn)連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu