声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1680|回复: 6

[编程技巧] 解一个复杂的的方程组

[复制链接]
发表于 2008-11-4 10:01 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
各位大侠,小弟刚学matlab,不会解这个方程,请大家看看,谢谢
求解一个未知数 kx

已知数:
a=1.51383^2;
b=-18.8997+i0.7937;
c=2.22;
d=1
d1=0.04;
d2=0.14;

(1)
kz0=sqrt(a*(2*pi/0.673)^2-kx^2);
kz1=sqrt(b*(2*pi/0.673)^2-kx^2);
kz2=sqrt(c*(2*pi/0.673)^2-kx^2);
kz3=sqrt(d*(2*pi/0.673)^2-kx^2);

(2)利用(1)的值
r01=(kz0/a-kz1/b)/(kz0/a+kz1/b);
r12=(kz1/b-kz2/c)/(kz1/b+kz2/c);
r23=(kz2/c-kz3/d)/(kz2/c+kz3/d);

(3)利用(2)的值
r123=[r12+r23*exp(i*2*kz2*d2)]/[1+r12*r23*exp(i*kz2*d2)]               %exp()中的首字母是i

(4)由上述所有从kx出发算的值可以得到下面一个等式
1+r01*r123*exp(i*2kz1*d1)=0


我不懂得如何利用上述式子求得kx的值,希望大家能给我解答一下,谢谢
回复
分享到:

使用道具 举报

发表于 2008-11-4 10:49 | 显示全部楼层
超过3次幂matlab很难解出来的,建议用mathematic试试。
发表于 2008-11-4 11:01 | 显示全部楼层

回复 楼主 zj_zhong 的帖子

help solve
 楼主| 发表于 2008-11-4 11:05 | 显示全部楼层

回复 板凳 sigma665 的帖子

我曾经在网上看到一些解方程组的例子,参照下用solve解,但是总是显示出错,也不知是语法错还是其他地方
发表于 2008-11-5 09:05 | 显示全部楼层
那就看错误信息是什么了
发表于 2008-11-5 09:37 | 显示全部楼层
不会啊























































                               
登录/注册后可看大图

绝世妖孽
宦海沉浮
发表于 2008-11-8 18:27 | 显示全部楼层
学习中,呵呵。莫怪。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-4 12:11 , Processed in 0.077988 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表