電氣類畢業(yè)論文-----基于matlab的逆變電源研究

1、<p><b>　　目錄</b></p><p><b>　　1 方案論證2</b></p><p>　　1.1 設計實現(xiàn)要求2</p><p>　　1.2 設計方案確定2</p><p><b>　　2 原理簡介3</b></p><p&

2、gt;　　2.1 升壓斬波電路3</p><p>　　2.1.1 升壓電路原理圖3</p><p>　　2.1.2 原理分析3</p><p>　　2.2 三相電壓型橋式逆變電路4</p><p>　　2.2.1 逆變電路原理圖4</p><p>　　2.2.2 逆變電路原理4</p><

3、;p>　　2.3 SPWM逆變器的工作原理5</p><p>　　2.4 Simulink仿真環(huán)境5</p><p><b>　　3 仿真建模5</b></p><p>　　3.1 斬波電路仿真建模6</p><p>　　3.2 逆變電路仿真建模7</p><p>　　3.3 逆變

4、電源仿真建模8</p><p><b>　　4 仿真實現(xiàn)9</b></p><p>　　4.1 斬波電路仿真實現(xiàn)9</p><p>　　4.2 逆變電路仿真實現(xiàn)10</p><p>　　4.3 逆變電源仿真實現(xiàn)10</p><p><b>　　5 心得體會12</b&g

5、t;</p><p><b>　　參考文獻12</b></p><p>　　附錄 </p><p><b>　　逆變電源研究</b></p><p><b>　　1 方案論證</b

6、></p><p>　　1.1 設計實現(xiàn)要求</p><p>　　本次設計要求對逆變電源進行Matlab仿真研究，輸入為100V，輸出為380V、50Hz三相交流電，采用PWM斬波控制技術，建立Matlab仿真模型并得到實驗結(jié)果。</p><p>　　1.2 設計方案確定</p><p>　　由于要求的輸出為380V、50Hz三相交流電

7、，顯然不能直接由輸入的100V直流電逆變產(chǎn)生，需將輸入的100V直流電壓通過升壓斬波電路提高電壓，再經(jīng)過逆變過程及濾波電路得到要求的輸出。根據(jù)課本所學的，可以采用升壓斬波電路和三相電壓型橋式逆變電路的組合電路，將升壓后的電壓作為逆變電路的直流側(cè)，得到三相交流電，同時采用PWM控制技術，使其頻率為50HZ。</p><p>　　根據(jù)對輸出電壓平均值進行調(diào)制的方式不同，斬波電路分為脈沖寬度調(diào)制（PWM）、頻率調(diào)制和混

8、合型三種。在此使用第一種方法，也是應用最多的方法。通過控制開關器件的通斷實現(xiàn)電能的儲存和釋放過程，輸出信號為方波，調(diào)節(jié)脈寬可以控制輸出的電壓的大小。</p><p>　　根據(jù)直流側(cè)電源性質(zhì)不同，逆變電路可分為電壓型逆變電路和電流型逆變電路。這里的逆變電路屬電壓型。采用等腰三角波作為載波，用SPWM進行雙極性控制。該電路的輸出含有諧波，除了使波形具有對稱性減少諧波和簡化控制外，還需要專門的濾波電路進行濾波。濾波電路

9、采用RLC濾波電路。</p><p>　　直流斬波電路采用PWM斬波控制，輸出的方波經(jīng)過濾波電路后變?yōu)橹绷麟娝屯孀冸娐?。逆變采用PWM逆變電路，采用SPWM作為調(diào)制信號，輸出PWM波形，再經(jīng)過濾波電路得到380V、50Hz三相交流電。系統(tǒng)總體框圖如圖1-1所示。</p><p>　　圖1-1 系統(tǒng)總體框圖</p><p>　　最后通過MATLAB集成環(huán)境simul

10、ink進行仿真，驗證結(jié)果。</p><p><b>　　2 原理簡介</b></p><p>　　2.1 升壓斬波電路</p><p>　　2.1.1 升壓電路原理圖</p><p>　　圖2-1 升壓斬波電路原理圖</p><p>　　2.1.2 原理分析</p><p>

11、;　　假設L值、C值很大，V通時，E向L充電，充電電流恒為I1，同時C的電壓向負載供電，因C值很大，輸出電壓uo為恒值，記為Uo。設V通的時間為ton，此階段L上積蓄的能量為EI1ton，V斷時，E和L共同向C充電并向負載R供電。設V斷的時間為toff，則此期間電感L釋放能量為 ，穩(wěn)態(tài)時，一個周期T中L積蓄能量與釋放能量相等，</p><p>　　,化簡得。輸出電壓高于電源電壓，故稱升壓斬波電路，也稱之為boos

