%%-----------------------------

syms t w dt real %定義符號(hào)變量，real表示符號(hào)變量定義在實(shí)數(shù)范圍內(nèi)。

dt=pi/180,w=1*5*pi/9, %%--角度要變換成弧度，w有三個(gè)值，改一下，計(jì)算一次。

for k=1:40 %由于本方程不容易得到解析解，所以，用循環(huán)變量，k從1增加到50，每步默認(rèn)增加1。

f=k, %把變化的循環(huán)變量賦值給f舵機(jī)頻響。

c=1/(4*t*f), %------------------------------------------

s=w*t*((c-1)+exp((-1*c)))-dt/2^0。5, %%%把關(guān)系式列出。

ss=subs(s), % 把s中的已知的變量(包括具體數(shù)值和符號(hào)表達(dá)式)，代入。

ff=solve(ss,'t') , %解方程，求ss=0的解，把舵機(jī)時(shí)間常數(shù)t看作變量解方程。

t1=double(ff(1)),% w=100時(shí)，循環(huán)到36，即f=36，方程無(wú)解，循環(huán)被終止，這時(shí)，只得到35個(gè)值，直接plot()畫(huà)圖即可。

%t2=double(ff(2)), %由于上式解方程中，有兩值，所以，分開(kāi)賦值給t1、t2并把它們變換為雙精度數(shù)值，以防止以分?jǐn)?shù)形式給出。

fx(k)=f,tx1(k)=t1,%tx2(k)=t2, %把變量賦值給能循環(huán)記錄的數(shù)組。

end %循環(huán)計(jì)算完畢，結(jié)束。

fx=double(fx),

hold on, %固定，以后的圖象都在一張圖上。

plot(fx,tx1), %畫(huà)出ff一個(gè)解的圖象。

▲▲▲▲▲▲▲▲▲▲▲▲

程序已經(jīng)修改完畢，花點(diǎn)時(shí)間自己看懂了，以后就可以隨你修改了。(4。4。9)中的有錯(cuò)誤，上面的等式中分母少了一個(gè)4，依據(jù)圖象下面的等式是對(duì)的。

全部

總結(jié)

以上是生活随笔為你收集整理的matlab程序怎么改,修改matlab程序的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。