QPSK modulation

%try cosine wave plot also
%a=.5;  
pi=3.14;
f=1;

x=[0 0 0 1 1 0 1 1 0 0 1 0  1 1 0 1];
nx=size(x,2);

i=1;
while i<nx+1
    t = i:0.001:i+2;
    if x(i)==0 && x(i+1)==0
        z=sin(2*pi*f*t-3*pi/4);
        dibit=0;
    elseif x(i)==0 && x(i+1)==1
        z=sin(2*pi*f*t-pi/4);
        dibit=1;
    elseif x(i)==1 && x(i+1)==0
        z=sin(2*pi*f*t+pi/4);
        dibit=2;
    else
        z=sin(2*pi*f*t+3*pi/4);
        dibit=3;
    end
    %subplot(2,1,1);
    plot(t,dibit,'r');
    %grid on;
    hold on;
    %subplot(2,1,2);
    plot(t,z);
    hold on;
    grid on;
    axis([1 20 -3 3]);
    i=i+2;

end

FM modulation

Ac=1;
fc=1;

Am=1;
fm=0.1;
kf=1.5;
B=kf*Am/fm;

t=[0:0.1:20];
ct=Ac*cos(2*pi*fc*t);
mt=Am*cos(2*pi*fm*t);
FMt=Ac*cos(2*pi*fc*t+B*sin(2*pi*fm*t));

subplot(2,2,1);
plot(mt);
ylabel('message signal');
grid on;

subplot(2,2,2);
plot(ct);
ylabel('carrier');
grid on;

subplot(2,2,3);
plot(FMt);
ylabel('FM signal');
grid on;

AM modulation

clc;
Ac=2;
fc=.9;

Am=.5;
fm=.05;
Fs=100;

ka=1;

t=[0:0.1:50];
ct=Ac*cos(2*pi*fc*t);
mt=Am*cos(2*pi*fm*t);
DSBAMt=ct.*(1+ka*mt);

subplot(2,2,1);
plot(mt);
ylabel('message signal');
grid on;

subplot(2,2,2);
plot(ct);
ylabel('carrier');
grid on;

subplot(2,2,3);
plot(DSBAMt);
ylabel('DSBAM signal');
grid on;