




版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、符號(hào)積分變換傅里葉變換及其反變換1. 傅里葉變換 f=f(x)F=F(w) syms x w u vf=sin(x)*exp(-x2);F1=fourier(f) F1 = transform:fourier(sin(x)/exp(x2), x, -w) >> f=x;F2=fourier(f) F2 = pi*dirac(w, 1)*2*i >> h=x*exp(-abs(x);F3=fourier(h) F3 = -(w*4*i)/(w2 + 1)2 >> h=x*exp(-abs(x);F3=fourier(h,u) F3 =-(u*4*i)/(u2
2、+ 1)22. 傅里葉反變換 syms w v x tg=exp(-abs(x);IF2=ifourier(g) IF2 = 1/(pi*(t2 + 1)拉普拉斯變換及其反變換1. 拉普拉斯變換 syms x s t v f1=sqrt(t);L1=laplace(f1) L1 =pi(1/2)/(2*s(3/2)2. 拉普拉斯反變換 syms a s t u v x f=exp(x/s2);IL1=ilaplace(f) IL1 =ilaplace(exp(x/s2), s, t)Z變換及其反變換方程的解析解線性方程組的解析解包括求解線性方程組和非線性方程組的函數(shù)solve(),也有求解常微
3、分方程組的函數(shù)dsolve() L1='x+y+z=10'L2='3*x+2*y+z=14'L3='2*x+3*y-z=1'; %L1、L2、L3分別是三個(gè)字符串g=solve(L1,L2,L3)g = x: 1x1 sym y: 1x1 sym z: 1x1 sym %表明g是一個(gè)結(jié)構(gòu)數(shù)組,其中每個(gè)元素為一>> g.x %符號(hào)類型的量,用如下方法查看方程解的具體值 ans =1一般求解方法: L1='x+y+z=10'1 / 6L2='3*x+2*y+z=14'L3='2*x+3*y-z=1&
4、#39;x y z=solve(L1,L2,L3) x =1 y =2 z =7線性方程組的解析解>> f=sym('a*x2+b*x+c=0');xf=solve(f) xf = -(b + (b2 - 4*a*c)(1/2)/(2*a) -(b - (b2 - 4*a*c)(1/2)/(2*a)>>syms y z u v w eq1=u*y2+v*z+w;eq2=y+z+w;y z=solve(eq1,eq2,y,z) y = (v + 2*u*w + (v2 + 4*u*w*v - 4*u*w)(1/2)/(2*u) - w (v + 2*u*w
5、 - (v2 + 4*u*w*v - 4*u*w)(1/2)/(2*u) - w z = -(v + 2*u*w + (v2 + 4*u*w*v - 4*u*w)(1/2)/(2*u) -(v + 2*u*w - (v2 + 4*u*w*v - 4*u*w)(1/2)/(2*u)常微分方程組的解析解在微分方程組的表達(dá)式equ中,大寫(xiě)字母D表示對(duì)自變量(設(shè)為x)的微分算子:D=d/dx,D2=d2/d2x.微分算子D后面的字母則表示為因變量,即帶求解的未知函數(shù)。>> y=dsolve('Dy+a*x=0','x') y =C2 - (a*x2)/2&g
6、t;> y=dsolve('D2y+2*x=2*y','x') y = x + C4*exp(2(1/2)*x) + C5/exp(2(1/2)*x)>> y=dsolve('D2y+2*x=2*y','y(2)=5','Dy(1)=2','x') y =x + (exp(2(1/2)*x)*(6*exp(2(1/2) + 2(1/2)/(2*exp(2(1/2)*(exp(2*2(1/2) + 1) - (2(1/2)*exp(2*2(1/2)*(exp(2(1/2) - 3*2
7、(1/2)/(2*exp(2(1/2)*x)*(exp(2*2(1/2) + 1)MATLAB程序設(shè)計(jì)全局變量 global A B C 變量名區(qū)分大小寫(xiě)腳本文件是m文件中最簡(jiǎn)單的一種輸入頓號(hào)輸出參數(shù),用命令語(yǔ)句可以控制MATLAB命令工作空間的所有數(shù)據(jù)。編程計(jì)算向量元素的平均值:x=input('輸入向量:x='); %average_1.m計(jì)算向量元素的平均值m,n=size(x);if (m=1|n=1) average=sum(x)/length(x)else error('必須輸入向量。')End如果m文件的第一個(gè)可執(zhí)行語(yǔ)句以function開(kāi)始,該文
8、件就是函數(shù)文件,每一個(gè)函數(shù)文件都定義一個(gè)函數(shù)。函數(shù)有自己獨(dú)立的工作空間,它與MATLAB的工作空間區(qū)分開(kāi)。Function y=average_2(x) %函數(shù)average_2(x)用以計(jì)算向量元素的平均值 %輸入?yún)?shù)x為輸入向量,輸出參數(shù)y為計(jì)算的平均值MATLAB的程序控制結(jié)構(gòu)For循環(huán)一般形式For 循環(huán)控制變量=表達(dá)式1:表達(dá)式2:表達(dá)式3 For x=-2.0:0.25:-0.75While循環(huán)選擇結(jié)構(gòu):MATLAB中的選擇結(jié)構(gòu)語(yǔ)句有if語(yǔ)句、switch語(yǔ)句和try語(yǔ)句A=;for k=1:5 for j=1:5 if k=j A(k,k)=5; elseif abs(k-j)=
9、1 A(k,j)=1; else A(k,j)=0; end endendA結(jié)果:A = 5 1 0 0 0 1 5 1 0 0 0 1 5 1 0 0 0 1 5 1 0 0 0 1 5程序流的控制1. break語(yǔ)句終止本層for或while循環(huán),跳轉(zhuǎn)到本層循環(huán)語(yǔ)句end的下一條語(yǔ)句。2. return語(yǔ)句終止被調(diào)用函數(shù)的運(yùn)行,返回到調(diào)用函數(shù)。3. pause語(yǔ)句Pause:暫停程序運(yùn)行,按任意鍵繼續(xù)Pause(n):程序暫停運(yùn)行n秒后繼續(xù)。Pause on/off:允許/禁止其后的程序暫停4. continue語(yǔ)句在for循環(huán)或while循環(huán)中遇到該語(yǔ)句,將跳過(guò)其后的循環(huán)體語(yǔ)句,進(jìn)行下一
10、次循環(huán)。數(shù)據(jù)的輸入與輸出鍵盤(pán)輸入語(yǔ)句 input x=input(prompt);x=input(prompt,s) %不至于將輸入的數(shù)字看成是數(shù)值型數(shù)據(jù)。屏幕輸出語(yǔ)句disp(x)M數(shù)據(jù)文件的存儲(chǔ)/加載(save/load)格式化文本文件的存儲(chǔ)/讀?。╢printf/fxcanf)二進(jìn)制數(shù)據(jù)文件的存儲(chǔ)/讀取(fwrite/fread)數(shù)據(jù)文件行存儲(chǔ)/讀取(fgetl/fgets)MATLAB文件操作MATLAB數(shù)據(jù)可視化使用函數(shù)figure來(lái)建立圖形窗口,最簡(jiǎn)單的調(diào)用方式figureClose關(guān)閉當(dāng)前圖形窗口基本二維圖形繪制(1) plot:x軸和y軸均為線性刻度(2) Loglog:x軸和
11、y軸均為對(duì)數(shù)刻度(3) Semilogx:x為對(duì)數(shù)刻度,y為線性刻度(4) Semilogy:x為線性刻度,y為對(duì)數(shù)刻度(5) Plotyy:繪制雙縱坐標(biāo)圖形,繪制兩條具有不同縱坐標(biāo)的曲線,調(diào)用格式 Plotyy(x1,y1,x2,y2)Plot是最基本的二維繪圖函數(shù):Plot(Y),plot(X,Y)>> x=0:0.05:5;y=sin(x.2);plot(x,y);>> x=0:0.05:5;y1=0.2*x-0.8;y2=sin(x.2);plot(x,y1,x,y2);>> x=0:pi/180:2*pi; %兩條不同縱坐標(biāo)的曲線,顏色不同加以區(qū)分
12、y1=exp(-0.3*x).*cos(2*x);y2=10*exp(-1.5*x);plotyy(x,y1,x,y2);條形圖:>> x=-2.9:0.2:2.9;bar(x,exp(-x.*x);極坐標(biāo): t=0:0.1:2*pi;polar(t,abs(cos(2*t);針狀圖:>> x=0:0.1:4;y=(x.0.8).*exp(-x);stem(x,y)階梯圖:>> x=0:0.25:10;stairs(x,sin(2*x)+sin(x)餅圖:>> x=43,78,88,43,21;pie(x)色彩和線型>> x=0:0.
13、2:8;y1=0.2+sin(-2*x);y2=sin(x.0.5);plot(x,y1,'g-+',x,y2,'r-d'); %y1采用綠色、實(shí)線、加號(hào)標(biāo)記; %y2采用紅色,虛線,菱形表示坐標(biāo)軸及標(biāo)注函數(shù)axis指定坐標(biāo)軸的刻度范圍其調(diào)用格式為Axis(xmin,xmax,ymin,ymax)hold on/off 保持原有圖形/刷新原有圖形axis on/off 顯示/取消坐標(biāo)軸xlabel(option) x軸加標(biāo)注,option表示任意選項(xiàng)ylabel(option) y軸加標(biāo)注Title(option) 圖形加標(biāo)題Legend(option) 圖形加
14、標(biāo)注Grid on/off 顯示/取消網(wǎng)格線Box on/off 給坐標(biāo)加/不加邊框線>>x=0:0.05:5;y1=exp(0.4.x)-1.5;y2=sin(x*4);plot(x,y1,x,y2,'r-.');line(0,5,0,0); %畫(huà)直線,代替橫坐標(biāo)xlabel('input');ylabel('output');title('two function');legend('y1=exp(0.4x)-1.5','y2=sin(x*4)')grid on %畫(huà)網(wǎng)格線在一個(gè)圖形
15、窗口中用函數(shù)subplot可以同時(shí)畫(huà)出多個(gè)子圖形x=linspace(0,2*pi,100);subplot(2,2,1);plot(x,sin(x);xlabel('x');ylabel('y');title('sin(x)');subplot(2,2,2);plot(x,cos(x);xlabel('x');ylabel('y');title('cos(x)');subplot(2,2,3);plot(x,exp(x);xlabel('x');ylabel('y');title('exp(x)');subplot(2,2,4);plot(x,exp(-x);xlabel('x');ylabel
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 企業(yè)內(nèi)部管理中的區(qū)塊鏈技術(shù)應(yīng)用
- 從數(shù)據(jù)安全到業(yè)務(wù)協(xié)同看區(qū)塊鏈在供應(yīng)中的應(yīng)用價(jià)值
- 從行業(yè)視角看醫(yī)療數(shù)據(jù)安全及患者信息泄露的預(yù)防措施研究報(bào)告
- 醫(yī)療設(shè)備創(chuàng)新中的知識(shí)產(chǎn)權(quán)保護(hù)方法
- 醫(yī)療健康產(chǎn)業(yè)發(fā)展中的資源規(guī)劃與布局
- AI技術(shù)在醫(yī)療教育中的倫理挑戰(zhàn)與對(duì)策
- 公司出資買房合同范例
- 醫(yī)療大數(shù)據(jù)在健康管理中的價(jià)值與隱私保護(hù)的平衡
- 代理招商授權(quán)合同范例
- 亞馬遜運(yùn)營(yíng)托管合同范例
- 中外比較文學(xué)研究專題智慧樹(shù)知到期末考試答案2024年
- 心臟康復(fù)護(hù)理專家共識(shí)
- CO2氣體保護(hù)焊-基本操作方法(焊接技能)
- (高清版)TDT 1012-2016 土地整治項(xiàng)目規(guī)劃設(shè)計(jì)規(guī)范
- 河南安陽(yáng)的紅色故事
- 基于PLC的蔬菜大棚設(shè)計(jì)
- 家政保姆護(hù)工培訓(xùn)課件
- 糖尿病足圍手術(shù)期護(hù)理
- 鄉(xiāng)村振興中的鄉(xiāng)村安全與穩(wěn)定維護(hù)
- 《醫(yī)院勞動(dòng)合同書(shū)》電子版
- 營(yíng)銷策劃 -菌小寶益生菌2023品牌介紹手冊(cè)
評(píng)論
0/150
提交評(píng)論