4

A=0.5;
f=3000;
N=4;
Fs1=100000;
Ts1=1/Fs1;
n=0:ceil(N*Fs1/f);
t=linspace(0,4/f,1000);
Fs=[4000,5000,6000,10000,12000];
msi=[0,0,0,0,0];
z=A*cos(2*pi*f*Ts1);
for i =1:length(Fs)
    Ts=1/Fs(i);
    n1=0:ceil(N*Fs(i)/f);    %==0:1:ceil(N*Fs/f)
    x_s= A*cos(2*pi*f*n1*Ts);
    b=0:max(n1);
    v=0:(max(n1)/max(n)):max(n1);
    x=0:max(n1);

    % rcons_signal=interp1(x,x_s,v,'linear');
    % subplot(4,1,1);
    % plot(n,rcons_signal);
    % xlabel("t");
    %ylabel("X");
    %title('Reconstructed signal linear');
    %mse=mean((rcons_signal-z).^2);


    rcons_signal=interp1(x,x_s,v,'cubic');
    subplot(5,1,i);
    plot(n,rcons_signal);
    xlabel("t");
    ylabel("X");
    title('Reconstructed signal cubic',num2str(Fs(i)));
    mse(i)=mean((rcons_signal-z).^2);

    %rcons_signal=interp1(x,x_s,v,'spline');
    % subplot(4,1,3);
    % plot(n,rcons_signal);
    % xlabel("t");
    % ylabel("X");
    % title('Reconstructed signal spline');
    %mse=mean((rcons_signal-z).^2);

    %rcons_signal=interp1(x,x_s,v,'pchip');
    %subplot(4,1,4);
    % plot(n,rcons_signal);
    % xlabel("t");
    %ylabel("X");
    % title('Reconstructed signal pchip');
    %mse=mean((rcons_signal-z).^2);
end