24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1181  |  回复: 7

514741735

银虫 (小有名气)

[求助] 请问我求解的数值解合理吗?谁能帮我验证一下? 已有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
回复此楼

» 猜你喜欢

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

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

Carpathia

银虫 (小有名气)

把解代进去看看不行吗
追求理想的道路上东张西望
2楼2015-01-08 21:19:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

不过需要说明的是,由于 正弦 余弦 的周期性,该方程组的解不唯一。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
4楼2015-01-08 21:41:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

514741735

银虫 (小有名气)

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

谢谢你,这样我就放心了
自助者天助
5楼2015-01-09 08:03:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

514741735

银虫 (小有名气)

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

请问,我在编程的时候是参照例子做的二元方程,变量a和b用的是确定值替换的,如果我想把a=5,b=200作为方程,列四元方程组,程序应该怎么编呢?谢谢了
自助者天助
6楼2015-01-09 08:16:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
6楼: Originally posted by 514741735 at 2015-01-09 08:16:13
请问,我在编程的时候是参照例子做的二元方程,变量a和b用的是确定值替换的,如果我想把a=5,b=200作为方程,列四元方程组,程序应该怎么编呢?谢谢了...

不需要列成四元方程组吧。按如下把a b 定义好即可。
function solve_nolin_eqs
clear all;clc
[x fval]=fsolve(@myfun,[3,4])

function q=myfun(p)
a=5;
b=200;
x=p(1);
y=p(2);
q(1)=abs((a*sin(x)+(b+a*x......;
q(2)=abs(((b+a*x)*sin(x)-a*......;
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
7楼2015-01-09 08:41:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

514741735

银虫 (小有名气)

引用回帖:
7楼: Originally posted by 月只蓝 at 2015-01-09 08:41:43
不需要列成四元方程组吧。按如下把a b 定义好即可。
function solve_nolin_eqs
clear all;clc
=fsolve(@myfun,)

function q=myfun(p)
a=5;
b=200;
x=p(1);
y=p(2);
q(1)=abs((a*sin(x)+(b+a*x......;
...

谢谢指点,看来是我的思路不对了
自助者天助
8楼2015-01-09 08:56:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 514741735 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 320求调剂 +6 深郊akm 2026-04-17 6/300 2026-04-20 18:57 by fs26jie
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +4 yexuqing 2026-04-19 4/200 2026-04-20 14:47 by brantleo
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 通信工程求调剂!!! +7 zlb770521 2026-04-14 7/350 2026-04-19 20:56 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂推荐 +9 小聂爱学习 2026-04-14 9/450 2026-04-19 17:03 by 中豫男
[考研] 291求调剂 +12 关忆北. 2026-04-14 13/650 2026-04-19 16:50 by 中豫男
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 297,工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 22408 312求调剂 +24 门路摸摸 2026-04-14 26/1300 2026-04-18 13:04 by wunaiy88
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
信息提示
请填处理意见