12、t變換器。T與的比值為升壓比，將升壓比的倒數(shù)記作β，由于α+β=1，故。升壓斬波電路能使輸出電壓高于電源電壓的原因 ：L儲能之后具有使電壓泵升的作用，并且電容C可將輸出電壓保持住。</p><p>　　2.2 三相電壓型橋式逆變電路</p><p>　　2.2.1 逆變電路原理圖</p><p>　　圖2-2 三相橋式逆變電路原理圖</p><p

13、>　　2.2.2 逆變電路原理</p><p>　　該電路采用雙極性控制方式，U、V和W三相的PWM控制通常公用一個三角載波，三相的調(diào)制信號、和一次相差120°。U、V和W各相功率開關器件的控制規(guī)律相同，現(xiàn)以U相為例來說明。當>時，給上橋臂以導通信號，給下橋臂以關斷信號，則U相相對于直流電源假想中點的輸出電壓。當<時，給以導通信號，給以關斷信號，則。和的驅(qū)動信號始終是互補的。當給（）加

14、導通信號時，可能是（）導通，也可能是二極管（）續(xù)流導通，這要由阻感負載中電流的方向來決定。V相和W相的控制方式都和U相相同。</p><p>　　2.3 SPWM逆變器的工作原理</p><p>　　由于期望的逆變器輸出是一個正弦電壓波形，可以把一個正弦半波分作N等分。然后把每一等分的正弦曲線與橫軸所包圍的面積都用個與此面積相等的等高矩形脈沖來代替，矩形脈沖的中點與正弦波每一等分的中點重合

15、。這樣，由N個等幅不等寬的矩形脈沖所組成的波形為正弦的半周等效。同樣，正弦波的負半周也可用相同的方法來等效。</p><p>　　這一系列脈沖波形就是所期望的逆變器輸出SPWM波形。由于各脈沖的幅值相等，所以逆變器可由恒定的直流電源供電，也就是說，這種交一直一交變頻器中的整流器采用不可控的二極管整流器就可以了。逆變器輸出脈沖的幅值就是整流器的輸出電壓。當逆變器各開關器件都是在理想狀態(tài)下工作時，驅(qū)動相應開關器件的信

16、號也應為與形狀相似的一系列脈沖波形，這是很容易推斷出來的。</p><p>　　2.4 Simulink仿真環(huán)境</p><p>　　Simulink是Matlab的仿真集成環(huán)境，是一個實現(xiàn)動態(tài)系統(tǒng)建模、仿真的集成環(huán)境。它使Matlab的功能進一步增強，主要表現(xiàn)為：①模型的可視化。在Windows環(huán)境下，用戶通過鼠標就可以完成模型的建立與仿真；②實現(xiàn)了多工作環(huán)境間文件互用和數(shù)據(jù)交換；③把理

17、論和工程有機結(jié)合在一起。利用Matlab下的Simulink軟件和電力系統(tǒng)模塊庫(SimPowerSystems)進行系統(tǒng)仿真是十分簡單和直觀的，用戶可以用圖形化的方法直接建立起仿真系統(tǒng)的模型，并通過Simulink環(huán)境中的菜單直接啟動系統(tǒng)的仿真過程，同時將結(jié)果在示波器上顯示出來。本文主要通過對逆變電源的Matlab仿真，研究逆變電路的輸入輸出及其特性，以及一些參數(shù)的選擇設置方法，從而為以后的學習和研究奠定基礎，同時也學習使用Matla

18、b軟件的Simulink集成環(huán)境進行仿真的相關操作。</p><p><b>　　3 仿真建模</b></p><p>　　根據(jù)系統(tǒng)總體框圖，可將其分為PWM升壓斬波電路和三相逆變電路電路（含濾波電路），下面分別對其進行仿真建模。</p><p>　　3.1 斬波電路仿真建模</p><p>　　斬波電路我采用了升壓斬波

19、電路，電路原理圖如3-1所示，原理前面也講得很清楚了。電路輸出的電壓還要經(jīng)逆變后濾波，故對波形的要求不是很高，與負載并聯(lián)的電容C取很大，就可以達到濾波的目的，因此不需另外添加濾波電路。</p><p>　　圖3-1 升壓斬波電路</p><p>　　該電路中開關器件用IGBT，控制IGBT的波形由PWM脈沖生成器Pulse Generator產(chǎn)生，Pulse Generator在Simul

