1. 设)(n x 和)(n y 是有限长的序列，]1.0,1,1.0,1[)(-=↑n x ，]1,1.0,1,1.0[)(-=↑
n y ，箭头所指位置表示n ＝0的序列值，箭头右边依次是n ＝1、2、3 ┉，箭头左边依次是n ＝－
1、－
2、－3 ┉。

求这两个序列的线性相关函数。

（分布使用直接计算法和表格法求解）
2. 试采用傅里叶变换对一段医学信号进行处理。

要求：有原信号波形、源代码和处理结果。

肌电原始信号：
做傅里叶变换：
N=10000;
M=1;
y1=fft(x1,N);
subplot(4,1,1) ;
plot(f(1:N/2),y1(1:N/2));
axis([0 500 0 20]);
grid on;
3.试采用频谱分析对一段医学信号进行处理。

要求：有原信号波形、源代码和处理结果。

肌电原始信号：
做频谱分析：
clear;
close all;
%fft 频率分析
a=load('EMG.txt');
y=fft(a,10000); %做10000点福利叶变换fs=1000;
N=length(y);
mag=abs(y);
f=(0:N-1)/N*fs;
figure;
plot(f,mag);%做幅频谱
xlabel('频率');
ylabel('幅值');
title('肌电幅频 N=10000');
figure;
plot(f,angle(y));% 做相频谱
xlabel('频率');
ylabel('相位');
title('肌电相频 N=10000');
grid on;
频率幅值
频率相位肌电相频 N=10000
做功率谱分析：
clear;
close all;
a=load('EMG.txt');
y=fft(a,10000); %做10000点傅里叶叶变换 fs=1000;
N=length(y);
mag=abs(y);
f=(0:N-1)/N*fs;
power1=(mag.^2)/10000;%周期图法求功率谱 figure;
plot(f,power1);
xlabel('频谱');
ylabel('功率谱');
title('肌电信号功率谱');
grid on;
频率功率谱肌电信号功率谱。