此微分方程組可用ode45（）函數(shù)求得其數(shù)值解 。
首先，建立自定義函數(shù)
function dy = rigid(t,y)
dy = zeros(3,1);
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
end
其二，用ode45（）函數(shù)求出x、y值
[x,y] = ode45(@rigid,[0 12],[0 1 1]);
其三，用plot（）函數(shù)繪出x—y，x—dy/dx，x—d2y/dx2曲線圖
plot(x,y(:,1),'-',x,y(:,2),'-.',x,y(:,3),'.')
運(yùn)行結(jié)果

猜你喜歡

• 清朝小說(shuō)家吳敬梓為什么寫(xiě)出了儒林外史 儒林外史作者是誰(shuí)
• 如何寫(xiě)出熱門(mén)短視頻標(biāo)題？熱門(mén)短視頻標(biāo)題模板
• 二階矩陣求逆，用伴隨矩陣方法。有點(diǎn)懵，求解答
• 廣汽匯理汽車(chē)金融怎么樣靠譜嗎？急求解答！！
• 求解析 高二地理會(huì)考
• 天珠變中都哪幾章有推倒？？求大神求解。
• 中政行測(cè)怎么樣啊 ，，求解。。。
• 沙姜與生姜有什么不同？求解
• 安裝CAD2012出現(xiàn)C++ 2005錯(cuò)誤 求解
• 句號(hào)打不出來(lái)了但是書(shū)名號(hào)卻能打出來(lái) 這是怎么回事啊 求解啊同樣都是一個(gè)按鍵 按鍵也沒(méi)壞啊！ 求大神幫忙