20、ink Library Browser的Simulink下拉菜單Sources類別中。繪制仿真圖時，打開Simulink Library Browser，可以在分類菜單中查找所需元件，也可以直接在查找欄中輸入元件名稱，如Pulse Generator，雙擊查找。找到元件后直接將其拖到新建Model文件窗口中即可。電路中其他元件按以上方法找出，放入Model文件窗口中。其中電阻、電感和電容元件，選擇SimPowerSystems下拉菜單E

21、lements類別中的Series RLC Branch，放入窗口后，雙擊該圖標，在Branch Type中選擇相應類型，如電阻選R，電感選L，選擇完畢后單擊OK按鈕。放齊元件后，按升降壓斬波電路原理圖連接電路，為了方便觀察輸出，應在輸出端加上電壓測量裝置Voltage Measurement，并在Simulink下拉菜單Commonly Used Blocks類別中選擇Scope，即示波器，以觀測輸出電壓波形。所</p>

22、<p>　　圖3-2 升壓斬波電路仿真模型</p><p>　　3.2 逆變電路仿真建模</p><p>　　三相橋式PWM型逆變電路的模型可參考Marlab—>Help—>Demos—>Simulink—>SimPowerSystems—>General Demos中的Three-Phase Two-Level PWM Voltage Sourc

23、e Converters。此電路采用了三相逆變橋集成塊Universal Bridge 3 arms，濾波電路也已由Three-Phasse Parallel RLC Load模塊構(gòu)成，不需另加濾波電路。在此電路的基礎上稍作修改，即構(gòu)成三相橋式PWM型逆變電路模型，如圖3-3所示。其中變壓器僅起隔離作用，不對電壓進行升降。</p><p>　　圖3-3 三相橋式PWM型逆變電路模型</p><

24、p>　　3.3 逆變電源仿真建模</p><p>　　將斬波電路的輸出接到逆變電路的輸入，就得到逆變電源仿真模型，如圖3-4所示。</p><p>　　圖3-4 逆變電源仿真模型</p><p><b>　　4 仿真實現(xiàn)</b></p><p>　　4.1 斬波電路仿真實現(xiàn)</p><p>

25、　　打開斬波電路窗口，根據(jù)參考資料設置初試參數(shù)，設置時雙擊元件圖標。輸入直流電設為100V，開關器件IGBT和二極管Diode使用默認參數(shù)。負載R=50 ，電感L=6e-04H（即0.6mH），電容C=3e-05F（即30uF）。設置PWM發(fā)生器周期Period為0.0001s，占空比Pulse Width(% of period)為84.1%，其他參數(shù)不變。單擊Start simulation按鍵，開始仿真，雙擊示波器Scope，觀察

26、輸出波形圖。此時輸出波形持續(xù)等副震蕩，且幅值太高，很不理想。分析知升降壓斬波電路中電感和電容值均應很大，將電容值改為500uF（C=5e-04F），電感值為3.8mH，觀察波形，如圖4-1 所示，輸出電壓約0.2s后穩(wěn)定在630V。</p><p>　　圖4-1逆變電路仿真波形</p><p>　　通過幾次調(diào)節(jié)各元件參數(shù)發(fā)現(xiàn)，改變電感和電容的值，輸出電壓穩(wěn)定值也在變化。電容的作用主要是使輸

27、出電壓保持住，電容值過小輸出波形會持續(xù)震蕩，應取較大，但過大的電容值會使輸出電壓穩(wěn)定的時間太長。根據(jù)以上規(guī)律反復改變各元件參數(shù)，直到得到滿意的結(jié)果。</p><p>　　4.2 逆變電路仿真實現(xiàn)</p><p>　　逆變電路中的參數(shù)主要有離散PWM生成器（Discrete PWM Generator）中的載波頻率（Carrier frequency）、采樣時間（Sample time）、調(diào)

28、制參數(shù)（Modulation index）和輸出電壓頻率（Frequency of output voltage），變壓器（Transformer）中的繞組參數(shù)（Winding parameters），Three-Phasse Parallel RLC Load模塊中的參考點相電壓（Nominal phase-to-phase voltage）和頻率（Nominal frequency）。由于要求輸出交流電頻率為50Hz，上述參數(shù)中的頻

29、率都應為50Hz。載波頻率過小時，輸出的波形有些跳動，和正弦波有一定的差距，如圖4-2 逆變電路仿真波形所示。</p><p>　　圖4-2逆變電路仿真波形</p><p>　　4.3 逆變電源仿真實現(xiàn)</p><p>　　首先應將斬波電路的輸出電壓調(diào)到630V左右，再對逆變電源進行仿真。反復調(diào)節(jié)參數(shù)知當斬波電路中PWM脈沖生成器的占空比達到84.1%時，輸出的直流

