 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 2272 | 回復(fù): 6 | |||
小劉8885木蟲 (正式寫手)
|
[求助]
Matalab 程序調(diào)試
|
|
下述程序運(yùn)行出錯(cuò)，求高手指點(diǎn)。這是一個(gè)M文件。 function f4=func_main(name,a,q1,q2,x2,flag) % flag=0 represent complex permittivity % flag=1 represent complex permeability data=load(name); xf=data(:,1); if (flag==0) y1=data(:,2)+data(:,3)*i; elseif (flag==1) y1=data(:,4)+data(:,5)*i; else disp('0 is complex perimittivity, 1 is complex permeability'); return; end x11=func1(a,q1,y1,x2); %調(diào)用函數(shù)func1，求解x1 %Hanai formula 的求解結(jié)果 %求解的是x1，為絕對(duì)介電常數(shù) x1=x11; a=real(x1); b=imag(x1); c=[a,b]; x0=zeros(length(data),1); F=fsolve(@(y2)func2(a,q2,x1,x2,y2),x0); %調(diào)用前述func1所求得的x1，代入函數(shù)func2 %求解超越方程的y2值。 f4=F; r = [real(f4), imag(f4)]; plotyy(xf,a,xf,b); save('data1.txt', 'c', '-ascii'); save('data2.txt', 'r', '-ascii'); function x1=func1(a,q1,y1,x2) x1=(q1.*x2-x2+y1.*(x2./y1).^a)./(q1+(x2./y1).^a-1); %求解x1 function F=func2(a,q2,x1,x2,y2) F=(((1-q2)*(x2-x1))./(x2^a)).*(y2.^a)-y2+x1; %求解y2 另外的調(diào)用代碼如下： % Exop4 clear up; clc; a=1/3; q1=0.06; q2=0.107; x2=3.7+0.703*i; flag=0; func_main('3.txt',a,q1,q2,x2,flag); 復(fù)制到command窗口就可以運(yùn)行。 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : func_main.m
- 附件 2 : 3.txt
- 附件 3 : Exop4.m
2012-05-27 11:38:14, 945 bytes
2012-05-27 11:38:25, 3.1 K
2012-05-27 11:39:11, 119 bytes
» 猜你喜歡
0703化學(xué) 305求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿211，0703化學(xué)310分求調(diào)劑
已經(jīng)有3人回復(fù)
-
廣西大學(xué)材料導(dǎo)師推薦
已經(jīng)有5人回復(fù)
-
化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
資源與環(huán)境 調(diào)劑申請(qǐng)(333分)
已經(jīng)有4人回復(fù)
-
311求調(diào)劑
已經(jīng)有14人回復(fù)
-
269專碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
085600材料與化工306
已經(jīng)有3人回復(fù)
-
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑
已經(jīng)有12人回復(fù)
-
材料求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 無(wú)法找到調(diào)試信息
已經(jīng)有3人回復(fù)
- 求助一下程序調(diào)試過程出現(xiàn)的錯(cuò)誤
已經(jīng)有4人回復(fù)
- 螺旋波波頭 程序調(diào)試 除零錯(cuò)誤
已經(jīng)有12人回復(fù)
- matlab調(diào)試問題
已經(jīng)有3人回復(fù)
- GDB單步調(diào)試問題
已經(jīng)有8人回復(fù)
- 對(duì)"推板式隧道電阻爐“調(diào)試精通的請(qǐng)進(jìn)！！
已經(jīng)有31人回復(fù)
- 【求助】用C#編寫的Windows應(yīng)用程序中，用到了數(shù)組，不能成功調(diào)試
已經(jīng)有4人回復(fù)
- 【求助】厭氧調(diào)試？
已經(jīng)有8人回復(fù)
- 【求助】VS2005數(shù)值計(jì)算程序調(diào)試【急求�。。�！加金幣】
已經(jīng)有9人回復(fù)
- 【求助】VC++數(shù)值算法集應(yīng)用到MFC中報(bào)錯(cuò)，程序調(diào)試求助
已經(jīng)有8人回復(fù)
- 【求助】VC++程序算法集調(diào)試報(bào)錯(cuò)
已經(jīng)有8人回復(fù)
- 【求助】請(qǐng)問怎么用MATALAB編寫一個(gè)空間幾何體的受力變形程序？？
已經(jīng)有3人回復(fù)
- 【求助】KRUSS K10表面張力儀調(diào)試方法【有效期至2010年7月11日】
已經(jīng)有3人回復(fù)
- 【討論】關(guān)于針對(duì)CASTEP收斂問題的參數(shù)調(diào)試總結(jié)。
已經(jīng)有79人回復(fù)
csgt0
榮譽(yù)版主 (著名寫手)
彩色掛圖
-
專家經(jīng)驗(yàn): +2 - 應(yīng)助: 367 (碩士)
- 貴賓: 0.457
- 金幣: 7807
- 散金: 8897
- 紅花: 20
- 沙發(fā): 7
- 帖子: 2809
- 在線: 515.8小時(shí)
- 蟲號(hào): 449046
- 注冊(cè): 2007-11-02
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
- 管轄: 計(jì)算模擬
2樓2012-05-28 09:10:18
3樓2012-05-28 13:11:52
小劉8885
木蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 10042.8
- 散金: 200
- 帖子: 468
- 在線: 160.5小時(shí)
- 蟲號(hào): 735175
- 注冊(cè): 2009-03-30
- 性別: MM
- 專業(yè): 材料物理化學(xué)
|
您可以幫我看看，改正點(diǎn)乘之后，運(yùn)行之后顯示下面的錯(cuò)誤： Error using trustnleqn (line 28) Objective function is returning undefined values at initial point. FSOLVE cannot continue. Error in fsolve (line 399) [x,FVAL,JACOB,EXITFLAG,OUTPUT,msgData]=... Error in func_main (line 27) F=fsolve(@(y2)func2(a,q2,x1,x2,y2),x0); |
4樓2012-05-28 18:43:31
小劉8885
木蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 10042.8
- 散金: 200
- 帖子: 468
- 在線: 160.5小時(shí)
- 蟲號(hào): 735175
- 注冊(cè): 2009-03-30
- 性別: MM
- 專業(yè): 材料物理化學(xué)
|
這句語(yǔ)句就是畫一個(gè)雙y軸的的圖，好像沒有錯(cuò)啊。 function f4=func_main(name,a,q1,q2,X,flag) % flag=0 represent complex permittivity % flag=1 represent complex permeability data=load(name); xf=data(:,1); x2=X*ones(length(data),1); if (flag==0) y1=data(:,2)+data(:,3)*i; elseif (flag==1) y1=data(:,4)+data(:,5)*i; else disp('0 is complex perimittivity, 1 is complex permeability'); return; end x11=func1(a,q1,y1,x2); %調(diào)用函數(shù)func1，求解x1 %Hanai formula 的求解結(jié)果 %求解的是x1，為絕對(duì)介電常數(shù) x1=x11; a=real(x1); b=imag(x1); c=[a,b]; plotyy(xf,a,xf,b); x0=zeros(length(data),1); F=fsolve(@(y2)func2(a,q2,x1,x2,y2),x0); %調(diào)用前述func1所求得的x1，代入函數(shù)func2 %求解超越方程的y2值。 f4=F; r = [real(f4), imag(f4)]; save('data1.txt', 'c', '-ascii'); save('data2.txt', 'r', '-ascii'); function x1=func1(a,q1,y1,x2) x1=(q1.*x2-x2+y1.*(x2./y1).^a)./(q1+(x2./y1).^a-1); %求解x1 function F=func2(a,q2,x1,x2,y2) F=(((1-q2)*(x2-x1))./(x2.^a)).*(y2.^a)-y2+x1; %求解y2 我把 plotyy(xf,a,xf,b); 這一語(yǔ)句提到前面去是可以執(zhí)行的，還畫出了圖。運(yùn)行之后提示我問題出在 fsolve 那個(gè)語(yǔ)句。 |
5樓2012-05-28 18:46:10
csgt0
榮譽(yù)版主 (著名寫手)
彩色掛圖
-
專家經(jīng)驗(yàn): +2 - 應(yīng)助: 367 (碩士)
- 貴賓: 0.457
- 金幣: 7807
- 散金: 8897
- 紅花: 20
- 沙發(fā): 7
- 帖子: 2809
- 在線: 515.8小時(shí)
- 蟲號(hào): 449046
- 注冊(cè): 2007-11-02
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
- 管轄: 計(jì)算模擬
小劉88856樓2012-05-29 09:21:27
小劉8885
木蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 10042.8
- 散金: 200
- 帖子: 468
- 在線: 160.5小時(shí)
- 蟲號(hào): 735175
- 注冊(cè): 2009-03-30
- 性別: MM
- 專業(yè): 材料物理化學(xué)
7樓2012-05-29 10:05:46
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 小劉8885 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué) 305求調(diào)劑 +6 | FY_yy 2026-03-14 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +13 | 冬十三 2026-03-15 | 14/700 |
|
|
|
[考研] 269專碩求調(diào)劑 +4 | 金恩貝 2026-03-21 | 4/200 |
|
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
|
[考研] 326求調(diào)劑 +4 | mlpqaz03 2026-03-15 | 4/200 |
|
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +20 | 不知道叫什么！ 2026-03-15 | 22/1100 |
|
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5 | 守候夕陽(yáng)CF 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無(wú)極 2026-03-20 | 5/250 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|