차례 제 1 장 MATLAB 연습 MATLAB에서사용되는기호들

Size: px

Start display at page:

Download "차례 제 1 장 MATLAB 연습 MATLAB에서사용되는기호들"

병만 유
5 years ago
Views:

2 차례 제 1 장 MATLAB 연습 MATLAB에서사용되는기호들 연산자 (operators) MATLAB에서기본적으로사용하는유용한함수및명령어들 데이터입, 출력방법 명령어흐름제어 그래프그리는방법 실습 제 2 장디지털신호처리기본개념 디지털신호처리란? 디지털신호처리사용분야 신호의분류 디지털시스템 샘플링정리 Aliasing 양자화 (quantization) 와 A/D 변환기 실습 제 3 장디지털신호표현 디지털기본신호표현 단위임펄스신호 (Unit Impulse signal) 단위계단신호 (Unit Step signal) 실수지수신호 (Real exponential signal) 복소지수신호 (Complex exponential signal) 랜덤신호 (Random signal) 임의신호발생기 White Gaussian 잡음발생기 이산신호연산을위한몇가지편리한함수들 실습 제 4 장선형시스템 이산시스템 (Discrete System) 콘볼루션 (Convolution) 과상관 (Correlation)

3 4.2.1 콘볼루션 상관 (Correlation) 실습 제 5 장임펄스응답과선형차분방정식 임펄스응답 (Impulse Response) 간단한몇가지임펄스응답예 선형차분방정식 차분방정식을통하여주파수응답을구하는방법 실습 제 6 장 z- 변환 (z-transform) z- 변환의정의 DTFT의일반화에의한정의 z-변환의정의 전달함수에의한정의 극점과영점 수렴영역 (Region of Convergence : ROC) 안정도판별 z-변환의성질 (property) 역 Z-변환 (Inverse Z-Transform) 실습 제 7 장이산푸리에변환 (Discrete Fourier Transform) 아날로그주기신호 아날로그비주기신호 이산시간푸리에변환 (DTFT) 이산푸리에변환 (DFT) 이산푸리에변환정의 푸리에분석의중요성 고밀도 (high-density) 스펙트럼과고분해능 (high-resolution) 스펙트럼 7.6 이산푸리에변환의응용 고속푸리에변환 (Fast Fourier Transform : FFT) DIT 방식

4 실습 제 8 장디지털필터 H/W 구조 디지털필터시스템구성을위한기본소자표현 FIR 필터구조 직접형 (Direct Form) 구조 전치형 (Transpose Form) 구조 직렬연결형 (Cascade Form) 구조 IIR 필터구조 직접형 I(Direct Form I) 구조 직접형 II(Direct Form II) 또는케노닉형 (Canonic Form) 구조 직렬연결형 (Cascade Form) 구조 병렬형 (Parallel Form) 구조 실습 제 9 장 FIR 필터설계 디지털필터설계단계 디지털필터의사양결정 FIR 필터나 IIR 필터의결정 FIR 필터설계 주파수샘플링법 창함수법 등리플필터를위한최적근사화법 실습 제 10 장 IIR 필터설계 아날로그필터설계 아나로그필터에서디지털필터로의변환방법 임펄스불변변환 (impulse invariance transformation) 쌍일차변환 (bilinear transformation) 법 주파수대역변환 실습 제 11 장적응필터 (Adaptive Filter)

5 실습 제 12 장디지털신호처리응용 과제 제 13 장디지털신호처리응용 과제 제 14 장디지털신호처리응용 과제 참고문헌

6 1 MATLAB Mathworks MATLAB. MATLAB MATLAB Toolbox. MATLAB MATLAB C C++., MATLAB C C++, DSP TMS320 Mathworks MATLAB.. 1. C C++ 2. TI TMS320C6000 Motolora M Xilinx ALTERA EPLD(FPGA) HDL 4. x-pc MABLAB Emulation Simulink Labview MATLAB. MATLAB. MATLAB, MATLAB. MATLAB 1. Steven W. Smith, "The Scientist and Engineer's Guide to Digital Signal Processing" 2. Vinay K. Ingle, John G. Proakis, "Digital Signal Processing using the MATLAB" 3.