30、電壓約為630V，此時的波形如圖4-5 所示，輸出電壓先大幅震蕩，大約0.2s后，穩(wěn)定在630V左右。</p><p>　　圖4-3逆變電源斬波輸出波形</p><p>　　圖4-4 逆變電源輸出波形</p><p>　　改變逆變電源仿真模型中的參數(shù)到要求值，單擊Start simulation按鍵開始仿真，圖4-4為逆變電源輸出波形。有圖可知，逆變電源輸出三相交流

31、電相電壓波形幅值為311V，各相電壓互差120°，周期為0.02s即頻率為50Hz，滿足任務要求。</p><p><b>　　5 心得體會</b></p><p>　　通過本次電力電子技術課程設計，我學會了使用Matlab軟件仿真集成環(huán)境Simulink進行仿真的基本操作方法，熟悉了用protel畫電路圖，也對斬波電路、逆變電源的原理有了進一步的理解。電力

32、電子技術是我們非常喜歡的一門課程，學到的知識非常實用，這次課設讓我把書本上的理論知識和實際電路仿真聯(lián)系起來，在增強實踐能力的同時，我也學到了運用matlab仿真電力電子電路，以后有很多電路我都可以做出仿真來。</p><p>　　在做此次課程設計中，我還是遇到了很多問題的，比如畫圖時找不到IGBT元件，設定電感電容參數(shù)時花了好多時間，仿真時還找不到所需的元件，最后還是通過問同學和上網(wǎng)查資料才完成，還真有點辛苦。M

33、atlab很強大，以后還要繼續(xù)研究??！</p><p><b>　　參考文獻</b></p><p>　　1]王兆安等.電力電子技術.北京：機械工業(yè)出版社，2007</p><p>　　[2]薛定宇.基于Matlab/Simulink的系統(tǒng)仿真技術與應用.北京：清華大學出版社，2006</p><p>　　[3]王丹力等

34、.Matlab控制系統(tǒng)設計、仿真、應用.北京:中國電力出版社,2007</p><p>　　[4]周建興等. MATLAB從入門到精通.北京：人民郵電出版社，2008</p><p>　　[5]陳國呈.PWM逆變技術及應用.北京：中國電力出版社,2007</p><p><b>　　附錄一：</b></p><p>&l

35、t;b>　　程序代碼（1）：</b></p><p>　　clear all;</p><p>　　t=-2*pi:0.01:2*pi;</p><p>　　y=3*sin(2*t);</p><p>　　plot(t,y);</p><p>　　xlabel ('時間');ylabe

36、l('信號強度');</p><p>　　title('正弦信號');</p><p><b>　　grid on;</b></p><p>　　axis([-7 7 -3.5 3.5]);</p><p><b>　　實驗結(jié)果：</b></p><

37、;p><b>　　正弦序列：</b></p><p>　　clear all;</p><p><b>　　n=-20:20;</b></p><p>　　y=sin(0.1*pi*n);</p><p>　　stem(n,y);</p><p>　　xlabel (&

38、#39;n');ylabel('信號強度');</p><p>　　title('正弦信號序列');</p><p>　　axis([-20 20 -3.5 3.5]);</p><p><b>　　grid on;</b></p><p><b>　　輸出結(jié)果：<

39、/b></p><p><b>　　程序代碼（2）：</b></p><p>　　fs=100;%設定采樣頻率</p><p><b>　　N=256;</b></p><p><b>　　n=0:N-1;</b></p><p><b>

40、;　　t=n/fs;</b></p><p>　　x=sinc(2*pi*t-8); %生成sa信號</p><p>　　subplot(311);</p><p>　　plot(t,x);%作sa信號的時域波形</p><p>　　xlabel('t');</p><p>　　ylabel

41、('y');</p><p>　　title('sa信號sinc(2*pi*t-8)時域波形');</p><p><b>　　grid on;</b></p><p>　　y=fft(x,N);%進行N點fft變換</p><p>　　mag=abs(y);%求幅值</p>

42、<p>　　f=(0:length(y)-1)'*fs/length(y);%進行對應的頻率轉(zhuǎn)換</p><p>　　subplot(312);</p><p>　　plot(f,mag);%做頻譜圖</p><p>　　axis([0,100,0,20]);</p><p>　　xlabel('頻率(Hz)

43、9;);</p><p>　　ylabel('幅值');</p><p>　　title('sa信號sinc(2*pi*t-8)時域波形幅頻譜圖N=256');</p><p><b>　　grid on;</b></p><p>　　z=angle(y);</p><

