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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 软件使用与操作 » 求高手解方程组（应该不难）
233/3返回列表上一页123
查看: 2646  |  回复: 22
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

Littlecaesar

银虫 (小有名气)
引用回帖:
20楼: Originally posted by waterque at 2014-06-15 10:57:24
应该可以吧，试试就知道了。...

大侠你好，可否继续请教一下，还是原来的方程，但这次问题不一样了。为了方便，我把问题重新描述一遍。求高手解方程组（应该不难）
赞一下 回复此楼
21楼2014-09-11 02:58:21
已阅   回复此楼   关注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楼2014-09-12 12:08:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Littlecaesar 的主题更新
233/3返回列表上一页123
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-17 04:26 by k06BKNblrGRH
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-17 04:04 by k06BKNblrGRH
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-17 03:40 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-17 03:16 by Tql5LQhh5rLK
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-17 03:14 by Tql5LQhh5rLK
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 8/400 2026-09-17 02:52 by Tql5LQhh5rLK
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-13 8/400 2026-09-17 02:04 by Tql5LQhh5rLK
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 02:02 by Tql5LQhh5rLK
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +7 3n8v2C8RimXI 2026-09-13 8/400 2026-09-17 01:38 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-16 17:27 by i3KBRD19ERkG
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:00 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:15 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 02:14 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