24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Littlecaesar

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 7.7
散金: 110
帖子: 65
在线: 55.4小时
虫号: 3225011
注册: 2014-05-22
性别: GG
专业: 生物材料

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

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

Untitled.png

回复此楼

» 猜你喜欢

反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有183人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

1楼 2014-06-12 22:25:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Littlecaesar

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 7.7
散金: 110
帖子: 65
在线: 55.4小时
虫号: 3225011
注册: 2014-05-22
性别: GG
专业: 生物材料

引用回帖:

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)