44、p>　　subplot(313);</p><p>　　plot(f,z);%做相譜圖</p><p>　　axis([0,100,-5,5]);</p><p>　　xlabel('頻率(Hz)');</p><p>　　ylabel('角度');</p><p>　　title

45、('sa信號sinc(2*pi*t-8)相頻圖N=256');</p><p><b>　　grid on;</b></p><p><b>　　波形：</b></p><p><b>　　程序代碼（3）：</b></p><p>　　function vara

46、rgout = two(varargin)</p><p>　　% Begin initialization code - DO NOT EDIT</p><p>　　gui_Singleton = 1;</p><p>　　gui_State = struct('gui_Name', mfilename, ...</p>

47、<p>　　'gui_Singleton', gui_Singleton, ...</p><p>　　'gui_OpeningFcn', @two_OpeningFcn, ...</p><p>　　'gui_OutputFcn', @two_OutputFcn, ...</p><p>　　'

48、gui_LayoutFcn', [] , ...</p><p>　　'gui_Callback', []);</p><p>　　if nargin && ischar(varargin{1})</p><p>　　gui_State.gui_Callback = str2func(varargin{1});</

49、p><p><b>　　end</b></p><p>　　if nargout</p><p>　　[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});</p><p><b>　　else</b></p><p

50、>　　gui_mainfcn(gui_State, varargin{:});</p><p><b>　　end</b></p><p>　　% End initialization code - DO NOT EDIT</p><p>　　% --- Executes just before erzhi is made visible

51、.</p><p>　　function two_OpeningFcn(hObject, eventdata, handles, varargin)</p><p>　　% This function has no output args, see OutputFcn.</p><p>　　--------------------------------------

52、-------------------</p><p>　　%接下來是菜單的程序</p><p>　　---------------------------------------------------------</p><p>　　function file_Callback(hObject, eventdata, handles)</p>&l

53、t;p>　　--------------------------------------------------------</p><p><b>　　%文件打開程序</b></p><p>　　function open_Callback(hObject, eventdata, handles)</p><p>　　[name,pat

54、h]=uigetfile({'*.bmp'},'載入圖像');</p><p>　　if isequal(name,0)|isequal(path,0)</p><p>　　errordlg('沒有選中文件','出錯');</p><p><b>　　return;</b><

55、/p><p><b>　　else</b></p><p>　　x=imread([path,name]);</p><p>　　axes(handles.axes1);</p><p>　　imshow(x);</p><p>　　handles.img=x;</p><p>

56、;　　handles.noise_img=x;</p><p>　　guidata(hObject,handles)</p><p><b>　　end</b></p><p>　　---------------------------------------------------</p><p>　　function

57、 save_Callback(hObject, eventdata, handles)</p><p><b>　　%文件保存</b></p><p>　　[filename,pathname] = uiputfile('*.bmp','圖片保存為');</p><p>　　if isequal([filenam

58、e,pathname],[0,0])</p><p>　　errordlg('沒有保存','出錯');</p><p><b>　　return;</b></p><p><b>　　else</b></p><p>　　file=strcat(pathname,fi

59、lename);</p><p>　　(handles.axes2);</p><p>　　i=getimage(gca);</p><p>　　imwrite(i,file);</p><p><b>　　end</b></p><p>　　---------------------------

60、--------------------------</p><p>　　function exit_Callback(hObject, eventdata, handles)</p><p><b>　　clc;</b></p><p>　　close all;</p><p>　　close(gcf);</p&

61、gt;<p>　　---------------------------------------------------------</p><p><b>　　%編輯菜單的程序</b></p><p>　　------------------------------------------------------</p><p>

62、;　　function edit_Callback(hObject, eventdata, handles)</p><p>　　-------------------------------------------------</p><p>　　function huidu_Callback(hObject, eventdata, handles)</p><p&g

63、t;<b>　　%灰度處理</b></p><p>　　axes(handles.axes2);</p><p>　　if isrgb(handles.img) </p><p>　　y=rgb2gray(handles.img); %RGB</p><p>　　imshow(y);</p>

64、;<p><b>　　else</b></p><p>　　msgbox('這已經(jīng)是灰度圖像','轉(zhuǎn)換失敗');</p><p><b>　　end </b></p><p>　　------------------------------------------------

65、------</p><p>　　function liangdu_Callback(hObject, eventdata, handles)</p><p><b>　　%亮度處理</b></p><p>　　prompt={'輸入?yún)?shù)1','輸入?yún)?shù)2','輸入gamma'};</p&g

