应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助：matlab求解一个方程组，最后一个方程有三角函数和虚数单位i
11/1返回列表
查看: 509  |  回复: 0
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 猜你喜欢
1楼 2013-10-07 16:59:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

科研从小木虫开始，人人为我，我为人人
相关版块跳转 我要订阅楼主 wh1125 的主题更新
11/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