[求助] 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊

各位前辈好，我想通过解这样一个复数方程，得到k_z的值，但是我这样编写的程序，总是出错，哪位大侠能帮我看看是那个地方出错了呢？望帮忙修改一下，先谢过各位了！ 下面是我自己写的程序： clc; clear; syms k_z i=sqrt(-1); a=60*10^(-9); mu_air=4*pi*10^(-7); epsl_air=8.854187817e-12; epsl_met=-18.36+0.46i; c=3.0*10^8; wavelength=632.8*10^(-9); k_0=2*pi/wavelength; f=@(k_z) (-epsl_met*k_0^2*(8-(epsl_met*k_0^2-k_z^2)^(0.5)*a))/16*(1-0.25*(epsl_met*k_0^2-k_z^2)*a^2)+mu_air*epsl_air*k_0^2/((mu_air*epsl_air*k_0^2-k_z^2)*a^2*log(1.1229/(mu_air*epsl_air*k_0^2-k_z^2)^(0.5)*a)) ferror=1e-5; k_z1=(0.95+0.0008*i)*k_0; k_z2=(1.2+0.005*i)*k_0; fplot(f,[k_z1,k_z2]) xlabel('k_z') ylabel('f(k_z)') hold on %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% k_z0=(k_z1+k_z2)/2; while abs(f(k_z0))>ferror     plot(k_z0,f(k_z0),'+')     if f(k_z1)*f(k_z0)>0         k_z1=k_z0;     else         k_z2=k_z0;     end     k_z0=(k_z1+k_z2)/2; end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% format short disp('[           k_z0            f(k_z0)') disp([k_z0 f(k_z0)]) 那位大侠能帮忙修改一下我的程序呢？小弟在此谢过各位了 [ 来自科研家族 Ocean家族 ] 这是我要解的复数方程及已知参数 这是出现的警告 [ Last edited by zdzhang on 2011-11-16 at 17:29 ]

回复此楼

» 本帖已获得的红花（最新10朵）

wxliu1992

» 猜你喜欢

• 论文终于录用啦！满足毕业条件了 已经有10人回复
• 2025年遐想 已经有4人回复
• 投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球 已经有8人回复
• 求个博导看看 已经有18人回复

» 本主题相关价值贴推荐，对您同样有帮助:

• 求助前辈帮忙看一下AFM图 已经有3人回复
• 物理模型中的复数解 已经有12人回复
• 哪位前辈能够帮小弟一把，小弟在此谢过了。 已经有2人回复
• matlab解一元三次方程，得到的都是复数根。 已经有7人回复
• matlab fsolve得出的解全部为复数解 已经有15人回复
• 咨询下各位前辈，做水平井物理模拟装置的问题 已经有4人回复
• 这句要不要简写一下，谢谢前辈~~ 已经有3人回复
• 10次代数方程怎么找解析解 最好详细给出matlab程序 谢谢 已经有5人回复
• 【讲义】MATLAB课件（东南大学） 已经有154人回复
• 各位朋友，如何用matlab编写这个方程，并求解？ 已经有15人回复
• 拜求，这种带参数的方程怎么解比较好？ 已经有4人回复
• 求matlab 解方程的问题 已经有5人回复
• 帮忙既精通数学又精通VB编程的大神解指数方程 已经有4人回复
• 已知一随机信号的功率谱密度函数，如何利用Matlab产生其单次时域样本 已经有6人回复
• 求回归方程参数估计的最优化算法matlab代码 已经有9人回复
• 请各位前辈帮我看下这句话是什么意思，谢谢了先 已经有3人回复
• 请问高压下，模拟水的物性方程用哪个？谢谢大侠们 已经有3人回复
• MATLAB解方程出现??? Subscripted assignment dimension mismatch. 已经有3人回复
• matlab求积分，出现奇点怎么办？ 已经有8人回复
• matlab的command window的问题 已经有3人回复
• 【求助】matlab 解矩阵方程 重金悬赏解决方案 已经有10人回复
• 【求助】方程组有解该用单数还是用复数 已经有4人回复