版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

小劉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)容由用戶(hù)自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶(hù)本人，如對(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

» 猜你喜歡

085601材料工程找調(diào)劑已經(jīng)有11人回復(fù)
總分293求調(diào)劑已經(jīng)有9人回復(fù)
334分一志愿武理材料求調(diào)劑已經(jīng)有5人回復(fù)
327求調(diào)劑已經(jīng)有6人回復(fù)
329求調(diào)劑已經(jīng)有10人回復(fù)
275求調(diào)劑已經(jīng)有15人回復(fù)
291求調(diào)劑已經(jīng)有6人回復(fù)
309求調(diào)劑已經(jīng)有7人回復(fù)
297求調(diào)劑已經(jīng)有11人回復(fù)
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有6人回復(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的回帖

小劉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的回帖

查看全部 7 個(gè)回答

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é)

引用回帖:

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的回帖

查看全部 7 個(gè)回答

普通表情龍兔虎貓

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

[考研] 297求調(diào)劑 +11	田洪有 2026-03-26	11/550	2026-03-29 13:14 by peike
[考研] 352分-085602-一志愿985 +5	海納百川Ly 2026-03-29	5/250	2026-03-29 09:57 by Sjndkwm
[考研] 0856求調(diào)劑 +13	zhn03 2026-03-25	14/700	2026-03-29 08:13 by fmesaito
[考研] 316求調(diào)劑 +7	江辭666 2026-03-26	7/350	2026-03-28 21:28 by sanrepian
[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +3	AZMK 2026-03-27	5/250	2026-03-28 16:19 by xxxsssccc
[考研] 312，生物學(xué)求調(diào)劑 +3	小譯同學(xué)abc 2026-03-28	3/150	2026-03-28 15:32 by 落�？伤�
[考研] 328求調(diào)劑 +7	嗯滴的基本都 2026-03-27	7/350	2026-03-28 04:19 by fmesaito
[考研] 340求調(diào)劑 +5	jhx777 2026-03-27	5/250	2026-03-28 04:18 by fmesaito
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4	Ripcord上岸 2026-03-27	4/200	2026-03-27 17:00 by zhyzzh
[考研] 287求調(diào)劑 +10	land xuxu 2026-03-26	10/500	2026-03-27 15:33 by 帕爾馬拉特
[考研] 085601 材料工程 313分求調(diào)劑 +5	Ong3 2026-03-27	5/250	2026-03-27 12:24 by goldfish51
[考研] 333求調(diào)劑 +3	question挽風(fēng) 2026-03-23	3/150	2026-03-27 11:29 by 不吃魚(yú)的貓
[考研] 286求調(diào)劑 +4	lim0922 2026-03-26	4/200	2026-03-27 10:28 by guoweigw
[考研] 332求調(diào)劑 +6	032500 2026-03-25	6/300	2026-03-25 22:45 by 418490947
[考研] 290分調(diào)劑求助 +3	吉祥止止陳 2026-03-25	3/150	2026-03-25 19:58 by barlinike
[考研] B區(qū)考研調(diào)劑 +4	yqdszhdap－ 2026-03-22	5/250	2026-03-25 08:51 by baoball
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7	小羊睡著了? 2026-03-23	10/500	2026-03-25 02:24 by 清風(fēng)拂揚(yáng)。 m
[考研] 材料考研調(diào)劑生 +3	黃粱一夢(mèng)千年 2026-03-24	3/150	2026-03-24 17:00 by barlinike
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋?xiě)?/a>