| 5 | 1/1 | 返回列表 |
| 查看: 2282 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
小劉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)行。 |
木蟲 (正式寫手)
榮譽(yù)版主 (著名寫手)
彩色掛圖
![]() |
專家經(jīng)驗(yàn): +2 |

木蟲 (正式寫手)
|
您可以幫我看看,改正點(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ā)表 | |
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +11 | 蛋黃咸肉粽 2026-03-26 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +6 | 小熊~~ 2026-03-31 | 6/300 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-31 | 4/200 |
|
|
[考研] 生物考研337分求調(diào)劑 +3 | cgxin 2026-03-30 | 5/250 |
|
|
[考研] 282求調(diào)劑 +3 | 呼吸都是減肥 2026-03-31 | 3/150 |
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +18 | 北道巷 2026-03-24 | 18/900 |
|
|
[考研] 284求調(diào)劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|