66、t;<p>　　defans={'[0 0.7]','[0 1]','1'};</p><p>　　p=inputdlg(prompt,'輸入?yún)?shù)',1,defans);</p><p>　　p1=str2num(p{1});</p><p>　　p2=str2num(p{2});<

67、/p><p>　　p3=str2num(p{3});</p><p><b>　　gamma=p3;</b></p><p>　　x=(handles.img);</p><p>　　y=imadjust(x,p1,p2,gamma);</p><p>　　axes(handles.axes2);&l

68、t;/p><p>　　imshow(y);</p><p>　　------------------------------------------------------</p><p>　　function jietu_Callback(hObject, eventdata, handles)</p><p><b>　　%截圖&l

69、t;/b></p><p>　　set(handles.axes2,'HandleVisibility','ON');</p><p>　　axes(handles.axes2);</p><p>　　y=imcrop(handles.img);</p><p>　　imshow(y); </p&g

70、t;<p>　　handles.Timage=y;</p><p>　　------------------------------------------------------</p><p>　　function fangda_Callback(hObject, eventdata, handles)</p><p>　　------------

71、------------------------------------------</p><p>　　function lingjz_Callback(hObject, eventdata, handles)</p><p><b>　　%鄰近插值放大</b></p><p>　　axes(handles.axes2);</p>

72、<p>　　prompt={'輸入放大倍數(shù):'};</p><p>　　defans={'2'};</p><p>　　p=inputdlg(prompt,'輸入放大倍數(shù)',1,defans);</p><p>　　p1=str2num(p{1});</p><p>　　y=imr

73、esize(handles.img,p1,'nearest'); %最近鄰插值法放大</p><p>　　imshow(y);</p><p>　　------------------------------------------------------</p><p>　　function shuangxianfa_Callback

74、(hObject, eventdata, handles)</p><p><b>　　%雙線放大</b></p><p>　　axes(handles.axes2);</p><p>　　prompt={'輸入放大倍數(shù):'};</p><p>　　defans={'2'};</p&

75、gt;<p>　　p=inputdlg(prompt,'輸入放大倍數(shù)',1,defans);</p><p>　　p1=str2num(p{1});</p><p>　　y=imresize(handles.img,p1,'bilinear'); %最近鄰插值法放大</p><p>　　imshow

76、(y);</p><p>　　------------------------------------------------</p><p>　　function shuangsanci_Callback(hObject, eventdata, handles)</p><p><b>　　%雙三次放大</b></p><

77、p>　　axes(handles.axes2);</p><p>　　prompt={'輸入放大倍數(shù):'};</p><p>　　defans={'2'};</p><p>　　p=inputdlg(prompt,'輸入放大倍數(shù)',1,defans);</p><p>　　p1=str2n

78、um(p{1});</p><p>　　y=imresize(handles.img,p1,'bicubic'); </p><p>　　imshow(y);</p><p>　　-------------------------------------------------------</p><p>

79、;　　function suoxiao_Callback(hObject, eventdata, handles)</p><p><b>　　%鄰近插值縮小</b></p><p>　　axes(handles.axes2);</p><p>　　prompt={'輸入放大倍數(shù):'};</p><p>

80、　　defans={'0.2'};</p><p>　　p=inputdlg(prompt,'輸入放大倍數(shù)',1,defans);</p><p>　　p1=str2num(p{1});</p><p>　　y=imresize(handles.img,p1,'nearest'); %最近鄰插值法縮小&

81、lt;/p><p>　　imshow(y);</p><p>　　------------------------------------------------------</p><p>　　function shuangxian_Callback(hObject, eventdata, handles)</p><p><b>　

82、　%雙線性縮小</b></p><p>　　axes(handles.axes2);</p><p>　　prompt={'輸入放大倍數(shù):'};</p><p>　　defans={'0.2'};</p><p>　　p=inputdlg(prompt,'輸入放大倍數(shù)',1,defa

83、ns);</p><p>　　p1=str2num(p{1});</p><p>　　y=imresize(handles.img,p1,'bilinear'); %最近鄰插值法縮小</p><p>　　imshow(y); ------------------------------------------------------

