untitled.m

实现正弦波生成的matlab程序

clear all;
clc
N=50;
w=pi/6
b=cos(w);
A=[b b+1;b-1 b];
S=zeros(N,2);
S(1,1:2)=[0.4 0.3];
for i=2:N
     S(i,:)=A*S(i-1,:)';
end
subplot(2,1,1)
stem(S(:,1));
xlabel('时间序号')
ylabel('振幅')
title('s1[n]')
subplot(2,1,2)
stem(S(:,2))
xlabel('时间序号')
ylabel('振幅')
title('s2[n]')