小劉8885

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

應(yīng)助: 3 (幼兒園)
金幣: 10042.8
散金: 200
帖子: 468
在線: 160.5小時(shí)
蟲(chóng)號(hào): 735175
注冊(cè): 2009-03-30
性別: MM
專(zhuān)業(yè): 材料物理化學(xué)

[求助] 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)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : func_main.m

2012-05-27 11:38:14, 945 bytes

附件 2 : 3.txt

2012-05-27 11:38:25, 3.1 K

附件 3 : Exop4.m

2012-05-27 11:39:11, 119 bytes

» 猜你喜歡

本人考085602 化學(xué)工程專(zhuān)碩已經(jīng)有23人回復(fù)
求調(diào)劑院校信息已經(jīng)有4人回復(fù)
085600材料與化工306 已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有10人回復(fù)
328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷已經(jīng)有9人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有4人回復(fù)
286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有8人回復(fù)
資源與環(huán)境調(diào)劑申請(qǐng)(333分) 已經(jīng)有5人回復(fù)
280求調(diào)劑已經(jīng)有12人回復(fù)
269專(zhuān)碩求調(diào)劑已經(jīng)有5人回復(fù)

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

無(wú)法找到調(diào)試信息已經(jīng)有3人回復(fù)
求助一下程序調(diào)試過(guò)程出現(xiàn)的錯(cuò)誤已經(jīng)有4人回復(fù)
螺旋波波頭程序調(diào)試除零錯(cuò)誤已經(jīng)有12人回復(fù)
matlab調(diào)試問(wèn)題已經(jīng)有3人回復(fù)
GDB單步調(diào)試問(wèn)題已經(jīng)有8人回復(fù)
對(duì)"推板式隧道電阻爐“調(diào)試精通的請(qǐng)進(jìn)��！已經(jīng)有31人回復(fù)
【求助】用C#編寫(xiě)的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)問(wèn)怎么用MATALAB編寫(xiě)一個(gè)空間幾何體的受力變形程序？？已經(jīng)有3人回復(fù)
【求助】KRUSS K10表面張力儀調(diào)試方法【有效期至2010年7月11日】已經(jīng)有3人回復(fù)
【討論】關(guān)于針對(duì)CASTEP收斂問(wèn)題的參數(shù)調(diào)試總結(jié)。已經(jīng)有79人回復(fù)

1樓 2012-05-27 11:41:35

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

csgt0

榮譽(yù)版主 (著名寫(xiě)手)

彩色掛圖

專(zhuān)家經(jīng)驗(yàn): +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲(chóng)號(hào): 449046
注冊(cè): 2007-11-02
專(zhuān)業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計(jì)算模擬

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

感謝參與，應(yīng)助指數(shù) +1

F里的x2^a應(yīng)該點(diǎn)乘吧？有什么問(wèn)題出什么錯(cuò)也不說(shuō)，有幾個(gè)能幫你

贊一下

回復(fù)此樓

showmethemoney

2樓2012-05-28 09:10:18

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

lijie169

銅蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 159 (高中生)
金幣: 161.5
散金: 2018
紅花: 7
帖子: 1338
在線: 491.5小時(shí)
蟲(chóng)號(hào): 1180408
注冊(cè): 2010-12-30
專(zhuān)業(yè): 分離過(guò)程

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

★
感謝參與，應(yīng)助指數(shù) +1
小劉8885: 金幣+1, ★有幫助 2012-05-29 10:09:49

plotyy(xf,a,xf,b);這句不對(duì)吧！�。�

贊一下

回復(fù)此樓

3樓2012-05-28 13:11:52

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

小劉8885

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

應(yīng)助: 3 (幼兒園)
金幣: 10042.8
散金: 200
帖子: 468
在線: 160.5小時(shí)
蟲(chóng)號(hào): 735175
注冊(cè): 2009-03-30
性別: MM
專(zhuān)業(yè): 材料物理化學(xué)

引用回帖:

2樓: Originally posted by csgt0 at 2012-05-28 09:10:18
F里的x2^a應(yīng)該點(diǎn)乘吧？有什么問(wèn)題出什么錯(cuò)也不說(shuō)，有幾個(gè)能幫你

您可以幫我看看，改正點(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);

贊一下

回復(fù)此樓

4樓2012-05-28 18:43:31

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

小劉8885

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

應(yīng)助: 3 (幼兒園)
金幣: 10042.8
散金: 200
帖子: 468
在線: 160.5小時(shí)
蟲(chóng)號(hào): 735175
注冊(cè): 2009-03-30
性別: MM
專(zhuān)業(yè): 材料物理化學(xué)

引用回帖:

3樓: Originally posted by lijie169 at 2012-05-28 13:11:52
plotyy(xf,a,xf,b);這句不對(duì)吧�。�！

這句語(yǔ)句就是畫(huà)一個(gè)雙y軸的的圖，好像沒(méi)有錯(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í)行的，還畫(huà)出了圖。運(yùn)行之后提示我問(wèn)題出在 fsolve 那個(gè)語(yǔ)句。

贊一下

回復(fù)此樓

5樓2012-05-28 18:46:10

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

csgt0

榮譽(yù)版主 (著名寫(xiě)手)

彩色掛圖

專(zhuān)家經(jīng)驗(yàn): +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲(chóng)號(hào): 449046
注冊(cè): 2007-11-02
專(zhuān)業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計(jì)算模擬

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小劉8885: 金幣+29, ★★★★★最佳答案, 果然是牛人，大才。 2012-05-29 10:09:20

x0=zeros(length(data),1);
改成x0=ones(length(data),1);
如果x0=0，y2=0,F2就無(wú)窮大了所以報(bào)錯(cuò)

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

小劉8885

showmethemoney

6樓2012-05-29 09:21:27

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

小劉8885

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

應(yīng)助: 3 (幼兒園)
金幣: 10042.8
散金: 200
帖子: 468
在線: 160.5小時(shí)
蟲(chóng)號(hào): 735175
注冊(cè): 2009-03-30
性別: MM
專(zhuān)業(yè): 材料物理化學(xué)

送鮮花一朵

引用回帖:

6樓: Originally posted by csgt0 at 2012-05-29 09:21:27
x0=zeros(length(data),1);
改成x0=ones(length(data),1);
如果x0=0，y2=0,F2就無(wú)窮大了所以報(bào)錯(cuò)

果然犀利，被你點(diǎn)中了。問(wèn)題解決了，謝謝！

贊一下

回復(fù)此樓

7樓2012-05-29 10:05:46

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

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

[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +5	張vvvv 2026-03-15	7/350	2026-03-21 19:23 by ColorlessPI
[考研] 313求調(diào)劑 +4	肆叁貳壹22 2026-03-19	4/200	2026-03-21 17:33 by ColorlessPI
[考研] 求調(diào)劑 +3	白QF 2026-03-21	3/150	2026-03-21 13:12 by zhukairuo
[考研] 316求調(diào)劑 +6	梁茜雯 2026-03-19	6/300	2026-03-21 06:32 by Ecowxq666！
[考研] 299求調(diào)劑 +6	△小透明* 2026-03-17	6/300	2026-03-21 02:42 by JourneyLucky
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3	innocenceF 2026-03-17	3/150	2026-03-21 02:35 by JourneyLucky
[考研] 一志愿中國(guó)石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) +3	石能偉 2026-03-17	3/150	2026-03-21 02:22 by JourneyLucky
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 一志愿西南交大，求調(diào)劑 +5	材化逐夢(mèng)人 2026-03-18	5/250	2026-03-21 00:26 by JourneyLucky
[考研] 材料專(zhuān)碩英一數(shù)二306 +7	z1z2z3879 2026-03-18	7/350	2026-03-20 23:48 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 350求調(diào)劑 +5	weudhdk 2026-03-19	5/250	2026-03-20 22:04 by luoyongfeng
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 085601專(zhuān)碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 一志愿蘇州大學(xué)材料工程（085601）專(zhuān)碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專(zhuān)利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111