用MATLAB语言，分两种情况f s≥2f max和f s<2f max描述采样定理

在MATLAB语言中，可以使用以下代码来描述采样定理的两种情况：

当采样频率 f_s 大于等于信号最高频率 2f_max 时，即 f_s ≥ 2f_max，可以使用以下代码进行采样：

% 生成信号
t = 0:0.001:1; % 定义时间范围为 0 到 1 秒，采样率为 1000 Hz
fmax = 50; % 定义信号的最高频率为 50 Hz
x = sin(2*pi*fmax*t); % 生成正弦信号

% 对信号进行采样
fs = 200; % 定义采样率为 200 Hz
Ts = 1/fs; % 计算采样周期
n = 0:Ts:1; % 定义采样时间范围
xn = sin(2*pi*fmax*n); % 采样信号

% 绘制信号和采样信号的图形
subplot(2,1,1);
plot(t,x);
title('原始信号');
xlabel('时间（秒）');
ylabel('幅度');
subplot(2,1,2);
stem(n,xn);
title('采样信号');
xlabel('时间（秒）');
ylabel('幅度');

当采样频率 f_s 小于信号最高频率 2f_max 时，即 f_s < 2f_max，会出现混叠现象。可以使用以下代码进行采样：

% 生成信号
t = 0:0.001:1; % 定义时间范围为 0 到 1 秒，采样率为 1000 Hz
fmax = 50; % 定义信号的最高频率为 50 Hz
x = sin(2*pi*fmax*t); % 生成正弦信号

% 对信号进行采样
fs = 80; % 定义采样率为 80 Hz
Ts = 1/fs; % 计算采样周期
n = 0:Ts:1; % 定义采样时间范围
xn = sin(2*pi*fmax*n); % 采样信号

% 绘制信号和采样信号的图形
subplot(2,1,1);
plot(t,x);
title('原始信号');
xlabel('时间（秒）');
ylabel('幅度');
subplot(2,1,2);
stem(n,xn);
title('采样信号');
xlabel('时间（秒）');
ylabel('幅度');

在这种情况下，采样信号中会出现频率大于 f_s/2 的混叠频率。可以使用数字信号处理技术进行滤波，去除混叠频率。