【答案】应助回帖

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有62人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助高手解方程（100金币不够可以再加）已经有8人回复
matlab如何求解多元多次方程，包括绘图已经有7人回复
关于病态方程组求解的疑问已经有5人回复
用matlab 求解参数方程组已经有8人回复
求高手帮我解一阶常微分方程组，非常感谢！！已经有5人回复
求助解复杂非线性方程组的好的方法已经有24人回复
用matlab求解一个非线性方程组的解已经有4人回复
求助-帮我用mathematics计算下下面的方程组已经有13人回复
求助，用matlab解微分方程组，希望高手能给予指点已经有17人回复
这个方程组怎么解啊已经有14人回复
matlab求解高阶方程组已经有3人回复
微分方程组求参数问题，求高人指教，感谢已经有11人回复
求高手解答Matlab求解微分方程的初始值和边界值问题已经有5人回复
求高手解一个二元二次方程组的解析解已经有9人回复
求Matlab高手解决线性方程组的迭代求解问题已经有17人回复
求助 MATLAB解方程组-fslove 已经有7人回复
求助（已经-无量纲化）的非线性方程组（已经仔细检查，能算，但exitflag=0求高手已经有4人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
请教Matlab高手解多元非线性方程组中的参数已经有6人回复
【求助】同余方程组的解已经有11人回复
【求助】求教matlab解非线性方程组已经有9人回复

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

回帖支持 ( 显示支持度最高的前 50 名 )

waterque

金虫 (正式写手)

应助: 59 (初中生)
金币: 1989.1
散金: 96
红花: 10
帖子: 759
在线: 191.5小时
虫号: 2095273
注册: 2012-10-30
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★
感谢参与，应助指数 +1
fegg7502: 金币+1, 鼓励交流 2014-09-11 07:00:12

设x = L2 - L1,
f(x) = Acos(ax) + Bsin(ax) % a ---> alpha
则f'(x) = -Aasin(ax) + Bacos(ax)
那么你的问题就变成：
f(x) = 0;
f'(x) = epsilon_m;

并不能得到具体的L1, L2.

赞一下(1人)

行自己的路：waterque@hotmail.com

12楼2014-06-13 18:23:44

waterque

金虫 (正式写手)

应助: 59 (初中生)
金币: 1989.1
散金: 96
红花: 10
帖子: 759
在线: 191.5小时
虫号: 2095273
注册: 2012-10-30
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★
fegg7502: 金币+1, 鼓励交流 2014-09-11 07:00:20

引用回帖:

12楼: Originally posted by waterque at 2014-06-13 18:23:44
设x = L2 - L1,
f(x) = Acos(ax) + Bsin(ax) % a ---> alpha
则f'(x) = -Aasin(ax) + Bacos(ax)
那么你的问题就变成：
f(x) = 0;
f'(x) = epsilon_m;

并不能得到具体的L1, L2.

matlab代码：
clear;
A = -0.30714;
B = 0.12966;
a = 0.0206;
epm = linspace(0.07,0.1,10);
x = -1500:0.01:1500;
f = A*cos(a*x)+B*sin(a*x);
h = -B*a*cos(a*x)+A*a*sin(a*x);
plot(x,f)
hold on;
plot(x,x*0,'r');

for k = 1:length(epm)
plot(x,h+epm(k),'g');
end
hold off;

说明：分别画出两个函数的曲线来，其中一个是f(x), 另一个是 h(x) = -f'(x) + epsilon_m;
方程组的解是这两个曲线相交于横轴的点。但从图中看的话，h(x)最小的值也是在横轴上方的，所以无实数解。

curves.jpg

赞一下(1人)

行自己的路：waterque@hotmail.com

13楼2014-06-13 18:54:15

普通回帖

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

感谢参与，应助指数 +1

B值没给？还有出了图片最好再给出文本。

2楼2014-06-12 23:10:28

Littlecaesar

银虫 (小有名气)

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

谢谢楼上，刚才是没有给出B值，重新写了一下。BW，什么文本？不懂。。。

Untitled.png

3楼2014-06-12 23:27:27

Littlecaesar

银虫 (小有名气)

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

引用回帖:

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

改过来了，，，BW，什么文本？

4楼2014-06-12 23:29:44

NicoleLas

木虫 (正式写手)

应助: 66 (初中生)
金币: 4288.7
红花: 26
帖子: 514
在线: 177.2小时
虫号: 2872321
注册: 2013-12-15
性别: GG
专业: 计算机软件

引用回帖:

4楼: Originally posted by Littlecaesar at 2014-06-12 23:29:44
改过来了，，，BW，什么文本？...

文本格式代码...............
你贴个图片难道还指望可能帮助你的人对照着图片输入公式一遍？
为自己，为他人，都得考虑一下，dingd版主算是很热心的了，如果你直接给出了文本代码，他可以很方便地用1stOpt最新版给你计算一下，很可能很快就有结果了。

Be water, my friend.

5楼2014-06-12 23:57:18

Littlecaesar

银虫 (小有名气)

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

引用回帖:

5楼: Originally posted by NicoleLas at 2014-06-12 23:57:18
文本格式代码...............
你贴个图片难道还指望可能帮助你的人对照着图片输入公式一遍？
为自己，为他人，都得考虑一下，dingd版主算是很热心的了，如果你直接给出了文本代码，他可以很方便地用1stOpt ...

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

6楼2014-06-13 00:08:33