7 1.1 MATLAB. 1-1 ' ' ' [ ] ( ) ; :. (matrix) (transpose). ) a', A'. ) title('magnitude'),. ) a=[ ], (function). ) a(1,2), max(a). ) a=sqrt(2)+exp(3);. ) t=[0:0.01:1];. ) c=a.^2, c=a./b % (comment line).... >> (variable) = (expression)

8 . 1-2 pi i j inf NaN (, ),,, i j. (infinity). ( nonzeros ) / 0 (Not-a-number). 0 / 0 inf / inf 1) >>A=[ ] % A = >>A=[1; 4; 8; 6;] A = % >>a=[1:4:1] % a = >>a=[9:-2:1] % a =

9 2) >>A=[1 1 1;1 2 3;1 3 6] >>A= A = % %,. 3) >> a=[1 2 3]; b=[4 5 6]; >> c= a.* b c = >> a=[ ]; b=[ ]; >> y=a.^2 + b.^2 - a + b y = >>A=[1 1 1;1 2 3;1 3 6]; >>sum(a) % sum(a). ans = ) >> x=[-1 sqrt(3) 4/5] x =

10 'x' '( )' (subscript). >>x(3) ans = >> a=[1 2 3; 4 5 6; 7 8 9] a = >> a(1,3) ans = 3 >> a(2,2) ans = 5 >> a(3,2) ans = 8 >> a(1,:) % 1 ans = >> a(:,2) % 2 ans = >> a(3,2:3) % 3 2,3 ans = 8 9

11 >> a(2:3,2) % 2 2,3 ans = 5 8,.... >> P=[1, 2, 4, 6, 8]+[pi, 4, exp(1), 0, -1]+... [cos(0.1*pi), sin(pi/3), tan(3), atan(2), sqrt(pi)] P = ) >> a=[ ; ; ] a = >> a(1,3)= % a = >> a(2,6)= % % 6 6. a = >> a(3,5)=-55 % a =

12 >> b=2*a % a 2 b = ) (transpose : i j j i ) x " z' " (complex conjugate). x " z.' ". ". (dot) ". >>z=[1+2i 3+4i] z = i i >>z' ans = i i >>z.' ans = i i

13 1.2 (operators) 1) (Arithmetic Operators) 1. transpose( ".' " ), power ( ".^ " ), transpose( " ' " ), power( " ^ " ) 2. unary plus( " + " ), unary minus( " - " ) 3. ( " * " ), ( ". / " ), ( ". \ " ), ( " * " ), ( " / " ), ( " \ " ) 4. ( " + " ), ( " - " ) 5. colon ( " : " ) >>A=[3 9 5]; >>B=[2 1 5]; >>C=A./B.^2 % B.^2, A./. C = ) (Relational operators) (less), (not equal to)" < " (Less than)" <= " (Less than or equal to)" > " (Greater than)" >= " (Greater than or equal to)" == "(Equal to)" ~= " (Not equal to)"

14 >>A=[2 7 6;9 0-1;3.5 6]; % MATLAB 0. %, >>B=[ ;-3 2 5;4-1 7]; >>C= A<B C = " find ".,.,,, " find ",. >>A=magic(4) % magic. A = >>index=find(a > 8); >>A(index)=100 A =

15 3) (Logical Operators) "AND", "OR" TRUE", "FALSE" & AND OR ~ NOT 1.3 MATLAB sin(x) / cos(x) / tan(x) / asin(x) / acos(x) / atan(x) / exp(x) / log(x) / log10(x) / sqrt(x) / abs(x) / phase(x) / real(x) / imag(x) / conj(x) / mean(x) / std(x) / sum(x) / max(x) / min(x) / ceil(x) / floor(x) / round(x) / rem(x,y) / fix(x) / length(x) / size(x) / ones(rn,cn) / zeros(rn,cn) / eye(rn,cn) / linspace(sn,incn,en) / pi(= ) x=3+3j or x=3+3i or x=3+3*j or x=3+3*i clear all; close all; clc; disp('dsp').

16 1.4, 1) load save MATLAB, save.m, load.m. save.m *.mat file, load.m *.mat file. >> x=eye(3); y='dsp'; z=2+2i; >> save Data_file x y z, x, y, z DataFile.mat file. >>clear all >>load Data_file % base workspace., Data_file.mat file x, y, z base workspace. >> save output_filename variables (xxx.mat ) >> load data_filename.mat ASCII, >> save output_filename.dat variables -ascii( /ascii : ascii ) >> load data_filename.txt. 2) Low_level file I/O fopen.m fclose.m.

