| 查看: 2422 | 回复: 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 |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有214人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 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
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3916
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24222
- 在线: 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
- 帖子: 499
- 在线: 432.2小时
- 虫号: 417863
- 注册: 2007-07-02
- 性别: GG
- 专业: 集成电路设计与测试
3楼2013-10-07 19:47:58
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2013-10-07 21:49:52
