24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 532  |  回复: 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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 wh1125 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-05 9/450 2026-04-05 22:24 by dongzh2009
[考研] 316求调剂 +5 yyx想调剂 2026-04-05 5/250 2026-04-05 22:22 by 咔咔咔咔9
[考研] 270分求调剂 +4 maxjxbsk 2026-04-01 4/200 2026-04-05 17:04 by yulian1987
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[考研] 298求调剂 +7 manman511 2026-04-05 7/350 2026-04-05 10:29 by 唐沐儿
[考研] 一志愿江南大学085501机械工程专硕326分,本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 11408,335分,本科211,求调剂,可转专业 +5 鳄梨大鳄鱼 2026-04-03 5/250 2026-04-04 22:49 by chongya
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] 材料调剂 +15 一样YWY 2026-04-01 15/750 2026-04-04 22:23 by hemengdong
[考研] 环境285分,过六级,求调剂 +10 xhr12 2026-04-02 10/500 2026-04-04 21:53 by bn53987
[考研] 化学调剂求助 +6 LULONG1 2026-04-03 6/300 2026-04-03 23:13 by qzxyhcsy
[考研] 334求调剂 +9 Trying] 2026-03-31 9/450 2026-04-03 15:18 by 琢珥丶
[考研] 085501一志愿天工大,机械专硕求调剂,跨材料 +3 33上 2026-04-03 3/150 2026-04-03 14:08 by 1753564080
[考研] 321求调剂 +17 y-yh 2026-04-01 20/1000 2026-04-03 12:57 by y-yh
[考研] 366求调剂 +7 sbdnd 2026-04-03 7/350 2026-04-03 12:40 by cymywx
[考研] 085600 295分求调剂 +19 W55j 2026-03-30 23/1150 2026-04-03 09:53 by 千千运气
[考博] 26年申博 +3 staryer 2026-03-30 4/200 2026-04-01 23:21 by ai4pharm
[考研] 285求调剂 +7 AZMK 2026-03-30 13/650 2026-04-01 17:00 by 七度不信任
[考研] 一志愿北交材料工程总分358 +5 cs0106 2026-04-01 7/350 2026-04-01 11:45 by wangjy2002
信息提示
请填处理意见