| 查看: 621 | 回复: 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 ] |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有99人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
找到一些相关的精华帖子，希望有用哦~
求用牛顿迭代法求非线性方程组
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
20