| 查看: 2266 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
小刘8885木虫 (正式写手)
|
[求助]
Matalab 程序调试
|
|
|
下述程序运行出错，求高手指点。这是一个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); %调用函数func1，求解x1 %Hanai formula 的求解结果 %求解的是x1，为绝对介电常数 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); %调用前述func1所求得的x1，代入函数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 另外的调用代码如下： % 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); 复制到command窗口就可以运行。 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：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
» 猜你喜欢
279分求调剂 一志愿211
已经有12人回复
-
297求调剂
已经有9人回复
-
085700资源与环境308求调剂
已经有13人回复
-
材料 336 求调剂
已经有4人回复
-
华东师范大学-071000生物学-293分-求调剂
已经有3人回复
-
280求调剂
已经有8人回复
-
324分 085600材料化工求调剂
已经有4人回复
-
一志愿华中科技大学，080502，354分求调剂
已经有5人回复
-
一志愿武理材料305分求调剂
已经有7人回复
-
271材料工程求调剂
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 无法找到调试信息
已经有3人回复
- 求助一下程序调试过程出现的错误
已经有4人回复
- 螺旋波波头 程序调试 除零错误
已经有12人回复
- matlab调试问题
已经有3人回复
- GDB单步调试问题
已经有8人回复
- 对"推板式隧道电阻炉“调试精通的请进！！
已经有31人回复
- 【求助】用C#编写的Windows应用程序中，用到了数组，不能成功调试
已经有4人回复
- 【求助】厌氧调试？
已经有8人回复
- 【求助】VS2005数值计算程序调试【急求！！！！加金币】
已经有9人回复
- 【求助】VC++数值算法集应用到MFC中报错，程序调试求助
已经有8人回复
- 【求助】VC++程序算法集调试报错
已经有8人回复
- 【求助】请问怎么用MATALAB编写一个空间几何体的受力变形程序？？
已经有3人回复
- 【求助】KRUSS K10表面张力仪调试方法【有效期至2010年7月11日】
已经有3人回复
- 【讨论】关于针对CASTEP收敛问题的参数调试总结。
已经有79人回复
小刘8885
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 10042.8
- 散金: 200
- 帖子: 468
- 在线: 160.5小时
- 虫号: 735175
- 注册: 2009-03-30
- 性别: MM
- 专业: 材料物理化学
4楼2012-05-28 18:43:31
2楼2012-05-28 09:10:18
3楼2012-05-28 13:11:52
小刘8885
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 10042.8
- 散金: 200
- 帖子: 468
- 在线: 160.5小时
- 虫号: 735175
- 注册: 2009-03-30
- 性别: MM
- 专业: 材料物理化学
|
这句语句就是画一个双y轴的的图，好像没有错啊。 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); %调用函数func1，求解x1 %Hanai formula 的求解结果 %求解的是x1，为绝对介电常数 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); %调用前述func1所求得的x1，代入函数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); 这一语句提到前面去是可以执行的，还画出了图。运行之后提示我问题出在 fsolve 那个语句。 |
5楼2012-05-28 18:46:10
10