24小时热门版块排行榜    

查看: 2199  |  回复: 22
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Littlecaesar

银虫 (小有名气)

[求助] 求高手解方程组(应该不难) 已有3人参与

下面的方程组,我用matlab里的fsolve求解的,但是每次结果都不准确,fval 值一个出来0, 另一个0.036,太大了。想把这个方程组解出来,求数学大神给点指点啊!

求高手解方程组(应该不难)
Untitled.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Littlecaesar

银虫 (小有名气)

引用回帖:
20楼: Originally posted by waterque at 2014-06-15 10:57:24
应该可以吧,试试就知道了。...

还有matlab代码:

K0=12e-5;
Eg=1;
tmax=0.52e-6;
delta_i=tmax/K0;
G=[0.1e-6 0.5e-6 1.0e-6 1.5e-6 2.0e-6 2.5e-6 3.0e-6];
for ii=1:length(G)
delta_f(ii)=2*G(ii)/tmax;
h=0.0004;
beta=sqrt(K0/(Eg*h));
L=125;
A(ii)=delta_i-delta_f(ii);
alpha(ii)=sqrt(tmax/(delta_f(ii)-delta_i)/(Eg*h));
fun=@(L1_c) A(ii)*cos(alpha(ii)*(L1_c-L))+(delta_i*beta*coth(beta*L1_c)/alpha(ii))*sin(alpha(ii)*(L-L1_c));
    L1_c0=125;
    opts = optimset('Diagnostics','off','Display','off');
    [L1_c,fval]=fzero(fun,L1_c0,opts)
    epsilon_m_c2 (ii)= -(A(ii)*alpha(ii)*sin(alpha(ii)*(L-L1_c))-(delta_i*beta*coth(beta*L1_c)/alpha(ii))*alpha(ii)*cos(alpha(ii)*(L-L1_c)));
    epsilon_m_c21(ii)=-A(ii)*alpha(ii)/sin(alpha(ii)*(L-L1_c));
end
figure(1)
plot(G*10^6,epsilon_m_c2*100,'ko-')
axis([0 3.5 0 14])
set(gca,'Ytick',[2 4 6 8 10 12])
xlabel('G_c(N/m)')
ylabel('\epsilon_c_2(%)')
set(gca,'FontSize',12)
set(get(gca,'XLabel'),'FontSize',12)
set(get(gca,'YLabel'),'FontSize',12)
22楼2014-09-11 03:00:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
B值没给?还有出了图片最好再给出文本。
2楼2014-06-12 23:10:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Littlecaesar

银虫 (小有名气)

谢谢楼上,刚才是没有给出B值,重新写了一下。BW,什么文本?不懂。。。
求高手解方程组(应该不难)-1
Untitled.png

3楼2014-06-12 23:27:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Littlecaesar

银虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2014-06-12 23:10:28
B值没给?还有出了图片最好再给出文本。

改过来了,,,BW,什么文本?
4楼2014-06-12 23:29:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见