Nummer Eins

%%% Part 3
pwd
cd s:\DSPLab\Lab1

clear all

a=1; b=2; c=3; d=4;
save variables.mat
clear all
whos
load variables.mat

%%% Part 4
clear all

a=3-4j, b=real(a), c=imag(a), d=abs(a), e=angle(a)

f=4; g=9; h=sqrt(f)+j*sqrt(g)

z=5+7*j;
m=real(z), n=imag(z), p=abs(z), q=angle(z)

%%% Part 5
clear all

a=[3 4; 2 1]
b=[1.5 -2.4 3.5 0.7; -6.2 3.1 -5.5 4.1; 1.1 2.2 -0.1 0]

b(1)
e=b(2,3), f=b([2 3],[1,3]), g=b(2,[3 4])

h=[1 2 3], k=[4; 7], m=[5 6; 8 9]

n=[h; k m]

clear all

a=[3 5 9], b=[3; 5; 9]
c=2:5, d=3:2:9

x=0.5:0.25:2.0;
y=sqrt(x);
x,y

f=[10 5 4 7 9 0]
g=[2 5 6]; h=f(g)
m=[1.5 -2.4  3.5 0.7; -6.2  3.1 -5.5  4.1; 1.1 2.2 -0.1  0]
n=m(1:2,2:4), o=m(:, 1:2), p=m(2, :)

clear all
roots ([2 4 10])
clear all

x= (-4 + sqrt(4^2-4*2*10))/(2*2)
x1= (-4 - sqrt(4^2-4*2*10))/(2*2)

%%% Part 6
clear all

'Signal and System Analysis'

M='MATLAB Character String'
C=M(8:16)

c='We learn to use MATLAB in EE 266 Laboratory'
A=strcat(c(17:22),c(33:43))

%%% Part 7
clear all

a=[1 2; 3 4]; b=[3 1; 7 8]; c=[2 4];
d=a+b, e=c*a, f=a^2, g=c'
h=a\b,  k=b/a
m=a.*b, n=b./a, o=b.^a

clear all
a=[1.5 3.3; 6 -4.5; -2.5 .7]
b=[.5 .3; -.1 .2; .4 -.3]
c=[1 2; 1 2]
d=[3.1 1.4 -.3; -.5 1.6 .1]
a-b*c^2+2*d'

clear all
b=[1 2 3; 4 5 6]
c=[3 2 1; 6 5 4]
d=[9 8 7; 1 2 3]
for(ii=1:1:2)
   for(jj=1:1:3)
      b(ii,jj)-c(ii,jj)*d(ii,jj)^4
   end
end

2-2*8^4

%%% Part 9
clear all
a=[1 3 2; 4 6 5], b=a>2&a<=5

c=[1 5 3 4 7 8], d=c>4

clear all
a=[1.2 -3.2 24; .6 -.3 -.5; -2.3 1.6 20]
c=a<2&a>=-1&a~=.6
b=a.*c
b(2,2)=2*b(2,2)*(b(2,2)<0);
b(2,3)=2*b(2,3)*(b(2,3)<0);
b

%%% Part 10
clear all
t=0.1;
x=2^t*sqrt(t)-sin(2*t)/3
y=2^(t*sqrt(t))-sin(2*t)/3
f=0:2:4; w=2*pi*f;
X=(3-j*0.1*w)./(1.5+j*0.2*w)

t= 0:0.5:2;
x=(t+1).*(t>=0&t<1)+2*(t>=1&t<=2)

clear all

x=[-3.6, -2.5, -1.4, -1, 0, 1.4, 2.5, 3.6]
round(x)
floor(x)
ceil(x)

clear all

x=(log(2+(sin(3))^2)+exp(-.2))/(sqrt(2^(1.6)+3^(-.5)))

t=-1.2
while(t~=1.6)
   w=3*t^3+2*t^2-t+sin(t)
   x=0*(t<0)+2*(t>=0)
   t=t+.4;
end

%%% Part 11
clear all

t=0:0.2:0.8; x=zeros(size(t));
for k=1:3;
    x=x+sqrt(k)*t.^sqrt(1.2*k);
end;

for m=1:3;
    for n=1:4;
        y(m,n)=m+n;
    end;
end;

clear all
t=0:.5:10; y=zeros(size(t));
for q=1:length(t);
    y=t+q;
end;

clear all
n=0;
d=5;
if n==2;
    y=10*d;
else;
    y=0;
end;

clear all

n=1;
while 2*n<5000; n=2*n;end;

for k=1:4;
    if k==1; x(k)=3*k;
        else if k==2|k==4;
            x(k)=k/2;
        else;
            x(k)=2*k;
        end;
    end;
end;

clear all

c='t'; n=2;
if c=='f'; c='false'; y=Nan; end;
d=0.1:0.1:0.4;
if c=='t';
    if n==2;
       y=10*d(n);
    else;
       y=0;
    end;
end;
if n==2;
    y=10*d(n);
else;
    y=0;
end;

%number2
clear all
t=0;
f=10;
x=[0 0 0; 0 0 0; 0 0 0; 0 0 0; 0 0 0];
for k=1:1:5
   for m=1:1:3
       x(k,m)=3*cosd(2*pi*f*t+.1)
       f=f+5;
   end;
   f=10;
   t=t+.1;
end;

%3
clear all
t=0;
for w=35:5:45
    while(exp(1.2)*cos(w*t)<10 && t^3 <10)
        t=t+.01;
    end;
    t=t-.01
    t=0;
end;

%4
clear all
t=-1;
x=zeros(1,10001);
for k=1:1:10001
    x(k)=exp(-abs(t));
    t=t+.0002;
end;

%%% Part 12
clear all

a=[1 0 2 3 0 4];
b=find(a)
n=find(a>2)
size(a)
c=zeros(size(a))
max(a)
min(a)
mean(a)
sum(a)

d=-0.1:0.1:0.2;
dm=meshgrid(d,1:3)

%2
clear all
t=0;
for k=1:1:11
    x(k)=4*cos(2*pi*t+.2)+3*sin(pi^2*t);
    t=t+.1;
end;
max(x)
min(x)
mean(x)
find(x>4)

%3
clear all
A=[1 4 3 2; 4 1 2 5; 3 3 5 1];
size(A)
max(max(A))
min(min(A))
max(A,[],2)
min(A,[],2)
mean(A)
mean(mean(A))

%%% Part 13
clear all
t=0:0.01:10;
f1=0.2;
f2=0.425;
[s1, s2, s3]=sumsin(t,f1,f2);
plot(t,s1);
hold on
plot(t,s2,'-r');
plot(t,s3,'-g');
xlabel('t');
ylabel('y');
title('AJP Lab 1 Sinusoids');
legend('s1','s2','s3')
hold off
figure
x1=t; x2=t; x3=t;
subplot(3,1,1); plot(x1,s1);
title('AJP Lab 1 Sinusoids');
ylabel('y1');
subplot(3,1,2);plot(x2,s2);
ylabel('y2');
subplot(3,1,3);plot(x3,s3);
xlabel('t');
ylabel('y3');