我要用一下这些代码写实验报告

%%%%实验一：%%%%
x=1:1:100;
y=5*ones(1,100);
y1=awgn(y,1);        %加入高斯白噪声
y2=smooth(y1,21);    %移动平均滤波
%画图
hold on
plot(y,'k','linewidth',1);
plot(x,y1,'r','linewidth',1);
plot(x,y2,'g','linewidth',1);
title('直流信号');
xlabel('时间');
ylabel('电流');

　　

%%%%实验二：%%%%
x=0:0.1:1;
y=[1,2,3.2,6.2,7,7.6,8,9.2,9.3,9.4,15];
%一次拟合
h1=polyfit(x,y,1);
z1=polyval(h1,x);
%二次拟合
h2=polyfit(x,y,2);
z2=polyval(h2,x);
%三次拟合
h3=polyfit(x,y,3);
z3=polyval(h3,x);
%画图部分
plot(x,y,'k*',x,z1,'r',x,z2,'g',x,z3,'b');
%求平均拟合误差
sum1=0;
sum2=0;
sum3=0;
for i=1:11
    s1=z1(i)-y(i);
    sum1=sum1+abs(s1);
    s2=z2(i)-y(i);
    sum2=sum2+abs(s2);
    s3=z3(i)-y(i);
    sum3=sum3+abs(s3);
end
error1=sum1/11;
error2=sum2/11;
error3=sum3/11;

　　

%%%% 实验三：%%%%
x=0:0.01:15/16;
y=cos(4.8*pi*x);          %原图

fs=16;                    %采样频率
N=16;                     %采样点数
n=-N/2:N/2;
t=n/fs;                   %采样时间
g=cos(4.8*pi*t);          %采样函数
G=fftshift(fft(g,N));     %傅里叶变换
m=abs(G(1:N))*2*pi/N;     %求真实幅值
f=n*fs/N;                 %频域间隔

%画时域图
figure(1);
subplot(1,2,1);
plot(x,y,'k','linewidth',1);            %原图
title('原图');
xlabel('时间/s');
ylabel('幅值');
subplot(1,2,2);
plot(t,g,'b','linewidth',1);             %采样后图
title('采样后图');
xlabel('时间/s');
ylabel('幅值');

%画频谱图
figure(2);
subplot(2,1,1);
stem(f(1:N),m(1:N));
title('DFT图');
xlabel('频率/Hz');
ylabel('幅值');
subplot(2,1,2);
plot(f(1:N),abs(G),'linewidth',1);
title('FFT图');
xlabel('频率/Hz');
ylabel('幅值');

　　

嘿嘿