| 24小時(shí)熱門(mén)版塊排行榜 |
| 查看: 2327 | 回復(fù): 22 | ||
Littlecaesar銀蟲(chóng) (小有名氣)
|
[求助]
求高手解方程組(應(yīng)該不難) 已有3人參與
|
|
|
下面的方程組,我用matlab里的fsolve求解的,但是每次結(jié)果都不準(zhǔn)確,fval 值一個(gè)出來(lái)0, 另一個(gè)0.036,太大了。想把這個(gè)方程組解出來(lái),求數(shù)學(xué)大神給點(diǎn)指點(diǎn)啊! Untitled.png |
金蟲(chóng) (正式寫(xiě)手)

金蟲(chóng) (正式寫(xiě)手)

鐵桿木蟲(chóng) (職業(yè)作家)
銀蟲(chóng) (小有名氣)
銀蟲(chóng) (小有名氣)
木蟲(chóng) (正式寫(xiě)手)

銀蟲(chóng) (小有名氣)
|
clear clc K0=9.5e-5; Eg=1; Em=0.003; tmax=0.52e-6; delta_i=tmax/K0; G=0.8e-6; delta_f=2*G/tmax; h=0.0004; beta=sqrt(K0/(Eg*h)); L=100; epsilon_m=linspace(0,0.1,11) % 1st critical strain epsilon_m_c1=delta_i*beta*coth(beta*L) % 2nd critical strain A=delta_i-delta_f alfa=sqrt(tmax/(delta_f-delta_i)/(Eg*h)) fun=@(L1_c) A*cos(alfa*(L1_c-L))+(delta_i*beta*coth(beta*L1_c)/alfa)*sin(alfa*(L-L1_c)); L1_c0=100; opts = optimset('Diagnostics','off','Display','off'); [L1_c,fval]=fzero(fun,L1_c0,opts) epsilon_m_c2 =-(A*alfa*sin(alfa*(L-L1_c))-(delta_i*beta*coth(beta*L1_c)/alfa)*alfa*cos(alfa*(L-L1_c))) x=linspace(0,L,100); for ii=1:11 if epsilon_m(ii)>epsilon_m_c2 % calculate L1 and L2 fun=@(Lc) [ A*alfa*sin(alfa*(Lc(2)-Lc(1)))-(delta_i*beta*coth(beta*Lc(1))/alfa)*alfa*cos(alfa*(Lc(2)-Lc(1)))+epsilon_m(ii); ... A*cos(alfa*(Lc(1)-Lc(2)))+(delta_i*beta*coth(beta*Lc(1))/alfa)*sin(alfa*(Lc(2)-Lc(1)))]; L0=[100,100]; opts = optimset('TolFun',1e-7,'TolX',1e-20); [Lc,fval,exitflag]=fsolve(fun,L0,opts) L1=Lc(1); L2=Lc(2); for jj=1:100 % region I if x(jj)<L1|x(jj)==L1 epsilon_g(jj)=epsilon_m(ii)-delta_i*beta*cosh(beta*x(jj))/sinh(beta*L1); % region II elseif L1<x(jj)&x(jj)<L2 B=delta_i*beta*coth(beta*L1)/alfa; epsilon_g(jj)=epsilon_m(ii)+A*alfa*sin(alfa*(x(jj)-L1))-B*alfa*cos(alfa*(x(jj)-L1)); elseif x(jj)>L2 epsilon_g(jj)=0; end end plot(x,epsilon_g,'g') hold on axis([0,200,0,0.12]) end end |
銀蟲(chóng) (小有名氣)
木蟲(chóng) (正式寫(xiě)手)

銀蟲(chóng) (小有名氣)
鐵桿木蟲(chóng) (職業(yè)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|