MATLAB에 의한 ANALOG 신호처리

소개글

MATLAB에 의한 ANALOG 신호처리에 대한 보고서 자료입니다.

목차

1. analog 신호의 표시 방법

2 중요한 신호 생성용 함수

본문내용

ct(t/2)+delta(t+1);
rt=cumsum(st)*0.001;
figure(1);
plot(t,rt);
% 미분에 의한 확인
[m,n]=size(st);
st2(2:n)=diff(st)/0.001;
st2(1)=0;
figure(2);
plot(t,st);
hold on; plot(t,st2,‘-r’);hold off;
위 MATLAB 프로그램을 실행하면 2 개의 그림이 생성되는 데, 첫 그림은 적분한 신호 파형이고 두 번째 그림은 원 신호 와 적분한 신호를 다시 미분 한 것을 겹쳐 그린 것이다. 2 번째 그림의 화면이 모두 붉은 색이면 2 신호가 일치 하는 것을 알 수 있다.
위 프로그램에서 st를 다양하게 바꾸어서 확인 해 볼 수 있다. 즉
(a) st=exp(-abs(t));
(b) st=cos(2*pi*t);
(c) st=sinc(t);
(d) st=exp(-2*t) .* unit(t);
(e) st=exp(j*2*pi*t)+exp(-j*2*pi*t);
의 경우에 대해 확인 해 보라. sinc 함수는 내장된 함수를 사용하면 되고, (e)의 경우
st 가 복소수 이지만 실수부만 취하여 그래프를 그린다.
(나) 다음 적분값을 구하라
(a)
(b)
(c)
(d)
(sol)
(a)
sigma=0.5;
x=-3:0.001:7;
s=exp(-(x-2) .^2/(2*sigma^2));
s=1/(sqrt(2*pi)*sigma)*s;
sum(s)*0.001
이론적 계산결과는 1이다. 위 프로그램 실행결과는 값이 0.5 와 1일때는 1이 되지만 값이 2 일때는 1.9877 이 된다. 이 경우는 x의 초기치와 최종치의 범위를 증가 시키면 2에접근한다.
(b) (a)의 셋째 줄을 다음과 같이 수정
s=x .* exp(-(x-2) .^2/(2*sigma^2));
이론적 계산 결과는 값에 관계없이 2 이다. 실행 결과는 (a) 의 경우와 비슷하다.
(c) (a)의 3 째줄을 다음과 같이 수정.
s=(x-2) .^2 .* exp(-(x-2) .^2/(2*sigma^2));
이론적 계산결과는 이 된다 . 실행 결과는 (a)의 경우와 비슷하다.
(d)
f2=1/(2*pi);
f=0:0.0002:f2;
s=1 ./(4+4*pi^2 *f .^2);
s1=sum(s)*0.0002
s2=1/(4*pi) * atan(0.5)
이론적 계산결과는 =0.036896 이고 실행결과는 0.03691로 거의 같은 값이 된다.

키워드

MATLAB,   에 의한 ,   ANALOG 신호처리