17 fprintf.m fscanf.m,. 3) input input. >> a=input(' a= ') a= 6 a = 6 >> x=a+a.^2; >> x x = ) FOR for i = 1:10 end ; x=0; for i = 1:10 x=x+i; end x x = 55 2) IF if ( ) ; elseif ;

18 else end ; x=0; y=0; z=0; for k=0:1:100; if (2*round(k/2)==k) x=x+1; elseif (3*round(k/3)==k) y=y+1; else z=z+1; end end x y z x = y = z = ) WHILE while ( ). ; end % n=1; while (prod(1:n)<1000) n=n+1; end n n = 7

19 4) SWITCH switch case 1, ; case 2, ;... otherwise ; end % 1 % 2 % case x=input(' x=') y=rem(x,2); switch y case 1, disp(' ') case 2, disp(' ') end x=45 x = MATLAB

20 MATLAB. "plot" stem". PLOT 1) >> t=[0:0.1:2*pi]; >> y=sin(t); >> plot(t,y)

22 >> t=[0:0.1:2*pi]; >> y=sin(t); >> z=cos(t); >> plot(t,y,t,z) 3) hold >> hold on, hold on >> hold off.

23 >> x=0:.1:2*pi; >> y=sin(x); >> plot(x,y) >> grid on >> hold on >> plot(x,exp(-2*x),'r:*') STEM >> n=[0:0.2:2*pi]; >> y=sin(n); >> stem(n,y) >> stem(n,y,'filled')

24 log x, y log,. >>w=logspace(-2,2,1000); >>H=1./(1+w.^2).^0.5; >>subplot(1,3,1),semilogx(w,h); >>subplot(1,3,2),semilogy(w,h); >>subplot(1,3,3),loglog(w,h); subplot(,, ). >> subplot(2,2,1), plot(1:10) >> x = 0:.1:2*pi; >> subplot(2,2,2), plot(x,sin(x)) >> subplot(2,2,3), plot(x,exp(-x),'r') >> subplot(2,2,4), plot(x+sin(x)+exp(-x),'k')

25 'option'. >> plot(x,y,'option') >> x = 0:.1:2*pi; >> plot(x,sin(x),'r--d')

26 . 1-5 option option option option b blue. point - solid g green o circle : dotted r red x x-mark -. dashdot c cyan + plus -- dashed m magenta * star y yellow s square k black d diamond v ^ triangle (down) triangle(up) < triangle(left) > p triangle(right) pentagram h hexagram

27 >> x = 0:.1:2*pi; >> plot(x,sin(x),'r--d'); >> xlabel('frequency(w)'); >> ylabel('amplitude[v]'); >> title('sin function'); 1) >> xlabel('^ {super } and _ {sub }') 2) >> ylabel('some \bfboldface\rm and some \it {italics }') 3) >> text(3.5, 20, 'y = e^ { \alphax}sin( \betax)')

28 . 1-6.

29 Property Editor 1) Figure Property 2) Axes Property

31 1 MATLAB. 1. MATLAB 2. MATLAB MATLAB MATLAB MATLAB, 3. (1) MATLAB meditor. (2) MATLAB. (3).

32 [ 1].,. [ 2] sin cos. y1 y2.

33 [ 3] sin cos save ASCII. load,. "o". [ 4] "input" a,b,c,d,e a) b) log c) d) e) MATLAB..

<4D6963726F736F667420506F776572506F696E74202D204D41544C4142B0ADC0C7B7CF28B9E8C6F7BFEB295F3031C0E55FBDC3C0DBC7CFB1E22E707074205BC8A3C8AF20B8F0B5E55D>

<4D6963726F736F667420506F776572506F696E74202D204D41544C4142B0ADC0C7B7CF28B9E8C6F7BFEB295F3031C0E55FBDC3C0DBC7CFB1E22E707074205BC8A3C8AF20B8F0B5E55D> MATLAB MATLAB 개요와 응용 1장 MATLAB 시작하기 10 5 0 황철호 -5-10 30 20 10 0 0 5 10 15 20 25 MATLAB 시작하기 이장의내용 MATLAB의여러창(window)들의 특성과 목적 기술 스칼라의 산술연산 및 기본 수학함수의 사용. 스칼라 변수들(할당 연산자)의 정의 및 변수들의 사용 방법 스크립트(script) 파일에