




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)
文檔簡介
1、【精品文檔】如有侵權(quán),請聯(lián)系網(wǎng)站刪除,僅供學(xué)習(xí)與交流matlab學(xué)習(xí)筆記.精品文檔.1).>> x=-2*pi:0.1:2*pi;>> y1=sin(x);>> y2=sin(x+pi/3)+2;>> y3=cos(x);>> plot(x,y1,'.-');>> hold on>> plot(x,y2,'.-')>> hold on>> plot(x,y3,'.-')>> plot(x,y1,'.-');>
2、;> hold on>> plot(x,y2,'.-')>> hold on>> plot(x,y3,'.-')>> h=legend('sin($x$)','sin($x+fracpi3$)','cos($x$)')h = 159.0056Warning: Unable to interpret TeX string "fracpi3$)".>> set(h,'Interpreter','latex
3、9;)>> xlabel('$x$','Interpreter','latex')>> ylabel('$y$','Interpreter','latex') %添加函數(shù)圖標(biāo)及縱橫坐標(biāo)軸名稱,如下圖顯示2) .常用畫圖命令:Surf(x,y,z)、ezsurf(z)、ezmesh(z)%使用surf畫圖,要求想x,y,z是同維數(shù)的矩陣%使用ezsurf和ezmesh時,曲面函數(shù)必須能寫成顯函數(shù)的形式,如%可以使用(x,y)+函數(shù)表達(dá)式 定義匿名函數(shù),如%surf(x,y,z,
4、x); %設(shè)置圖形顏色變化方向為沿x軸,此項可省,一般不要求。圖形屬性: axis屬性 Axis tight %將坐標(biāo)軸減小到可能的最小尺寸 Axis equal %使x,y軸的刻度設(shè)置相同 Axis square %使圖形長寬相同 Axis off %不顯示坐標(biāo)軸 Axis(x0 x1 y0 y1 z0 z1); %手動改變軸的范圍 View(n,m) %設(shè)定圖形觀察視角,n為在水平方向測量的角度,m為從水平面向上測量的仰俯角。 Subplot(n,m,x) %將圖像窗口分為nxm個子窗口,當(dāng)前指的是第x個,按行數(shù)。當(dāng)要對第y個圖做一定的設(shè)置時,先輸入命令subplot(n,m,y),然后再
5、輸入要做的設(shè)置命令。 Legend(xxx) %給圖形設(shè)置圖標(biāo),也可指定圖標(biāo)的位置,默認(rèn)為右上角。 Text(x,y,xxx) %給圖像添加注釋內(nèi)容,注釋位于圖形中(x,y)處。 Gtext(xxx); %給圖像添加注釋,位置隨意,將十字光標(biāo)點在哪兒就放哪兒。 delete(findobj('type','text') %刪除通過gtext()命令添加的注釋內(nèi)容 Close all %關(guān)閉所有圖像窗口 Title(xxx) %給圖形添加標(biāo)題 Material屬性,設(shè)置曲面的材料屬性,值有dull,shiny,metallic,該命令可結(jié)合用在surf,mesh,
6、pcolor,fill,fill3等函數(shù)中。 Plotyy(x,y1,x,y2,plot) %在同一個圖中畫出兩個圖形,公用一個橫軸,縱軸不同,分列兩邊,且縱軸刻度設(shè)置可不同,自動的。 Semilogx 和semilogy 表示橫軸和縱軸采用對數(shù)坐標(biāo)畫出x-y的圖形,用法為semilogx(x,y)和semilogy(x,y). Loglog(x,y) %橫縱軸均采用對數(shù)形式畫出x-y圖形 Grid on %打開網(wǎng)格,系統(tǒng)默認(rèn)為grid off ,即關(guān)閉網(wǎng)格。 A=Complex(x,y) %產(chǎn)生復(fù)數(shù)A=x+yi Real(A)=x,imag(A)=y. Fill(x,y,n) %產(chǎn)生一個有x
7、,y坐標(biāo)定義的填充顏色的圖形,填充顏色根據(jù)n對色圖進(jìn)行索引確定。 light('Color',選項一,'Style',選項二,'Position',選項三),選項一為表示光的顏色的矩陣,取RGB三元組或相應(yīng)的顏色字符。選項二可取為'infinite'和'local'兩個值,分別表示無窮遠(yuǎn)光和近光。選項三為三維坐標(biāo)點組成的向量形式x,y,z。對遠(yuǎn)光,它表示光穿過該點射向原點;對于近光,它表示光源所在位置。假如函數(shù)不包含任何參數(shù),則采用缺省設(shè)置:白光、無窮遠(yuǎn)光、穿過(1,0,1)射向坐標(biāo)原點。一般后兩個參數(shù)可省,默認(rèn)。
8、-例: x,y,z=peaks; surf(x,y,z); shading interp; light('Color',1 0 1,'Style','local','Position',-4,-4,10); %此命令表示在點-4,-4,10處有一處品紅色光源例1.畫出橢圓柱面圖:編寫的matlab程序如下: %因為沒能寫成顯函數(shù)的形式,所以調(diào)用命令函數(shù)ezmesh時使用格式ezmesh(x,y,z),當(dāng)能寫成顯函數(shù)時可以直接調(diào)用格式ezmesh(z)。例2.繪制二元函數(shù)三維表面圖:程序如下:Meshgrid的用法舉例:兩種情況3)
9、 積分命令: a).Int-符號積分法 應(yīng)用: int(s)符號表達(dá)式s的不定積分. int(s,v)符號表達(dá)式s關(guān)于變量v的不定積分. int(s,a,b)符號表達(dá)式s的定積分, a,b分別為上下限. int(s,v,a,b)符號表達(dá)式s關(guān)于變量v從 a到b的定積分. 當(dāng)int求不出符號解,會自動轉(zhuǎn)求數(shù)值解. 注:使用符號積分時需結(jié)合使用syms-符號定義命令,例如聲明兩個符號變量 x,y,matlab語句格式為:Syms x,y; 例3.求函數(shù)不定積分: Matlab程序為: b).一重數(shù)值積分命令:quad(f,a,b)、quadl(f,a,b) 其中f為積分函數(shù),a、b為積分上下限。
10、例4.求數(shù)值積分 Matlab程序如下: c).二重數(shù)值積分命令:dblquad(f,a,b,c,d) 其中a,b為外重積分限,c,d為內(nèi)重積分限,f為積分函數(shù)。 d).三重數(shù)值積分命令:triplequad(f,a,b,c,d,e,f) 其中a,b為外重積分限,c,d為中間重積分限,e,f為內(nèi)重積分限,f為積分函數(shù)。4) .幾個特殊矩陣: eye(n):n階單位矩陣; Ones(n):n階1矩陣; Zeros(n):n階0矩陣;5) .矩陣運算: 乘法:格式 a*b; %注意和數(shù)組矩陣乘法區(qū)別a.*b;注,乘以逆和標(biāo)準(zhǔn)除法等同右乘,后除可抵消對矩陣做的左乘效果。 求逆:格式 inv(a)和運
11、算符向后除,例inv(a)*c<=>ac;注,c/a表示 指數(shù)運算:,其中k為正數(shù),分別表示A連乘K次,和? rot90(A,n):將矩陣A逆時針旋轉(zhuǎn)n個90度 tril(A):保留A的下三角并將上三角填充為0. unit8,unit16(A):將A中數(shù)據(jù)轉(zhuǎn)化為固定位數(shù)的無符號整數(shù)類型。6) 幾個命令函數(shù):find():查找命令; 例:a(find(a=0)=inf ,表示把a矩陣中的0元素全部替換為無窮。Min():查找最小命令;Numel(a),計算矩陣a的元素個數(shù);y n=Max(x,2),查找最大元素及其位置,返回值y為最大值,n為對應(yīng)位置7) Matlab與記事本的數(shù)據(jù)交
12、互: A).matlab讀入記事本文件中的數(shù)據(jù): 如果記事本文件*.txt中內(nèi)容全部為數(shù)值是,可直接調(diào)用命令A(yù)=load(*.txt),由此matlab中會生成一個數(shù)值矩陣A。 如果文件中不全為數(shù)值,則不能用load命令,應(yīng)該用name,type,valuetextread(*.txt,%s %n %f,2),此格式代表文件中第一列為字符串,第二列為整數(shù),第三列為浮點數(shù),讀取命令只讀取文件中的前兩行,生成三個列矩陣name type value,具體情況具體分析。 注:要求所讀取的txt文件的位置必須在matlab中“當(dāng)前目錄”中包含,如當(dāng)前目錄為:D:tempwork,則txt文件必須保存在
13、該路徑下,否則命令執(zhí)行出錯。 Fscanf():更強大的一個讀取命令。B).將matlab數(shù)據(jù)寫入記事本文件中: 如果只想保存數(shù)值數(shù)據(jù),則可調(diào)用save命令,格式為:save filename a b .;代表把變量a ,b .等值存入名為filename的.mat格式的文件中,該文件的路徑為“當(dāng)前目錄”。且該文件只能用matlab打開,其他軟件無法打開。 如果想把各種數(shù)據(jù)保存為.txt文件,則應(yīng)該使用fprintf命令,格式為: Fid=fopen(file.txt,wt); %打開一個文件,并返回一個指標(biāo)fid Fprintf(fid,數(shù)據(jù)表1n); %可省 Name=sally;type
14、s=1;x=3;y=4; %賦值一次 Fprintf(fid,%s type%n %f %f,name,types,x,y); %寫入一次 Name=jim;types=2;x=1;y=3; %第二次賦值 Fprintf(fid,%s type%n %f %f,name,types,x,y); %第二次寫入 Fclose(fid); %關(guān)閉文件則生成一個txt文件,內(nèi)容為: 數(shù)據(jù)表1 Sally 3 4 Jim 1 38) 數(shù)據(jù)擬合:A).多項式擬合: 指令擬合:p=polyfit(x,y,n); %返回系數(shù)矩陣p; Yi=Polyval(p,xi); %計算多項式的擬合函數(shù)值 用命令poly
15、2sym(p)輸出擬合函數(shù)表達(dá)式 圖形窗口擬合:先畫出散點圖,plot(x,y,r*);然后在圖形窗口上點擊Tools下的basic fitting,打開擬合對話框進(jìn)行操作。B).指定函數(shù)類型擬合: 假設(shè)根據(jù)散點圖知道擬合函數(shù)類型為f(x)形式的,則可利用.m進(jìn)行擬合。編寫m文件: 然后保存M文件,命名為cfun.m;然后就可以在matlab命令窗口中輸入文件名執(zhí)行文件內(nèi)容了,產(chǎn)生擬合結(jié)果。C).曲線擬合工具箱: 首先輸入要擬合的數(shù)據(jù)x和y矩陣,在matlab命令窗口中直接輸入cftool(x,y)即可調(diào)出工具箱窗口界面,然后選擇擬合類型等參數(shù)進(jìn)行擬合即可。8).matlab中的各種數(shù)學(xué)函數(shù)
16、inv 是矩陣求逆 quad 是求積分運算,并且是定積分 roots 是多項式求根 Round(x)就近取整9) .線性規(guī)劃問題: Matlab中線性規(guī)劃的標(biāo)準(zhǔn)型為: ,當(dāng)約束條件不是標(biāo)準(zhǔn)型時必須轉(zhuǎn)化為標(biāo)準(zhǔn)型。 基本函數(shù)形式為linprog(c,A,b),它的返回值是向量x的值。還有其他一些調(diào)用格式,如:,其中,fval返回目標(biāo)函數(shù)的值,Aeq和beq對應(yīng)等式約束是的上下界,為的初始值,是控制參數(shù)。例6.解下列線性規(guī)劃問題: ,可以編寫M文件,如下: 注:約束條件中有等號的藥單獨拿出來,寫成aeq和beq的形式,目標(biāo)函數(shù)是最大值時要加符號變成最?。划?dāng)目標(biāo)函數(shù)中含有常數(shù)項時要添加一個變量,并在約
17、束條件中添加一個等式。10) 圖形處理:Matlab 可識別的圖片文件格式有:tiff,png,hdf,bmp,jpg,jpeg,pcx,xwd,cur,ico等;Matlab 讀取圖像并轉(zhuǎn)化為真彩圖像、灰度圖像、索引圖像中的任一種。三種圖形介紹略。讀取命令:p=imread(filename.jpg,jpg),matlab使用這一個命令讀取所有類型的圖形文件。返回值p為一個MxNx3的包含unit8類型像素顏色值的數(shù)組,第二個參數(shù)jpg給出了明確的文件格式,該參數(shù)可選,因為matlab一般能夠從文件的內(nèi)容正確推斷出該文件的格式。顯示命令:image(p),只是該函數(shù)是對對象進(jìn)行了縮放,若要得
18、到實際尺寸的圖形,可用命令函數(shù)imshow();寫入命令:imwrite(data,newfilename.jpg,jpg),第三個參數(shù)jpg指定輸出文件的格式,data為像素數(shù)據(jù)矩陣,該命令將data數(shù)據(jù)寫入硬盤,存儲格式為jpg格式。movie2avi:動畫存儲為avi格式命令分析例子:求助,如何將下列動畫存為avi格式。set(gca,'nextplot','replacechildren');a,b,c=sphere(30);mesh(2*a,2*b,2*c,2*c);hidden off;hold on;for j=1:24x,y=meshgrid(-
19、sqrt(2):0.01:sqrt(2)*j/24);z=x.2+y.2;mesh(x,y,z);axis equal;F=getframe;endmovie(F,2);用matlab制作簡單仿真動畫,并生成.avi格式的電影文件:第一種形式:利用for循環(huán),在一定時間內(nèi)控制圖形窗口圖像的顯示,產(chǎn)生一段動態(tài)的演示過程:如下:set(gcf,'color','green'); grid on; set(gca,'zlim',-10,10');set(gca,'xlim',-10,10);set(gca,'ylim
20、39;,-10,10);for i=1:10 set(gca,'view',-i*3.75,3*i ); pause(0.2);end第二種形式:利用moviein和movie函數(shù),現(xiàn)將生成的動畫存入一個由movien函數(shù)定義的數(shù)組中,每一幀為數(shù)組的一個元素,最后用movie重復(fù)演示,movie后面的數(shù)字代表演示次數(shù)。x=-30:0.2:30;y=-30:0.2:30;x,y=meshgrid(x,y);n=5;M = moviein(n);for i=1:nz=sin(sqrt(2*(x).2+2*(y).2)-2*pi*i/10);zz=plot3(x,y,z,'p
21、arent',gca);mesh(x,y,z);grid on;colormap(0,0.9,0.5); light('position',1,1,2,'style','local','color','white');material(0.5,0.4,0.3,10,0.3);set(gca,'zlim',-10,10');M(i)=getframe(gca);endmovie(M,20)創(chuàng)建電影剪輯文件,并存儲起來,如下:aviobj=avifile('文件名.avi'
22、;,'fps',3);%定義一個avi文件, %AVIOBJ = AVIFILE(FILENAME,'PropertyName',VALUE,'PropertyName',VALUE,.)%各屬性詳細(xì)說明見matlab幫助for i=1:n%在當(dāng)前窗體上生成一幀圖像frame=getframe(gca); %獲得一幀圖像aviobj=addframe(aviobj,frame);%并加到電影剪輯文件中endaviobj=close(aviobj);%關(guān)閉文件,結(jié)束數(shù)值仿真模擬過程。11) .聲音的處理Matlab提供兩種讀取聲音文件的函數(shù):wav
23、read(filename.wav)和auread(filename.au),非別讀取wav格式的文件盒au格式的文件;且兩個函數(shù)均返回三個變量:取值(-1,1)之間的聲音向量(n行1列)-代表振幅、以HZ為單位的采樣頻率(1個常數(shù))和記錄數(shù)據(jù)所用的位數(shù)(8位或16位,分辨率不同而已。%分辨率取值一般8位或16位%)。播放命令:sound(data,rate),其中data是聲音值向量,rate為回放頻率,一般等于錄音時的頻率。該函數(shù)直接將數(shù)據(jù)傳遞為計算機的聲卡。例:b16,f16,n=wavread(filename.wav); Sound(b16,f16); b8,f8,n=wavread
24、(filename.wav); Sound(b8,f8); %可分辨8位和16位音質(zhì)的差別聲音片段剪切:例 a,b,n=wavread(filename.wav); C=a(1:2000);a(3000:end); %將文件中2000-3000的片段剪切掉并將剩余片段組合起來存入c。 C=C*0.5; %將音頻音量降低一半 Sound(a,b); %播放原文件 Sound(c,b); %播放剪切后的文件,對比效果 Plot(a); %畫出原音頻圖像 Plot(c) %畫出剪切后的文件圖形快速傅里葉變換fft(x):x為時域采樣序列,生成取值為復(fù)數(shù)的頻譜文件;ifft(x):逆變換,重新生成時域
25、序列。12) 隨機排序 如:x=1,2,3,4,3,2,1,5,編寫程序?qū)?shù)據(jù)隨機排序,程序如下 X=X(randperm(numel(x)numel(x)獲得變量長度randperm(n)產(chǎn)生長度為n的隨機數(shù)列 如:randperm(4)ans= 3 2 1 4X(.) 取得X中相應(yīng)元素產(chǎn)生新數(shù)列附錄A:matlab特殊字符、保留字和函數(shù).*數(shù)組乘法/矩陣除法./數(shù)組除法矩陣后除,等同求逆效果矩陣指數(shù)運算.數(shù)組指數(shù)運算<小于>大于<=小于等于>=大于等于=等于=不等&元素級邏輯與(向量)&&短路邏輯與(標(biāo)量)|元素級邏輯或(向量)|短路邏輯或(
26、標(biāo)量)一元非邏輯函數(shù)All(a)若邏輯向量a中的所有值都為真,則返回值為真And(a,b)若a,b都為真,則返回值為真Any(a)若邏輯向量a中有元素為真,則返回值為真Not(a)非a,即若a為真,則返回值為假,若a為假,則返回值為真Or(a,b)若a或b為真,則返回值為真文件輸入輸出loadsavecsveradfprinftimreadimwriteFclose(file)格式控制符%q由雙引號界定的引用串%e指數(shù)表示法%s字符串%f定點或十進(jìn)制表示法n換行%g定點或指數(shù)表示法t制表 b退格Ans默認(rèn)儲存變量Close all關(guān)閉所有窗口Clc清屏Exit,quit關(guān)閉matlabClea
27、r清楚變量Load載入文件Clf清空當(dāng)前圖形Who,whos查看所有變量Eps兩個浮點數(shù)之間可能出現(xiàn)的最小差Pi圓周率False邏輯假True邏輯真Inf無窮大NaN非數(shù)值基本數(shù)學(xué)函數(shù)Log(x)自然對數(shù)Exp(x)指數(shù)函數(shù)Log10(x)以10為底的對數(shù)Fix(x)將x向0方向舍入取整Mod(x,a)計算x/a的余數(shù)Floor(x)將x向負(fù)無窮舍入取整Rem(x,a)同上Round(x)將x就近取整Sqrt(x)開方Sin(x)正弦Asin(x)反正弦Cos(x)余弦Acos(x)反余弦Tan(x)正切Atan(x)反正切Atan2(y,x)給定x和y,計算反正切向量數(shù)組矩陣操作Cumsum
28、(v)計算v中元素的累積和Deal()給變量分配返回的元胞數(shù)組型結(jié)果Det(a)計算矩陣行列式的值Diag(a)抽取aj矩陣的對角元素或以向量a為對角建立矩陣Eye()生成單位陣findFliplr()左右翻轉(zhuǎn)矩陣Inv(a)矩陣求逆Length(a)獲取a的最大維數(shù)Linspace(fr,t,n)定義一個線性間隔的向量Magic(c)生成魔方矩陣v,in=max(a)查找a中最大元素的值和位置Mean(a)計算a中元素的平均值Ones(n,m)生成1矩陣Prod(x)計算x中所有元素的乘積Rand(n,m)產(chǎn)生0,1之前均勻分布的對應(yīng)維數(shù)矩陣Randn(n,m)產(chǎn)生0,1之間正態(tài)分布的對應(yīng)維
29、數(shù)矩陣Size(a)得到數(shù)組a的維數(shù)大小sparse定義稀疏矩陣v,in=sort(v)對向量v進(jìn)行排序Sum(a)計算a中元素之和Zeros(n,m)產(chǎn)生對應(yīng)維數(shù)的0數(shù)組Meshgrid(x,y)將向量x,y分別映射成二維數(shù)組2-D繪圖bar產(chǎn)生條狀圖barh產(chǎn)生水平條狀圖contour產(chǎn)生等高線圖hist產(chǎn)生直方圖loglog產(chǎn)生x,y軸均為對數(shù)尺度的圖pie產(chǎn)生餅圖plot產(chǎn)生x-y圖polar產(chǎn)生極坐標(biāo)圖semilogx僅x軸為對數(shù)尺度產(chǎn)生x-y圖semilogy近y軸為對數(shù)尺度產(chǎn)生x-y圖3-D繪圖Bar3三維條狀圖Barh3水平三維條狀圖gplot按圖論法則繪圖Mesh(x,y,z
30、)繪制曲面的網(wǎng)格圖Meshc(x,y,z)繪制帶有等高線的曲面網(wǎng)格圖Meshz(x,y,z)繪制帶有擴展垂直線的曲面網(wǎng)格圖meshgrid(n,m)產(chǎn)生三維繪圖用的網(wǎng)格peaks產(chǎn)生采樣矩陣用來要是繪圖函數(shù)Pie3繪制三維餅狀圖Plot3()繪制三維曲線圖sphere演示繪圖的例程函數(shù)surf繪制曲面圖surfc繪制帶有等高線圖的曲面圖waterfall產(chǎn)生僅在x軸具有擴展垂直線的網(wǎng)格曲面圖plotyy圖形外觀指示符- 虛線+ 加號< 朝左三角符- 實線. 圓點* 星號> 朝右三角符: 點線O 圓圈S 方塊符P 五星角符-. 點劃線X x符D 菱形符H 六星角符色圖值,與color
31、map命令搭配使用hsvspringcoolautumncopperjetsummerboneflagpinkwhitecolorcubehotprismwinterCftool擬合工具箱的用法:1.在Matlab 6.5以上的環(huán)境下,在左下方有一個"Start"按鈕,如同Windows的開始菜單,點開它,在目錄"Toolboxes"下有一個"Curve Fitting",點開"Curve Fitting Tool",出現(xiàn)數(shù)據(jù)擬合工具界面,基本上所有的數(shù)據(jù)擬合和回歸分析都可以在這里進(jìn)行。下面給你簡單介紹一下它的使用
32、方法。首先在Matlab的命令行輸入兩個向量,一個向量是你要的x坐標(biāo)的各個數(shù)據(jù),另外一個是你要的y坐標(biāo)的各個數(shù)據(jù)。輸入以后假定叫x向量與y向量,可以在workspace里面看見這兩個向量,要確保這兩個向量的元素數(shù)一致,如果不一致的話是不能在工具箱里面進(jìn)行擬合的。例如在命令行里輸入下列數(shù)據(jù):x=(0:0.02:0.98)'y=sin(4*pi*x+rand(size(x);此時x-y之間的函數(shù)近似的為正弦關(guān)系,頻率為2,但是存在一個誤差項。可以通過作圖看出它們的大體分布:plot(x,y,'*','markersize',2);打開曲線擬合共工具界面,點擊
33、最左邊的"Data."按鈕,出現(xiàn)一個Data對話框,在Data Sets頁面里,在X Data選項中選取x向量,Y Data選項中選取y向量,如果兩個向量的元素數(shù)相同,那么Create data set按鈕就激活了,此時點擊它,生成一個數(shù)據(jù)組,顯示在下方Data Sets列表框中。關(guān)閉Data對話框。此時Curve Fitting Tool窗口中顯示出這一數(shù)據(jù)組的散點分布圖。點擊Fitting.按鈕,出現(xiàn)Fitting對話框,F(xiàn)itting對話框分為兩部分,上面為Fit Editor,下面為Table of Fits,有時候窗口界面比較小,F(xiàn)it Editor部分會被收起來
34、,只要把Table of Fits上方的橫條往下拉就可以看見Fit Editor。在Fit Editor里面點擊New Fit按鈕,此時其下方的各個選框被激活,在Data Set選框中選中剛才建立的x-y數(shù)據(jù)組,然后在Type of fit選框中選取擬合或回歸類型,各個類型的擬合或回歸相應(yīng)的分別是:Custom Equations 用戶自定義函數(shù)Expotential e指數(shù)函數(shù)Fourier 傅立葉函數(shù),含有三角函數(shù)Gaussian 正態(tài)分布函數(shù),高斯函數(shù)Interpolant 插值函數(shù),含有線性函數(shù),移動平均等類型的擬合Polynomial 多項式函數(shù)Power 冪函數(shù)Rational 有
35、理函數(shù)(不太清楚,沒有怎么用過)Smooth Spline ?(光滑插值或者光滑擬合,不太清楚)Sum of sin functions正弦函數(shù)類Weibull 威布爾函數(shù)(沒用過)不好意思,沒有學(xué)過數(shù)理統(tǒng)計,所以很多東西都是用了才知道,翻譯也就不太準(zhǔn)確。不過在Type of fit選框下方有一個列表框,基本上各個函數(shù)類里的函數(shù)都寫成解析式列在下方以供選擇,所以找合適的函數(shù)還是比較容易的。在這個Type of fit選框中選擇好合適的類型,并選好合適的函數(shù)形式。于是點擊Apply按鈕,就開始進(jìn)行擬合或者回歸了。此時在Curve Fitting Tool窗口上就會出現(xiàn)一個擬合的曲線。這就是所要的
36、結(jié)果。在上面的例子中,選擇sum of sin functions中的第一個函數(shù)形式,點擊Apply按鈕,就可以看見擬合得到的正弦曲線。在Fitting對話框中的Results文本框中顯示有此次擬合的主要統(tǒng)計信息,主要有General model of sin1:. (函數(shù)形式)Coefficients (with 95% conffidence range) (95%致信區(qū)間內(nèi)的擬合常數(shù))a1=. ( . .) (等號后面是平均值,括號里是范圍).Godness of fit: (統(tǒng)計結(jié)果)SSE: . (方差)R-squared: . (決定系數(shù),不知道做什么的)Adjusted R-sq
37、uared: . (校正后的決定系數(shù),如何校正的不得而知)RMSE: . (標(biāo)準(zhǔn)差)上面的例子中經(jīng)過擬合得到的函數(shù)最后為y=0.9354*sin(12.36x+6.886)頻率為1.98加減0.03,和原來設(shè)置的頻率為2符合,相對誤差為1.5%。這是曲線擬合工具箱的一個最簡單的使用方法,上面還有很多功能,寫是寫不完的,自己參照這個基本的思路,翻著英漢詞典,看著幫助,然后一個按鈕一個按鈕的試吧。另外要說的是,如果想把這個擬合的圖像導(dǎo)出的話,在Curve Fitting Tool窗口的File菜單下選Print to Figure,此時彈出一個新的圖像窗口,里面是你要導(dǎo)出的圖像,在這個figure
38、窗口的File菜單里再選Export,選擇好合適的格式,一般是jpeg,選擇好路徑,點擊OK就可以了。出來的圖像可以在Word等編輯環(huán)境中使用,就不多說了。要修改圖像的性質(zhì),如數(shù)據(jù)點的大小、顏色等等的,只需要在對象上點右鍵,就差不多可以找到了。一些總結(jié):(1)help 命令 可以查其用法(2)“;”的用法 不加“;”顯示結(jié)果 加“;”時不顯示結(jié)果(3)數(shù)列之間的運算要加"."(4)在matlb中函數(shù)總是以離散數(shù)列形式給出的 (5)常數(shù): pi=圓周率; j=i=虛單元; ans=計算機結(jié)果的缺省變量名(6)數(shù)列:a=1 5 8 0 2矩陣:a=1 5 8;0 8 4改變數(shù)列
39、的某一個項的值:a(1)=11,則a=11 5 8 0 2/注意不是a(0) 一:matlb的序列的表示方法1:枚舉法 eg a=0 1 2 3 4 5 6 7 8 92:已知首相尾相,與等差 eg a=0:1:9 或者簡寫成 a="0:9"(等差1可省略)3:已知首相尾相,與序列總數(shù) eg a=linspace(0,9,10)應(yīng)用:x軸可寫為 x=0:n-1 二:畫圖1:畫出離散圖形:stem eg stem(y) 注意,此時的x由系統(tǒng)自定。 stem(x,y) 注意,x與y數(shù)列的點數(shù)一定要相同。2:畫出連續(xù)圖形:plot eg plot(x,y) plot(y)3:把多
40、個圖放在一塊:subplot(子圖) eg subplot(3,1,2)/一共三行一列,此圖在第2個 eg subplot(3,2,5)/一共三行兩列,此圖在第5個4:指定圖形標(biāo)題欄:title('') eg title('頻譜圖') 注意,要加'' 指定x(自變量)軸標(biāo)注文字:xlabel('') 指定y(自變量)軸標(biāo)注文字:ylabel('')5:坐標(biāo)軸的范圍的刻畫:axisx1,x2,y1,y2 以坐標(biāo)紙方式顯示: grid應(yīng)用: eg x=0:9; 定義x軸 y=sin(2*pi/30*x)./x; 定義y
41、軸 subplot(3,1,1); 在適當(dāng)位置處畫出一個空白的圖像 plot(x,y); 填充此圖像 title('一個例子'); 為此圖像加標(biāo)識6:把多個圖形畫于同一個坐標(biāo)中 eg plot(t, y, t, y0, t, -y0) plot(t, y, '-r', t, y0, ':b', t, -y0, ':b')/'-r'表示紅色線;':b'表示虛藍(lán)線 如果后續(xù)圖形也要畫在當(dāng)前的同一圖形窗內(nèi),并且希望二條線同時存在,可輸入:hold on;hold off關(guān)閉。7:不同的線的畫法具體的直角坐
42、標(biāo)圖繪圖格式:plot(x,y,'cm')式中x為自變量,y為因變量,c為顏色標(biāo)記,m為線型標(biāo)記符 b 藍(lán)色 g 綠色 r 紅色 c 青色 m 洋紅 y 黃色 k 黑色 w 白色 . 小黑點(線型)o 小圈號 x 叉號 + 十字號 * 星號 s 方形 d 菱形 v 三角形 三角形(向上) < 三角形(向左) > 三角形(向右) p 五角星 h 六角星 - 實連線 : 虛點連線 -. 點劃連線 - 雙劃連線注意:先寫subplot,再寫plot,最后為title axis等figure窗口不消失:grid on; 三:提取性質(zhì)1:提取出數(shù)列的長度:length(數(shù)列)
43、 四:運算數(shù)列之間的運算要加"."1:四則運算:+ ;- ;* ;/ ;以及數(shù)列之間的運算:.+ ;.- ;.* ;./;注意當(dāng)分母為一個數(shù)列時,只能用./不能用/取余運算:mod(a,b)Eg1mod(13,5)ans = 3Eg2mod(1:5,3)ans = 1 2 0 1 2取整:Fix(a)2:三角函數(shù):sin();cos();tan();cot();3:幕運算:以及數(shù)列之間的運算:. 指數(shù)運算:exp()應(yīng)用:eg n="0:50" x="-25:25" W=(pi/12.5)*k; /w=(2*pi/25)*k y=(e
44、xp(-j*pi/12.5).(n.*x); eg x="-25:25" w=(pi/12.5); y="exp"(-j*w*x);/所得為一系列復(fù)數(shù)注意有數(shù)列時則一定用.不用 eg:a=0:9; b="a".5; c="5".a;(自己記:符號兩邊不對稱的,一旦有一個數(shù)列,則用.)4:復(fù)數(shù)相關(guān)abs()復(fù)數(shù)大小(幅頻特性)angle()復(fù)角(相頻特性)real()實部imag()虛部 5:累加一維數(shù)組的話直接sum(b)要疊加行或列用sum(b,1)或sum(b,2)二維數(shù)組要疊加所有元素用sum(b(:) 五:
45、語句1: ifelse語句 if 條件 運行語句1; 運行語句2; 運行語句n; else/可以省略 運行語句1; 運行語句2; 運行語句n; end 也可以 if else if else end end2: for循環(huán)語句 for 數(shù)列/多用0:n-1去執(zhí)行n次;或應(yīng)用length取長度 運行語句1; 運行語句2; 運行語句n; End3:while 邏輯 語句1; 語句2;End 六:信號的產(chǎn)生1:0矩陣的創(chuàng)立:zeros(m行,n列)/行是一的矩陣也就是一個序列2:1矩陣的創(chuàng)立:ones(m行,n列)/行是一的矩陣也就是一個序列應(yīng)用:單位抽樣序列:1 zeros(1,n-1) =1 0
46、 0 0 0 0 0 0.(一共n個數(shù)字) 單位階躍序列:ones(1,n)=1 1 1 1 1 1 1 . 七:編寫*.m文件,自定義的函數(shù)過程:書寫方法:定義為*.m文件。多要利用循環(huán),分支語句。運行方法:進(jìn)入含有自定義函數(shù)的*.m文件。在命令行里鍵入:m文件的名字,回車。Eg: clear all;close all;clc;cd lib函數(shù):書寫方法:和過程相似,但要有:function 返回值=函數(shù)名(參數(shù)1,參數(shù)2,參數(shù)3。)主意:可以不定義返回值,參數(shù)n。m文件名=函數(shù)名局部變量:全局變量:利用global命令Eg:global num;num=1000;運行方法: 與一般的函數(shù)一樣,返回值,參數(shù)n在函數(shù)外部定義。 Eg: function re="text"(f,n)w=2*pi*f;t=0:n-1;re=sin(w*t);plot(re);多值返回ra,rb,rc,.=或者在m編譯器中debug>>run
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 組織變革與發(fā)展核心內(nèi)容
- 農(nóng)殘檢測培訓(xùn)
- 酒店培訓(xùn)后的心得匯報
- 腫瘤患者居家監(jiān)護(hù)制度
- 痿證的中醫(yī)診療體系
- 中層管理人員培訓(xùn)開發(fā)計劃
- 培訓(xùn)管理與開發(fā)
- 護(hù)理保健基本知識
- 聲音親切度培訓(xùn)
- 急性支氣管哮喘護(hù)理查房
- 運營崗位轉(zhuǎn)正述職報告
- 員工離職協(xié)議書模板
- 小講課闌尾炎病人的護(hù)理護(hù)理病歷臨床病案
- 2024年電阻陶瓷基體項目可行性研究報告
- 大國兵器學(xué)習(xí)通超星期末考試答案章節(jié)答案2024年
- 24秋國家開放大學(xué)《馬克思主義基本原理》專題測試參考答案
- 苗木采購施工方案
- 【大單元】1《北京的春節(jié)》第2課時 教學(xué)設(shè)計
- 教育現(xiàn)象及問題分析
- 2024年新華東師大版七年級上冊數(shù)學(xué)全冊教案(新版教材)
- 2024年一級健康管理師考前沖刺必會試題庫300題(含詳解)
評論
0/150
提交評論