514741735

應(yīng)助: 2 (幼兒園)
金幣: 102.7
散金: 33
帖子: 225
在線(xiàn): 60.4小時(shí)
蟲(chóng)號(hào): 2307572
注冊(cè): 2013-02-28
性別: GG
專(zhuān)業(yè): 制造系統(tǒng)與自動(dòng)化

[求助] 請(qǐng)問(wèn)我求解的數(shù)值解合理嗎？誰(shuí)能幫我驗(yàn)證一下？已有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

回復(fù)此樓

» 猜你喜歡

一志愿西南交通專(zhuān)碩材料355 本科雙非求調(diào)劑已經(jīng)有3人回復(fù)
295復(fù)試調(diào)劑已經(jīng)有4人回復(fù)
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有8人回復(fù)
材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦已經(jīng)有4人回復(fù)
能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有14人回復(fù)
0856調(diào)劑，是學(xué)校就去已經(jīng)有8人回復(fù)
材料080500調(diào)劑求收留已經(jīng)有5人回復(fù)
271材料工程求調(diào)劑已經(jīng)有7人回復(fù)
265求調(diào)劑已經(jīng)有8人回復(fù)
085410人工智能專(zhuān)碩317求調(diào)劑（0854都可以）已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

ode數(shù)值求解中重要的問(wèn)題已經(jīng)有9人回復(fù)
數(shù)值積分的求解已經(jīng)有7人回復(fù)
求大神幫忙！我血小板只有正常值的10%，嚴(yán)重嗎，我該怎么辦��？誰(shuí)能幫幫我，給跪了！已經(jīng)有21人回復(fù)
利用matlab求一階變系數(shù)微分方程組的數(shù)值解。已經(jīng)有7人回復(fù)
Matlab數(shù)值求解二階常微分方程已經(jīng)有9人回復(fù)
如何驗(yàn)證數(shù)值解的正確性已經(jīng)有7人回復(fù)
一個(gè)橢圓型偏微分方程的數(shù)值解已經(jīng)有19人回復(fù)
如何用數(shù)值解法求二階微分方程的本征值已經(jīng)有18人回復(fù)
數(shù)值求解N-S方程一般用什么編程��？？？已經(jīng)有9人回復(fù)
求高手解答Matlab求解微分方程的初始值和邊界值問(wèn)題已經(jīng)有5人回復(fù)
拋物型偏微分方程的數(shù)值解求解問(wèn)題，謝謝！已經(jīng)有14人回復(fù)
求和符號(hào)相乘的方程怎么編程求數(shù)值解已經(jīng)有7人回復(fù)
已知微分方程求解數(shù)值解已經(jīng)有7人回復(fù)
“數(shù)值分析”、“數(shù)值解法”等中常出現(xiàn)的“數(shù)值”一詞是什么意思？已經(jīng)有17人回復(fù)
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復(fù)
【求助】怎么用MATLAB求數(shù)值解�。ǜ呤终�(qǐng)幫忙�。� 已經(jīng)有6人回復(fù)
【求助】請(qǐng)教Matlab高手求定積分近似值問(wèn)題已經(jīng)有4人回復(fù)

自助者天助

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Carpathia

銀蟲(chóng) (小有名氣)

應(yīng)助: 35 (小學(xué)生)
金幣: 5.6
散金: 12
紅花: 2
帖子: 140
在線(xiàn): 24.9小時(shí)
蟲(chóng)號(hào): 1070633
注冊(cè): 2010-08-06
性別: GG
專(zhuān)業(yè): 動(dòng)力學(xué)與控制

把解代進(jìn)去看看不行嗎

贊一下

回復(fù)此樓

追求理想的道路上東張西望

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專(zhuān)家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線(xiàn): 3291.4小時(shí)
蟲(chóng)號(hào): 1122189
注冊(cè): 2010-10-14
專(zhuān)業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +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

兩個(gè)方程殘差 fval =

1.0e-010 *

0.2747 -0.1253
可見(jiàn)方程殘差是 1e-10 數(shù)量級(jí)，如果不考慮參數(shù)的具體物理意義，已經(jīng)算是很好的解了。

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專(zhuān)家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線(xiàn): 3291.4小時(shí)
蟲(chóng)號(hào): 1122189
注冊(cè): 2010-10-14
專(zhuān)業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

