function x=modulation2ask(a,fc,r) n=2000; i=linspace(0,1,n); x1=sin(2*pi*i*fc/r); an=length(a); x0=zeros(size(an*n)); for j=1:an if a(j) y=x1; else y=x0; end x(1,(j-1)*n+1:j*n)=y; end nn=(1:length(x))/n; plot(nn,x);grid on;xlabel('2ASK已调制信号'); title('20211260562胡盈盈'); deat=x; for j=1:an deat(1,(j-1)*n+1:j*n)=x(1,(j-1)*n+1:j*n).*x1; end subplot(233); plot(nn,deat);grid;xlabel('2ASK解调 1-与载波相乘'); 接着写2ASK解调 2-低通滤波；2ASK解调 恢复的信号代码

2ASK解调 2-低通滤波的代码如下：

fcut = fc*2; % 设置低通滤波器截止频率为载波频率的2倍
[b,a] = butter(6, fcut/(r*n/2), 'low'); % 计算6阶巴特沃斯低通滤波器系数
y = filter(b,a,deat); % 对解调信号进行滤波

2ASK解调 恢复的信号代码如下：

x_rec = zeros(1,an); % 初始化恢复的数据
for j = 1:an
x_rec(j) = mean(y((j-1)*n+1:j*n))/mean(x1); % 计算每个符号的平均值，并除以载波幅值
end

其中，x_rec即为解调后恢复的数据。