| 查看: 2556 | 回复: 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
» 猜你喜欢
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
言之凿凿，8月21号(明天)放榜
已经有14人回复
-
时间戳变了，能看出什么问题？
已经有17人回复
-
只有每年这种时候来逛逛小木虫
已经有21人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有125人回复
-
今天基金会出结果吗？20260819
已经有15人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 无法找到调试信息
已经有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人回复
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
- 专业: 材料物理化学
4楼2012-05-28 18:43:31
小刘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
小刘88856楼2012-05-29 09:21:27
小刘8885
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 10042.8
- 散金: 200
- 帖子: 468
- 在线: 160.5小时
- 虫号: 735175
- 注册: 2009-03-30
- 性别: MM
- 专业: 材料物理化学
7楼2012-05-29 10:05:46
