| 查看: 2213 | 回复: 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人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有5人回复
-
求个博导看看
已经有17人回复
-
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
自荐读博
已经有5人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 无法找到调试信息
已经有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
6楼2012-05-29 09:21:27
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