24小时热门版块排行榜    

查看: 1870  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 0703化学调剂 +8 妮妮ninicgb 2026-03-15 12/600 2026-03-19 05:55 by anny19840123
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
信息提示
请填处理意见