【答案】應(yīng)助回帖

不過(guò)需要說(shuō)明的是，由于正弦余弦的周期性，該方程組的解不唯一。

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

514741735

銀蟲(chóng) (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 102.7
散金: 33
帖子: 225
在線(xiàn): 60.4小時(shí)
蟲(chóng)號(hào): 2307572
注冊(cè): 2013-02-28
性別: GG
專(zhuān)業(yè): 制造系統(tǒng)與自動(dòng)化

引用回帖:

4樓: Originally posted by 月只藍(lán) at 2015-01-08 21:41:28
不過(guò)需要說(shuō)明的是，由于正弦余弦的周期性，該方程組的解不唯一。

謝謝你，這樣我就放心了

贊一下

回復(fù)此樓

自助者天助

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

514741735

銀蟲(chóng) (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 102.7
散金: 33
帖子: 225
在線(xiàn): 60.4小時(shí)
蟲(chóng)號(hào): 2307572
注冊(cè): 2013-02-28
性別: GG
專(zhuān)業(yè): 制造系統(tǒng)與自動(dòng)化

引用回帖:

4樓: Originally posted by 月只藍(lán) at 2015-01-08 21:41:28
不過(guò)需要說(shuō)明的是，由于正弦余弦的周期性，該方程組的解不唯一。

請(qǐng)問(wèn)，我在編程的時(shí)候是參照例子做的二元方程，變量a和b用的是確定值替換的，如果我想把a(bǔ)=5,b=200作為方程，列四元方程組，程序應(yīng)該怎么編呢？謝謝了

贊一下

回復(fù)此樓

自助者天助

6樓2015-01-09 08:16:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專(zhuān)家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線(xiàn): 3291.4小時(shí)
蟲(chóng)號(hào): 1122189
注冊(cè): 2010-10-14
專(zhuān)業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

引用回帖:

6樓: Originally posted by 514741735 at 2015-01-09 08:16:13
請(qǐng)問(wèn)，我在編程的時(shí)候是參照例子做的二元方程，變量a和b用的是確定值替換的，如果我想把a(bǔ)=5,b=200作為方程，列四元方程組，程序應(yīng)該怎么編呢？謝謝了...

不需要列成四元方程組吧。按如下把a(bǔ) 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*......;

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

7樓2015-01-09 08:41:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

514741735

銀蟲(chóng) (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 102.7
散金: 33
帖子: 225
在線(xiàn): 60.4小時(shí)
蟲(chóng)號(hào): 2307572
注冊(cè): 2013-02-28
性別: GG
專(zhuān)業(yè): 制造系統(tǒng)與自動(dòng)化

引用回帖:

7樓: Originally posted by 月只藍(lán) at 2015-01-09 08:41:43
不需要列成四元方程組吧。按如下把a(bǔ) 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......;
...

謝謝指點(diǎn)，看來(lái)是我的思路不對(duì)了

贊一下

回復(fù)此樓

自助者天助

8樓2015-01-09 08:56:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +3	@taotao 2026-03-20	3/150	2026-03-20 09:01 by 每天只擺一小會(huì)
[考研] 317求調(diào)劑 +3	申子申申 2026-03-19	7/350	2026-03-20 08:53 by 申子申申
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 【考研調(diào)劑】化學(xué)專(zhuān)業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑 +5	吃吃吃才有意義 2026-03-19	5/250	2026-03-19 16:18 by 30660438
[考研] 一志愿天津大學(xué)化學(xué)工藝專(zhuān)業(yè)（081702）315分求調(diào)劑 +11	yangfz 2026-03-17	11/550	2026-03-19 15:06 by houyaoxu
[考研] 317求調(diào)劑 +3	申子申申 2026-03-19	6/300	2026-03-19 14:16 by 申子申申
[考研] 一志愿天大材料與化工（085600）總分338 +5	蔡大美女 2026-03-13	5/250	2026-03-19 10:44 by 是小劉呀～
[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 311求調(diào)劑 +6	26研0 2026-03-15	6/300	2026-03-18 14:43 by haxia
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 0854，計(jì)算機(jī)類(lèi)招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 材料工程專(zhuān)碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 0854控制工程 359求調(diào)劑可跨專(zhuān)業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82