24小时热门版块排行榜

返回列表

zdzhang

木虫 (正式写手)

应助: 4 (幼儿园)
金币: 2685.9
散金: 668
红花: 3
帖子: 921
在线: 829.7小时
虫号: 1208276
注册: 2011-02-21
性别: GG
专业: 计算机体系结构

[求助] 求各位前辈帮忙！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

» 猜你喜欢

A期刊撤稿已经有3人回复
职称评审没过，求安慰已经有34人回复
垃圾破二本职称评审标准已经有17人回复
回收溶剂求助已经有6人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
申请26博士已经有5人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复

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

求助前辈帮忙看一下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人回复

1楼 2011-11-16 15:12:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wangyuewu08

木虫 (著名写手)

应助: 13 (小学生)
金币: 2794
红花: 4
帖子: 1141
在线: 279.1小时
虫号: 1290073
注册: 2011-05-09
性别: GG
专业: 制造系统与自动化

【答案】应助回帖

★
zdzhang(金币+2): 2011-12-06 09:32:03
zdzhang(金币+2): 2011-12-06 13:23:44
xiegangmai(金币+1): 谢谢应助！ 2011-12-06 22:18:21

错误出现在24行，我数了一下，好像是if f(k_z1)*f(k_z0)>0
这条语句出了问题，两个复数相乘得到的结果应该还是一个复数，除非它们是共轭的。当你把一个复数同0相比较时，程序会认为你是用实部与0比较，而忽略了虚步，警告就是这样显示的。
建议楼主求一下两复数乘积的模，然后再与0进行比较，希望对你有用。

赞一下(1人)

回复此楼

往事并不如烟

2楼2011-12-06 08:47:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wangyuewu08

木虫 (著名写手)

应助: 13 (小学生)
金币: 2794
红花: 4
帖子: 1141
在线: 279.1小时
虫号: 1290073
注册: 2011-05-09
性别: GG
专业: 制造系统与自动化

★ ★
xiegangmai(金币+2): 谢谢参与！ 2011-12-07 18:38:22

我把楼主的代码考到自己的电脑上，运行之后，发现了问题所在。
24行是输出plot(k_z0,f(k_z0),'+')的这条语句，在工作空间中发现k_z0和f(k_z0)都是复数，当x和y都是复数时，系统无法绘制相应的图形，所以提出警告。
我在command window中输入这样的语句plot( 3+4i, 4+3i )，系统显示相同的错误。
知道楼主的意图是想表示两个变量之间的关系，输入一个点，但如果在复数平面上表示的话，会输出两个点，无法反应自变量与应变量之间的关系。