514741735

銀蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 102.7
散金: 33
帖子: 225
在線: 60.4小時
蟲號: 2307572
注冊: 2013-02-28
性別: GG
專業(yè): 制造系統(tǒng)與自動化

[求助] 請問我求解的數(shù)值解合理嗎？誰能幫我驗證一下？已有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ù)此樓

» 猜你喜歡

085700資源與環(huán)境308求調(diào)劑已經(jīng)有7人回復(fù)
265求調(diào)劑已經(jīng)有5人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有3人回復(fù)
求調(diào)劑已經(jīng)有3人回復(fù)
一志愿吉林大學材料學碩321求調(diào)劑已經(jīng)有12人回復(fù)
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有3人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù) 已經(jīng)有4人回復(fù)
材料學碩318求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中國海洋大學，生物學，301分，求調(diào)劑已經(jīng)有6人回復(fù)

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

ode數(shù)值求解中重要的問題已經(jīng)有9人回復(fù)
數(shù)值積分的求解已經(jīng)有7人回復(fù)
求大神幫忙！我血小板只有正常值的10%，嚴重嗎，我該怎么辦啊？誰能幫幫我，給跪了！已經(jīng)有21人回復(fù)
利用matlab求一階變系數(shù)微分方程組的數(shù)值解。已經(jīng)有7人回復(fù)
Matlab數(shù)值求解二階常微分方程已經(jīng)有9人回復(fù)
如何驗證數(shù)值解的正確性已經(jīng)有7人回復(fù)
一個橢圓型偏微分方程的數(shù)值解已經(jīng)有19人回復(fù)
如何用數(shù)值解法求二階微分方程的本征值已經(jīng)有18人回復(fù)
數(shù)值求解N-S方程一般用什么編程啊？？？已經(jīng)有9人回復(fù)
求高手解答Matlab求解微分方程的初始值和邊界值問題已經(jīng)有5人回復(fù)
拋物型偏微分方程的數(shù)值解求解問題，謝謝！已經(jīng)有14人回復(fù)
求和符號相乘的方程怎么編程求數(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ù)值解�。ǜ呤终垘兔Π。� 已經(jīng)有6人回復(fù)
【求助】請教Matlab高手求定積分近似值問題已經(jīng)有4人回復(fù)

自助者天助

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

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

Carpathia

銀蟲 (小有名氣)

應(yīng)助: 35 (小學生)
金幣: 5.6
散金: 12
紅花: 2
帖子: 140
在線: 24.9小時
蟲號: 1070633
注冊: 2010-08-06
性別: GG
專業(yè): 動力學與控制

把解代進去看看不行嗎

贊一下

回復(fù)此樓

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

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

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

月只藍

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

專家經(jīng)驗: +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(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

兩個方程殘差 fval =

1.0e-010 *

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

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

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

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

月只藍

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

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

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

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

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

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

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

514741735

銀蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 102.7
散金: 33
帖子: 225
在線: 60.4小時
蟲號: 2307572
注冊: 2013-02-28
性別: GG
專業(yè): 制造系統(tǒng)與自動化

引用回帖:

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

謝謝你，這樣我就放心了

贊一下

回復(fù)此樓

自助者天助

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

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

514741735

銀蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 102.7
散金: 33
帖子: 225
在線: 60.4小時
蟲號: 2307572
注冊: 2013-02-28
性別: GG
專業(yè): 制造系統(tǒng)與自動化

引用回帖:

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

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

贊一下

回復(fù)此樓

自助者天助

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

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

月只藍

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

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

引用回帖:

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

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

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

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

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

514741735

銀蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 102.7
散金: 33
帖子: 225
在線: 60.4小時
蟲號: 2307572
注冊: 2013-02-28
性別: GG
專業(yè): 制造系統(tǒng)與自動化

引用回帖:

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......;
...

謝謝指點，看來是我的思路不對了

贊一下

回復(fù)此樓

自助者天助

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

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

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

[考研] 265求調(diào)劑 +5	梁梁校校 2026-03-19	5/250	2026-03-20 05:54 by laoshidan
[考研] 288求調(diào)劑 +15	于海海海海 2026-03-19	15/750	2026-03-19 22:41 by 學員8dgXkO
[考研] 0856調(diào)劑，是學校就去 +6	sllhht 2026-03-19	7/350	2026-03-19 19:50 by 制度的
[考研] 0703化學調(diào)劑，六級已過，有科研經(jīng)歷 +12	曦熙兮 2026-03-15	12/600	2026-03-19 19:42 by maocaozhuxi
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-18	3/150	2026-03-19 14:42 by peike
[考研] 317求調(diào)劑 +3	申子申申 2026-03-19	6/300	2026-03-19 14:16 by 申子申申
[考研] 346求調(diào)劑[0856] +3	WayneLim327 2026-03-16	6/300	2026-03-19 11:21 by WayneLim327
[考研] 材料080500調(diào)劑求收留 +4	一顆meteor 2026-03-13	4/200	2026-03-19 10:32 by 30660438
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風 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
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 070300化學319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 302求調(diào)劑 +10	呼呼呼。。。。 2026-03-17	10/500	2026-03-18 12:45 by Linda Hu
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[基金申請] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky