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论文润色/翻译怎么收费? 已经有289人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

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

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

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)最小的值也是在横轴上方的，所以无实数解。