| 查看: 631 | 回复: 0 | ||
wh1125木虫 (正式写手)
|
[求助]
求助：matlab求解一个方程组，最后一个方程有三角函数和虚数单位i
|
|
代码如下，求解：K00，K01，K02，K03，beta0五个变量。可能我的x0初值猜的不合适吧，总不能得到解。（说明：解中K01和K03均不能等于0） 求大神看下如何能正确求解！ function Cal_K0i_beta0 clear;clc; global k0 n0 n1 n2 n3 t g lambda=0.88e-6; %单位m k0=2*pi/lambda; n0=3.4; n1=3.6; n2=3.4; n3=3.55 t=1.0e-6;%单位m g=0.2e-6;%单位m x0=[8.45i 7.5 8.45i 7 25.7]*1e6; options=optimset('MaxIter',1000,'MaxFunEvals',5000); x=fsolve(@nonl_K0i_beta0,x0,options) end function f=nonl_K0i_beta0(x) global k0 n0 n1 n2 n3 t g K00=x(1); K01=x(2); K02=x(3); K03=x(4); beta0=x(5); f(1)=K00^2-k0^2*n0^2+beta0^2; f(2)=K01^2-k0^2*n1^2+beta0^2; f(3)=K02^2-k0^2*n2^2+beta0^2; f(4)=K03^2-k0^2*n3^2+beta0^2; f(5)=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))); end [ Last edited by wh1125 on 2013-10-7 at 19:36 ] |
回复此楼» 猜你喜欢
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有124人回复
-
物理学I论文润色/翻译怎么收费?
已经有56人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
求用牛顿迭代法求非线性方程组
已经有4人回复
- 求助：matlab求解一个方程组，其中一个方程有三角函数
已经有10人回复
- Matlab数值积分结果是虚数
已经有6人回复
- matlab无约束非线性方程求最小值的问题
已经有7人回复
- 求助matlab 非线性方程组全部解！
已经有6人回复
- 哪位大神来指导下数字信号处理怎么学，感激不尽
已经有30人回复
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- 求助matlab求解三角函数方程
已经有9人回复
- 请教下：如何求一个函数在一定范围内的全部解
已经有10人回复
- matlab 求解隐式常微分方程时用ode15i 函数出现的问题求助
已经有3人回复
- 【致H兄】数列题的万能多项式拟合求法
已经有28人回复
- matlab求助2-关于欧拉公式
已经有11人回复
- 【求助】利用matlab求解微分方程得出解函数中含0的意思
已经有3人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
200