版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

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

514741735

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 102.7
散金: 33
帖子: 225
在线: 60.4小时
虫号: 2307572
注册: 2013-02-28
性别: GG
专业: 制造系统与自动化

[求助] 请问我求解的数值解合理吗？谁能帮我验证一下？已有1人参与

function q=myfun(p)
x=p(1);
y=p(2);
q(1)=abs((5*sin(x)+(200+5*x)*cos(x))*(5*cos(x)-(200+5*x)*sin(x))^(-1)*((200+5*x)*cos(x)-(200+5*y)*cos(y))-(200+5*x)*sin(x)+(200+5*y)*sin(y))*(((5*sin(x)+(200+5*x)*cos(x))*(5*cos(x)-(200+5*x)*sin(x))^(-1))^2+1)^(-1/2)-30;
q(2)=abs(((200+5*x)*sin(x)-5*cos(x))*(5*sin(x)+(200+5*x)*cos(x))^(-1)*((200+5*y)*cos(y)-(200+5*x)*cos(x))-(200+5*y)*sin(y)+(200+5*x)*sin(x))*((((200+5*x)*sin(x)-5*cos(x))*(5*sin(x)+(200+5*x)*cos(x))^(-1))^2+1)^(-1/2)-140;

x=fsolve('myfun',[3,4],optimset('Display','off'))

x =

28.1992 28.6208

bw263h2307572_1420714392_234.png

回复此楼

» 猜你喜欢

心脉受损已经有5人回复
博士读完未来一定会好吗已经有15人回复
Springer期刊投稿求助已经有4人回复
读博已经有3人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复

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

ode数值求解中重要的问题已经有9人回复
数值积分的求解已经有7人回复
求大神帮忙！我血小板只有正常值的10%，严重吗，我该怎么办啊？谁能帮帮我，给跪了！已经有21人回复
利用matlab求一阶变系数微分方程组的数值解。已经有7人回复
Matlab数值求解二阶常微分方程已经有9人回复
如何验证数值解的正确性已经有7人回复
一个椭圆型偏微分方程的数值解已经有19人回复
如何用数值解法求二阶微分方程的本征值已经有18人回复
数值求解N-S方程一般用什么编程啊？？？已经有9人回复
求高手解答Matlab求解微分方程的初始值和边界值问题已经有5人回复
抛物型偏微分方程的数值解求解问题，谢谢！已经有14人回复
求和符号相乘的方程怎么编程求数值解已经有7人回复
已知微分方程求解数值解已经有7人回复
“数值分析”、“数值解法”等中常出现的“数值”一词是什么意思？已经有17人回复
【求助】matlab 求解微分方程中的未知参数已经有20人回复
【求助】怎么用MATLAB求数值解啊（高手请帮忙啊）已经有6人回复
【求助】请教Matlab高手求定积分近似值问题已经有4人回复

自助者天助

1楼2015-01-08 21:11:13

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

不过需要说明的是，由于正弦余弦的周期性，该方程组的解不唯一。

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

4楼2015-01-08 21:41:28

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

查看全部 8 个回答

Carpathia

银虫 (小有名气)

应助: 35 (小学生)
金币: 5.6
散金: 12
红花: 2
帖子: 140
在线: 24.9小时
虫号: 1070633
注册: 2010-08-06
性别: GG
专业: 动力学与控制

把解代进去看看不行吗

赞一下

回复此楼

追求理想的道路上东张西望

2楼2015-01-08 21:19:31

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
514741735: 金币+30, ★★★★★最佳答案 2015-01-09 08:04:57

CODE:

function solve_nolin_eqs

clear all;clc

[x fval]=fsolve(@myfun,[3,4])

function q=myfun(p)

x=p(1);

y=p(2);

q(1)=abs((5*sin(x)+(200+5*x)*cos(x))*(5*cos(x)-(200+5*x)*sin(x))^(-1)*((200+5*x)*cos(x)-(200+5*y)*cos(y))-(200+5*x)*sin(x)+(200+5*y)*sin(y))*(((5*sin(x)+(200+5*x)*cos(x))*(5*cos(x)-(200+5*x)*sin(x))^(-1))^2+1)^(-1/2)-30;

q(2)=abs(((200+5*x)*sin(x)-5*cos(x))*(5*sin(x)+(200+5*x)*cos(x))^(-1)*((200+5*y)*cos(y)-(200+5*x)*cos(x))-(200+5*y)*sin(y)+(200+5*x)*sin(x))*((((200+5*x)*sin(x)-5*cos(x))*(5*sin(x)+(200+5*x)*cos(x))^(-1))^2+1)^(-1/2)-140;

x =

28.1992 28.6208

两个方程残差 fval =

1.0e-010 *

0.2747 -0.1253
可见方程残差是 1e-10 数量级，如果不考虑参数的具体物理意义，已经算是很好的解了。

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

3楼2015-01-08 21:34:42

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

514741735

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 102.7
散金: 33
帖子: 225
在线: 60.4小时
虫号: 2307572
注册: 2013-02-28
性别: GG
专业: 制造系统与自动化

引用回帖:

4楼: Originally posted by 月只蓝 at 2015-01-08 21:41:28
不过需要说明的是，由于正弦余弦的周期性，该方程组的解不唯一。

谢谢你，这样我就放心了

赞一下

回复此楼

自助者天助

5楼2015-01-09 08:03:54

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

查看全部 8 个回答