84、--</p><p>　　function shuangsancisuox_Callback(hObject, eventdata, handles)</p><p><b>　　%雙三次縮小</b></p><p>　　axes(handles.axes2);</p><p>　　prompt={'輸入放大倍數(shù)

85、:'};</p><p>　　defans={'0.2'};</p><p>　　p=inputdlg(prompt,'輸入放大倍數(shù)',1,defans);</p><p>　　p1=str2num(p{1});</p><p>　　y=imresize(handles.img,p1,'bicu

86、bic'); </p><p>　　imshow(y); </p><p>　　------------------------------------------------------</p><p>　　function xuanzhuan_Callback(hObject, eventdata, handles)</p>

87、<p>　　-------------------------------------------------------</p><p>　　function updown_Callback(hObject, eventdata, handles)</p><p><b>　　%上下翻轉(zhuǎn)</b></p><p>　　axes(h

88、andles.axes2);</p><p>　　x=(handles.img);</p><p>　　if isrgb(handles.img) </p><p>　　for k=1:3 </p><p>　　y(:,:,k)=flipud(x(:,:,k));%上下翻轉(zhuǎn)函數(shù)</p><p><b>　　e

89、nd</b></p><p>　　imshow(y);</p><p><b>　　else</b></p><p>　　x=(handles.img);</p><p>　　y=flipud(x); </p><p>　　imshow(y);</p><p>

90、<b>　　end</b></p><p>　　---------------------------------------------------------</p><p>　　function leftright_Callback(hObject, eventdata, handles)</p><p><b>　　%左右翻轉(zhuǎn)

91、</b></p><p>　　axes(handles.axes2);</p><p>　　if isrgb(handles.img) </p><p>　　x=(handles.img);</p><p>　　for k=1:3 </p><p>　　y(:,:,k)=fliplr(x(:,:,k));%

92、左右翻轉(zhuǎn)函數(shù)</p><p><b>　　end</b></p><p>　　imshow(y);</p><p><b>　　else</b></p><p>　　x=(handles.img);</p><p>　　y=fliplr(x);</p><

93、p>　　imshow(y);</p><p><b>　　end</b></p><p>　　------------------------------------------------------</p><p>　　function zeft90_Callback(hObject, eventdata, handles)</

94、p><p><b>　　%左轉(zhuǎn)90度</b></p><p>　　axes(handles.axes2);</p><p>　　x=(handles.img);</p><p>　　y=imrotate(x,90); </p><p>　　imshow(y);</p>

95、<p>　　-------------------------------------------------------</p><p>　　function right90_Callback(hObject, eventdata, handles)</p><p><b>　　%右轉(zhuǎn)90度</b></p><p>　　axes(h

96、andles.axes2);</p><p>　　x=(handles.img);</p><p>　　y=imrotate(x,-90); </p><p>　　imshow(y);</p><p>　　------------------------------------------------------</p

97、><p>　　function other_Callback(hObject, eventdata, handles)</p><p><b>　　%任意角度旋轉(zhuǎn)</b></p><p>　　axes(handles.axes2);</p><p>　　prompt={'輸入?yún)?shù)1:'};</p>

98、<p>　　defans={'30'};</p><p>　　p=inputdlg(prompt,'輸入?yún)?shù)',1,defans);</p><p>　　p1=str2num(p{1});</p><p>　　y=imrotate(handles.img,p1); </p><p&g

99、t;　　imshow(y);</p><p>　　-----------------------------------------------------</p><p>　　function gs_Callback(hObject, eventdata, handles)</p><p><b>　　%加入高斯噪聲</b></p>

100、;<p>　　axes(handles.axes2);</p><p>　　prompt={'輸入?yún)?shù)1:','輸入?yún)?shù)2'};</p><p>　　defans={'0','0.02'};</p><p>　　p=inputdlg(prompt,'輸入?yún)?shù)',1,defa

101、ns);</p><p>　　p1=str2num(p{1});</p><p>　　p2=str2num(p{2}); y=imnoise(handles.img,'gaussian',p1,p2);</p><p>　　imshow(y);</p><p>　　handles.noise_img=y;<

102、/p><p>　　guidata(hObject,handles);</p><p>　　-------------------------------------------------------</p><p>　　function jy_Callback(hObject, eventdata, handles)</p><p><b

103、>　　%加入椒鹽噪聲</b></p><p>　　prompt={'輸入?yún)?shù)1:'}; %對話框的設置，用戶輸入的是字符串</p><p>　　defans={'0.02'}; %缺省值</p><p>　　p=inputdlg(prompt,'輸入?yún)?shù)',1,defans); &

104、lt;/p><p>　　p1=str2num(p{1}); %字符串轉(zhuǎn)化為數(shù)值</p><p>　　axes(handles.axes2);</p><p>　　x=(handles.img);</p><p>　　y=imnoise(x,'salt & pepper',p1);

105、 </p><p>　　imshow(y);</p><p>　　handles.noise_img=y;</p><p>　　guidata(hObject,handles);</p><p>　　--------------------------------------------------------</p><

106、p>　　function cx_Callback(hObject, eventdata, handles)</p><p><b>　　%加入乘性噪聲</b></p><p>　　axes(handles.axes2);</p><p>　　prompt={'輸入?yún)?shù)1:'};</p><p>　　

107、defans={'0.02'};</p><p>　　p=inputdlg(prompt,'輸入?yún)?shù)',1,defans);</p><p>　　p1=str2num(p{1});</p><p>　　y=imnoise(handles.img,'speckle',p1);</p><p>　　

108、imshow(y);</p><p>　　handles.noise_img=y;</p><p>　　guidata(hObject,handles);</p><p>　　-----------------------------------------------------</p><p>　　function zhifangtuto

109、ngji_Callback(hObject, eventdata, handles)</p><p>　　---------------------------------------------------------</p><p>　　function red_Callback(hObject, eventdata, handles)</p><p><

110、b>　　%R直方圖</b></p><p>　　set(handles.axes2,'HandleVisibility','ON');</p><p>　　axes(handles.axes2);</p><p>　　x=imhist(handles.img(:,:,1));

111、 %直方圖統(tǒng)計</p><p>　　x1=x(1:10:256);</p><p>　　horz=1:10:256;</p><p>　　bar(horz,x1);</p><p>　　set(handles.axes2,'xtick',0:50:255);</p><p>　　-----------

112、----------------------------------------------</p><p>　　function gray_Callback(hObject, eventdata, handles)</p><p><b>　　%G直方圖</b></p><p>　　set(handles.axes2,'HandleV

113、isibility','ON');</p><p>　　axes(handles.axes2);</p><p>　　if isrgb(handles.img) </p><p>　　x=imhist(handles.img(:,:,2)); %直方圖統(tǒng)計</p><p>　　x1=x(1:10:256);<

114、/p><p>　　horz=1:10:256;</p><p>　　bar(horz,x1);</p><p>　　set(handles.axes2,'xtick',0:50:255);</p><p><b>　　else</b></p><p>　　msgbox('這是灰

115、度圖像','旋轉(zhuǎn)失敗');</p><p><b>　　end</b></p><p>　　------------------------------------------------------</p><p>　　function blue_Callback(hObject, eventdata, handles

116、)</p><p><b>　　%B直方圖</b></p><p>　　set(handles.axes2,'HandleVisibility','ON');</p><p>　　axes(handles.axes2);</p><p>　　if isrgb(handles.img) &l

117、t;/p><p>　　x=imhist(handles.img(:,:,3)); %直方圖統(tǒng)計</p><p>　　x1=x(1:10:256);</p><p>　　horz=1:10:256;</p><p>　　bar(horz,x1);</p><p>　　%axis([

118、0 255 0 150000]);</p><p>　　set(handles.axes2,'xtick',0:50:255);</p><p>　　%set(handles.axes2,'ytick',0:2000:15000);</p><p><b>　　else</b></p><p&

119、gt;　　msgbox('這是灰度圖像','旋轉(zhuǎn)失敗');</p><p><b>　　end</b></p><p>　　---------------------------------------------------------</p><p>　　function junheng_Callback(h

120、Object, eventdata, handles)</p><p><b>　　%直方圖均衡</b></p><p>　　set(handles.axes2,'HandleVisibility','ON');</p><p>　　axes(handles.axes2);</p><p>

121、;　　if isrgb(handles.img) </p><p>　　a=histeq(handles.img(:,:,1));</p><p>　　b=histeq(handles.img(:,:,2));</p><p>　　c=histeq(handles.img(:,:,3));</p><p>　　k(:,:,1)=a;</

122、p><p>　　k(:,:,2)=b;</p><p>　　k(:,:,3)=c;</p><p>　　imshow(k);</p><p><b>　　else</b></p><p>　　h=histeq(handles.img); %直方圖均衡</

123、p><p>　　imshow(h);</p><p><b>　　end</b></p><p>　　--------------------------------------------------------</p><p><b>　　%頻譜分析</b></p><p>

124、　　-----------------------------------------------------</p><p>　　function pinpu_Callback(hObject, eventdata, handles)</p><p>　　--------------------------------------------------------</p>

125、<p>　　function pinputu_Callback(hObject, eventdata, handles)</p><p><b>　　%顯示頻譜圖</b></p><p>　　axes(handles.axes2);</p><p>　　x=(handles.img);</p><p>　　i

126、f isrgb(x)</p><p>　　m=fft2(x(:,:,1));</p><p>　　y=fftshift(m);</p><p>　　imshow(log(abs(y)),[]);</p><p><b>　　else</b></p><p>　　m=fft2(x);</p&g