| 5 | 1/1 | 返回列表 |
| 查看: 1125 | 回復: 3 | ||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
jyypzfswty新蟲 (初入文壇)
|
[求助]
求怎么用matlab解這個方程 已有2人參與
|
|||
|
要解個方程,但我對matlab不熟,請各位大神幫寫個程序代碼? 0.5sin(3.6x)cos(x)+3.6cos(3.6x)sin(x)=0 (x>0) 有無窮多解,求100個正的解。 |
主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗: +1059 |

木蟲 (正式寫手)
木蟲 (正式寫手)
|
%%這是非線性方程求解,普通的solve無法求得準確的解,fsolve無法解得100個解; %%我這里用的是數(shù)值解法,精確度在為10^(-4),可改變x來提高精確度,看個人需要; %%思路:求出y=0.5*sin(3.6*x).*cos(x)+3.6*cos(3.6*x).*sin(x)=0的100個解X; %可畫出其圖像,很容易看出y會在解處變號!。∮纱丝珊Y選出解 x=0:0.0001:70;%間隔越小,精確度越高 y=0.5*sin(3.6*x).*cos(x)+3.6*cos(3.6*x).*sin(x);%函數(shù) plot(x,y,'b') hold on n=1; for i=1 length(x)-1)if y(i)*y(i+1)<=0%如果小于0,知道此處y發(fā)生變號,那么此處對應的x即解 X(n)=x(i);%將找到的解賦值給另一個矩陣 n=n+1; end end X=X(2:101);%因為問題說前100個正解,不包括0 Y=0.5*sin(3.6*X).*cos(X)+3.6*cos(3.6*X).*sin(X); plot(X,Y,'.r')%畫出解的點 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
[考研] 322求調(diào)劑 +4 | 舊吢 2026-03-24 | 4/200 |
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +4 | 0857求調(diào)劑. 2026-03-24 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | question挽風 2026-03-23 | 3/150 |
|
|
[考研] 材料求調(diào)劑 +8 | @taotao 2026-03-21 | 8/400 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研] 一志愿北化求調(diào)劑 +3 | Jsman 2026-03-22 | 3/150 |
|
|
[考研] 279求調(diào)劑 +6 | 紅衣隱官 2026-03-21 | 6/300 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當牛 2026-03-24 | 4/200 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|