| 查看: 1772 | 回复: 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 |
回复此楼» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有284人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stopt能进行隐函数非线性拟合吗？谁能帮拟合一下新型液体的粘度混合计算公式
已经有5人回复
- 求大神用1stOpt帮忙模拟一个复数拟合。
已经有3人回复
- 急求，常微分方程的参数计算问题，有1stopt的帮忙解答一下
已经有5人回复
- 求大神帮忙用1stOpt 拟合复数数据
已经有13人回复
- 求助，用1stopt解方程组，两次得到的解不唯一？什么原因，大神帮忙
已经有6人回复
- sos!!1stopt解非线性方程组，一直出现常数定义错误的问题！求大神帮助！
已经有12人回复
- 求大侠帮助跑一下1stopt
已经有12人回复
- 哪位朋友有注册版1stOpt软件？急需
已经有5人回复
- 用1stopt拟合非线性方程结果与其他软件拟合结果差异大
已经有5人回复
- 求高版本1stopt，解决一个问题，微分方程组参数拟合
已经有3人回复
- 1stopt运行程序，总是出现字符函数错误的提示，请高手帮助
已经有5人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 1stopt4阶非线性常微分方程，帮小弟看一下，跪求。
已经有7人回复
- 1stOpt高版本求助！求高版本的帮忙走一下 谢了！
已经有16人回复
- 急切求助！一个mathematica求解一个简单一元非线性方程的问题
已经有4人回复
- 1stOpt求解常微分方程边值问题
已经有5人回复
- 哪位大哥能帮我用1stopt3.0以上版本跑一个拟合程序
已经有11人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
wh1125: 金币+5, ★有帮助, 谢谢你！按照你的办法能够求解，但这样得出的是实数解，没有复数。结果不是我期望的。 2013-10-07 19:45:33
dingd: 金币+1, 谢谢交流。 2013-10-07 21:45:05
感谢参与，应助指数 +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
wh1125
木虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 3078.9
- 红花: 5
- 帖子: 498
- 在线: 431.2小时
- 虫号: 417863
- 注册: 2007-07-02
- 性别: GG
- 专业: 集成电路设计与测试
3楼2013-10-07 19:47:58
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2013-10-07 21:49:52
wh1125
木虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 3078.9
- 红花: 5
- 帖子: 498
- 在线: 431.2小时
- 虫号: 417863
- 注册: 2007-07-02
- 性别: GG
- 专业: 集成电路设计与测试
5楼2013-10-07 22:32:19
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★
wh1125: 金币+5, ★有帮助, 谢谢！我的解要求K01和K03不为零，因为后面还有计算，他们要做分母。 2013-10-07 23:13:11
csgt0: 金币+1, 谢谢 2013-10-09 10:30:42
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