| 查看: 2233 | 回复: 5 | ||
zdzhang木虫 (正式写手)
|
[求助]
求各位前辈帮忙！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» 猜你喜欢
一志愿085404，总分291，四级已过，求调剂
已经有9人回复
-
285求调剂
已经有6人回复
-
266求调剂
已经有8人回复
-
材料专硕调剂
已经有9人回复
-
化学工程调剂289
已经有21人回复
-
机械调剂
已经有3人回复
-
368化学求调剂
已经有4人回复
-
0703调剂，一志愿天津大学319分
已经有21人回复
-
0854电子信息319求调剂（接受跨专业调剂）
已经有6人回复
-
调剂
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助前辈帮忙看一下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人回复
wangyuewu08
木虫 (著名写手)
- 应助: 13 (小学生)
- 金币: 2794
- 红花: 4
- 帖子: 1141
- 在线: 279.1小时
- 虫号: 1290073
- 注册: 2011-05-09
- 性别: GG
- 专业: 制造系统与自动化
2楼2011-12-06 08:47:11
wangyuewu08
木虫 (著名写手)
- 应助: 13 (小学生)
- 金币: 2794
- 红花: 4
- 帖子: 1141
- 在线: 279.1小时
- 虫号: 1290073
- 注册: 2011-05-09
- 性别: GG
- 专业: 制造系统与自动化
★ ★
xiegangmai(金币+2): 谢谢参与！ 2011-12-07 18:38:22
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 )，系统显示相同的错误。 知道楼主的意图是想表示两个变量之间的关系，输入一个点，但如果在复数平面上表示的话，会输出两个点，无法反应自变量与应变量之间的关系。  上次回帖说错了，当复数和复数进行大小比较的时候，系统只比较实部，实数可以看做虚部为0的复数，所以比较语句并没有问题。 |
3楼2011-12-07 18:26:08
wangyuewu08
木虫 (著名写手)
- 应助: 13 (小学生)
- 金币: 2794
- 红花: 4
- 帖子: 1141
- 在线: 279.1小时
- 虫号: 1290073
- 注册: 2011-05-09
- 性别: GG
- 专业: 制造系统与自动化
4楼2011-12-07 18:38:37
zdzhang
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 2685.9
- 散金: 668
- 红花: 3
- 帖子: 921
- 在线: 829.7小时
- 虫号: 1208276
- 注册: 2011-02-21
- 性别: GG
- 专业: 计算机体系结构
5楼2011-12-08 09:33:19
6楼2014-10-31 22:48:14
