怎麼用MATLAB程式設計求解單自由度質量

2025-07-14 02:40:16 字數 1708 閱讀 4336

1樓:匿名使用者

在網上搜了一下,不知道這個跟你的問題有沒有聯絡,好歹不能讓你零答案啊。**分割法:#include#includefloatd,e,f;voidmain()else}if(f1

用matlab求解下面方程,請編寫程式,用m.檔案

2樓:書宬

function y=zhongy(x) %儲存為 檔案。

y=[2*x(1)-x(2)-exp(-x(1));

x(1)+2*x(2)-exp(-x(2))]

x,fval,exitflag]=fsolve('zhongy',[;

x 為所求的解。

fval為解代入方程的驗證結果。

exitflag=1表示函式收斂於解;反之,=0

3樓:網友

輸入[x1,x2]=solve('2*x1-x2=exp(-x1)','-x1+2*x2=exp(-x2)')就可以求解了。

解出來是。x1 =

x2 =即x1=x2=exp(-x1)

用matlab程式設計 求零輸入響應

4樓:網友

零輸入的含義是 f(t)=0,直接用dsolve求解即可(其中題目第三個方程第一項有誤,應為三階導數y'''(t)):

dsolve('d2y+5*dy+6*y=0','y(0)=1,dy(0)=-1')

ans =2*exp(-2*t)-exp(-3*t)

dsolve('d2y+2*dy+5*y=0','y(0)=2,dy(0)=-2')

ans =2*exp(-t)*cos(2*t)

dsolve('d3y+4*d2y+5*dy+2*y=0','y(0)=0,dy(0)=-1,d2y(0)=-1')

ans =3*exp(-2*t)+3*exp(-t)-4*exp(-t)*t

多自由度振動自振頻率用matlab怎麼求

5樓:_愛成繞指柔

假如:m=diag([1,2,3])%質量矩陣k=[1,1,2;1,3,0;2,0,5]%剛度矩陣%以上是根據你的實際問題寫出的質量剛度矩陣a=inv(m)*k;

v,d]=eig(a); 計算特版徵向量和特徵值權%計算固有頻率並排序。

la=diag(d); 提取特徵值ww=sqrt(la); 提取固有頻率w=sort(ww) %固有頻率排序,即為所求頻率%進一步的可以求振型。

提取特徵向量並排序。

n=length(m) ;矩陣維度for j=1:n

for i=1:n

if w(j)==ww(i);

x(:,j)=v(:,i)/max(v(:,i));

endend

end %x求出應為n行n列矩陣,其中每一列都是乙個振型,一共三個,而且已經排序。

x %輸出振型矩陣結果如下:

請點選輸入**描述。

6樓:碌滯藝吵碌

從你的程式中,baiv的列。

du向量是振動的固有模態zhi,每個元素屬於不同的自由dao度,因為你取了專9個,屬所以新增橫座標就可以了,我覺得直接用x=1:9就行,能表示出趨勢就可以 for j=1:9; plot(x,v(:

j));hold on; end;

怎麼用MATLAB程式設計求方程組的解

b 4 3 9 8 a1 a b r rref a1 r 1.0000 0 0 0 1.4841 0 1.0000 0 0 0.6816 0 0 1.0000 0 0.5337 0 0 0 1.0000 1.2429 最後一列為方程的抄解,跟樓上的那位仁兄的好像有點差別,不過我把我的解代入原方程檢驗...

LU法求解線性方程組,matlab程式設計

a 1 4 0 1 1 5 1 0 3 2 0 3 4 0 1 4 b 11 12 7 5 l,u lu a 計算上三角和下三角x u l b 求解x matlab中怎樣利用lu分解求解方程組 如果a,b是引數,那麼x,y是未知數一個方程兩個未知數,有無限多的解啊換句話說是二維座標系中符合x 2 y...

matlab怎麼用有限差分法求解方程

fungecho 刪除舉報此資訊 如圖中,需要解出力的大小。求 matlab 程式,用 有限差分法 解 橢圓偏微分方程.題目如圖.建議看看這個 利用有限差分和matlab矩陣運算直接求解二維泊松.和你要的原理是完全一樣的.matlab程式.自己寫吧.matlab求解偏微分方程組,查了一下,說要用有限...