24小时热门版块排行榜    

查看: 1871  |  回复: 5

wh1125

木虫 (正式写手)

[求助] 一个1stOpt的初级问题:复数方程组计算问题(有代码)

初学1stOpt,以下1stOpt代码中有5个方程,求解5个变量K00,K01,K02,K03,beta0,其中i已定义为虚数单位,但是运行结果中,怎么总把i当成一个变量,有结果i=0。怎么样才能正确计算我要求解的那5个变量呢?请大神指点啊!

ComplexStr = i;
Constant K=2*pi/(0.88e-6),n0=3.4,n1=3.6,n2=3.4,n3=3.55,t=1e-6,g=0.2e-6;
Function
K00^2=K^2*n0^2-beta0^2;
K01^2=K^2*n1^2-beta0^2;
K02^2=K^2*n2^2-beta0^2;
K03^2=K^2*n3^2-beta0^2;
K03*(K01^2+K00*K02)*tan(K01*(t-g))+K01*(K03^2+K00*K02)*tan(K03*g)+i*(K01*K03*(K00+K02)-(K00*K03^2+K01^2*K02)*tan(K03*g)*tan(K01*(t-g)))=0;

运行结果:
k00: 3.28705929208081
beta0: 4.83100198907778
k01: 2.24062444176525
k02: 4.41804389934987
k03: 3.40456867241301
i: 0
回复此楼

» 猜你喜欢

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

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

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wh1125: 金币+5, 有帮助, 谢谢你!按照你的办法能够求解,但这样得出的是实数解,没有复数。结果不是我期望的。 2013-10-07 19:45:33
dingd: 金币+1, 谢谢交流。 2013-10-07 21:45:05
删除ComplexStr = i;
拆分
K03*(K01^2+K00*K02)*tan(K01*(t-g))+K01*(K03^2+K00*K02)*tan(K03*g)+i*(K01*K03*(K00+K02)-(K00*K03^2+K01^2*K02)*tan(K03*g)*tan(K01*(t-g)))=0;
成两个式子:
K03*(K01^2+K00*K02)*tan(K01*(t-g))+K01*(K03^2+K00*K02)*tan(K03*g)=0;
(K01*K03*(K00+K02)-(K00*K03^2+K01^2*K02)*tan(K03*g)*tan(K01*(t-g)))=0;
因为负数为零等价于实部为0 & 虚部为0
2楼2013-10-07 19:01:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wh1125

木虫 (正式写手)

引用回帖:
2楼: Originally posted by feixiaolin at 2013-10-07 19:01:58
删除ComplexStr = i;
拆分
K03*(K01^2+K00*K02)*tan(K01*(t-g))+K01*(K03^2+K00*K02)*tan(K03*g)+i*(K01*K03*(K00+K02)-(K00*K03^2+K01^2*K02)*tan(K03*g)*tan(K01*(t-g)))=0;
成两个式子:
K03*(K01^2+K00*K0 ...

虽然能得出结果,全部是实数解,但是明显解的数量级都不对。
3楼2013-10-07 19:47:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
csgt0: 金币+1, 谢谢 2013-10-09 10:30:26
不知道一楼的问题中待求参数是实数还是虚数?如果都是实数,2楼的方法很好啊,只是似乎无解。方程是否正确?
4楼2013-10-07 21:49:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wh1125

木虫 (正式写手)

引用回帖:
4楼: Originally posted by dingd at 2013-10-07 21:49:52
不知道一楼的问题中待求参数是实数还是虚数?如果都是实数,2楼的方法很好啊,只是似乎无解。方程是否正确?

方程是正确的。解没有限制,可以是实数解也可以是复数解。
5楼2013-10-07 22:32:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★
wh1125: 金币+5, 有帮助, 谢谢!我的解要求K01和K03不为零,因为后面还有计算,他们要做分母。 2013-10-07 23:13:11
csgt0: 金币+1, 谢谢 2013-10-09 10:30:42
下面这组参考下:

k00.实部: 0
k00.虚部: 7290129.24733499
k01.实部: 4269089.10229991
k01.虚部: 0
k02.实部: 0
k02.虚部: -7290129.24733499
k03.实部: 0
k03.虚部: 0
beta0.实部: 25346940.7278267
beta0.虚部: 0
6楼2013-10-07 23:09:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wh1125 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085700资源与环境308求调剂 +6 墨墨漠 2026-03-18 6/300 2026-03-20 00:02 by 23Postgrad
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 梁成伟老师课题组欢迎你的加入 +9 一鸭鸭哟 2026-03-14 11/550 2026-03-19 17:22 by !本暗一次!
[考研] 一志愿南昌大学,327分,材料与化工085600 +3 Ncdx123456 2026-03-19 3/150 2026-03-19 13:18 by houyaoxu
[考研] 一志愿中海洋材料工程专硕330分求调剂 +7 小材化本科 2026-03-18 7/350 2026-03-19 10:46 by Linda Hu
[考研] 材料专硕英一数二306 +5 z1z2z3879 2026-03-18 5/250 2026-03-19 07:43 by BruceLiu320
[考研] 304求调剂 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 354求调剂 +4 Tyoumou 2026-03-18 7/350 2026-03-18 21:45 by Tyoumou
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +7 步川酷紫123 2026-03-13 7/350 2026-03-18 17:12 by 尽舜尧1
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
信息提示
请填处理意见