信号处理仿真(MATLAB)实验1-6

信号处理仿真(MATLAB)

实验指导书

青岛大学自动化工程学院电子工程系

2006年4月

MATLAB 实验一

一、实验目的：

1. Be familiar with MATLAB Environment

2. Be familiar with array and matrix

3. Be familiar with MATLAB operations and simple plot function

二、实验内容：

1. Be familiar with Matlab 6.5

Startup Matlab 6.5, browse the major tools of the Matlab desktop ? The Command Windows

? The Command History Windows ? Launch Pad

? The Edit/Debug Window ? Figure Windows

? Workspace Browser and Array Editer ? Help Browser

? Current Directory Browser

PART I:下列选择练习，不需提交实验报告

1． Give the answer of the following questions for the array

1.10.0

2.1

3.56

0.01.16.62.83.412.10.10.30.41.31.45.10.01.10.0

a r r a y -????-?

?=??-??-?? 1）What is the size of array1? size(array) 4 5

2）What is the value of array1(4,1)?

3）What is the size and value of array1(:,1:2)? 1.1000 0

0 1.1000 2.1000 0.1000 -1.4000 5.1000

4） What is the size and value of array1([1 3], end)? 6.0000 1.3000

2． Give the answer of the following commad 1) a=1:2:5;

a = 1 3 5

2) b=[a ’ a ’ a ’]; b =

1 1 1 3 3 3 5 5 5 3) c=b(1:2:3,1:2:3); c = 1 1 5 5 4) d=a+b(2,:)

d = 4 6 8

5) w=[zeros(1,3) ones(3,1)’ 3:5’]

w = 0 0 0 1 1 1 3 4 5 3． Give the answer of the sub-arrays

1.10.0

2.1

3.56

0.01.16.62.83.412.10.10.30.41.31.45.10.01.10.0

a r r a y -????-?

?=??-??-?? 1) array1(3,:); 2.1000 0.1000 0.3000 -0.4000 1.3000 2) array1(:,3); 2.1000 -6.6000 0.3000 0

3) array1(1:2:3,[3 3 4])

2.1000 2.1000 -

3.5000 0.3000 0.3000 -0.4000 4) array1([1 1],:)

1.1000 0

2.1000 -

3.5000 6.0000 1.1000 0 2.1000 -3.5000 6.0000 4． Give the answer of the following operations

2

21

11

,,,(2

)12022a b c d e y e --??????====?

?????--??????

1) a+b

3 -3 -1

4 2) a*d 2 -2 -1 2

3) a.*d 2 0 0 2 4) a*c 6 -5 5) a.*c wrong 6)a\b

1.0000 1.0000 0.5000 1.5000

7) a.\b

0.5000 0.5000 0 1.0000 8)a.^b

2.0000 -0.5000 1.0000 4.0000

PART I:下列选择练习，不需提交实验报告

1．Edit & Run the m-file

% test step response function

wn=6; kosi=[0.1:0.1:1.0 2]; figure(1); hold on for kos=kosi

num=wn^2; den=[1,2*kos*wn,wn.^2]; step(num,den);step 阶跃响应num 分子den 分母适量矩阵

end

hold off;

012345678910

0.20.40.60.811.2

1.41.6

1.8Step Response

Time (sec)

A m p l i t u d e

2．Edit & Run the m-file

% test plot function

x=0:pi/20:3*pi; y1=sin(x); y2=2*cos(2*x); plot(x,y1,'rv:',x,y2,'bo--');

title('Plot the Line of y=sin(2x) and its derivative'); xlabel('X axis'); ylabel('Y axis'); legend('f(x)','d/dx f(x)'); grid on;

Plot the Line of y=sin(2x) and its derivative

X axis

Y a x i s

3． Edit & Run the m-file

% test subplot and loglog function

x=0:0.1:10; y=x.^2-10.*x+26; subplot(2,2,1); plot(x,y); grid on;

subplot(2,2,2); semilogx(x,y); grid on; subplot(2,2,3); semilogy(x,y); grid on; subplot(2,2,4); loglog(x,y); grid on;

0510

1020

3010-110

10

1

010

20

30

0510

10

10

1

10

2

10

-110

10

1

10

10

1

10

2

4．Edit & Run the m-file

% test max and plot function

volts=120; rs=50; rl=1:0.1:100;

amps=volts./(rs+rl); pl=(amps.^2).*rl; [maxvol,index]=max(pl); plot(rl,pl,rl(index),pl(index),'rh'); grid on;

0102030405060708090100

10203040506070

MATLAB 实验二

一、实验目的：

1. Learn to design branch statements program

2. Be familiar with relational and logical operators

3. Practice 2D plotting

二、实验内容：

PART I: (选择练习，不需提交实验报告)

1. Hold command exercise x=-pi:pi/20:pi;

y1=sin(x); y2=cos(x); plot(x,y1, 'b-'); hold on; plot(x,y2, 'k--'); hold off; legend ('sinx', 'cosx')

-4

-3

-2

-1

1

2

3

4

-1-0.8-0.6-0.4-0.200.20.40.60.8

2. Figure command exercise

figure(1);

subplot(2,1,1);

x=-pi:pi/20:pi; y=sin(x); plot(x,y); grid on; title('Subplot 1 Title'); subplot(2,1,2);

x=-pi:pi/20:pi; y=cos(x); plot(x,y); grid on; title('Subplot 2

Title');

-4

-3-2-10

1234

-1-0.500.5

1Subplot 1 Title

-4

-3

-2

-1

1

2

3

4

-1-0.500.5

1Subplot 2 Title

3． Polar Plots exercise

g=0.5;

theta=0:pi/20:2*pi;

gain=2*g*(1+cos(theta)); polar(theta,gain,'r-');

title('\fontsize{20} \bfGain versus angle \theta');

-4

-3-2-101234

-1-0.500.5

1Subplot 1 Title

Gain versus angle θ

4. Assume that a,b,c, and d are defined, and evaluate the following expression. a=20; b=-2; c=0; d=1;

(1) a>b; (2) b>d; (3) a>b&c>d; (4) a==b; (5) a&b>c; 6) ~~b; a=2; b=[1 –2;-0 10]; c=[0 1;2 0]; d=[-2 1 2;0 1 0]; (7) ~(a>b) (8) a>c&b>c (9) c<=d a=2; b=3; c=10; d=0;

(10) a*b^2>a*c (11) d|b>a (12) (d|b)>a a=20; b=-2; c=0; d=’Test’;

(13) isinf(a/b) (14) isinf(a/c) (15) a>b&ischar(d) (16) isempty(c)

5. Write a Matlab program to solve the function 1

()ln

1y x x

=-, where x is a number <1. Use an if structure to verify that the value passed to the program is legal. If the value of x is legal, caculate y(x). If not ,write a suitable error message and quit.

PART II: (需提交实验报告)

1.

3. Write out m. file and plot the figures with grids

Assume that the complex function f(t) is defined by the equation f(t)=(0.5-0.25i)t-1.0

Plot the amplitude and phase of function for 0 4.t ≤≤ t=0:0.1:4;

ft=(0.5-0.25i)*t-1.0; amp=abs(ft); phase=angle(ft);

subplot(211); plot(amp,t);

title('ft 的幅度相位图')

xlabel('t axis'); ylabel('ft 幅度 axis'); grid on subplot(212); plot(phase,t);

xlabel('t axis'); ylabel('ft 相位 axis'); grid on

0.4

0.60.81 1.2 1.4 1.6 1.8

0123

4

ft 的幅度相位图

t axis

f t 幅度 a x i s

-4

-3-2-1

01234

0123

4

t axis

f t 相位 a x i s

4. Write the Matlab statements required to calculate y(t) from the equation

22350

()350t t y t t t -+≥?=?

+

for value of t between –9 and 9 in steps of 0.5. Use loops and branches to perform this calculation.

5. Home work: 需提交实验报告

See the textbook page 73 : 2.12, 2.13, 2.15, 2.16

MATLAB实验三

一、实验目的：

1. Learn to design loop statements program

2. Be familiar with relational and logical operators

3. Practice 2D plotting

二、实验内容：

PART I: (选择练习，不需提交实验报告)

https://www.360docs.net/doc/c215377379.html,pare the 3 approaches follows (Loops and Vectorization)

%A. Perform calculation by For Loop with pre-initialize array

tic;

square=zeros(1,10000) %pre-initialize array

for ii=1:10000

square(ii)=ii^2;

square_root(ii)=ii^(1/2);

cube_root(ii)=ii^(1/3);

end

toc; t1=toc

%B. Perform calculation by For Loop without pre-initialize array

tic;

for ii=1:10000

square(ii)=ii^2;

square_root(ii)=ii^(1/2);

cube_root(ii)=ii^(1/3);

end

toc; t2=toc

%C. Perform calculation with vectors

tic;

ii=1:10000

square(ii)=ii.^2;

square_root(ii)=ii.^(1/2);

cube_root(ii)=ii.^(1/3);

end

toc; t3=toc

Home work: 需提交实验报告

Page 131: 3.2, 3.4, 3.8, 3.11

MATLAB 实验四

一、实验目的：

1. Learn to write MATLAB functions

2. Be familiar with complex data and character data

3. Practice 2D plotting

二、实验内容： Derivative of a function.

1. Write three Matlab functions to calculate the hyperbolic sine, cosine, and tangent functions:

s i n h (

),c o s h (

),t a n ()

2

2

x x

x x

x x

x x e e e e e e h e e

-----+-=+

then plot the shapes of hyperbolic sine, cosine, and tangent functions on one figure,

55x -≤≤.

2. Write a program use the function 32()552f x x x x =-+- and plot the line,

and search for the minimum and maximum in 200 steps over the range of 13x -≤≤, mark the minimum and maximum on the line figure.

3. Write a function to calculate the distance between two points 11(,)x y and 22(,)x y , that the points

should be given by ‘input ’ function.

4. Write a function complex_to that accept a complex number var , and returns two output

arguments containing the magnitude mag and angle theta of the complex number. The output angle should be in degrees.

Write another function polar_to_complex that accepts two input arguments containing the magnitude mag and angle theta of the complex number in degrees, and returns the complex number var .

4. Write a program that accepts a series of strings from a user with the input function, sorts

the strings into ascending order, and prints them out.

Derivative of a function. The derivative of a continuous function f(x) is defined by the equation

x

x f x x f x f dx d x ?-?+=→?)()(lim )(0

In a sampled function, this definition becomes

x

x f x f x f i i ?-=

+)

()()(1'

where

Home works: 需提交实验报告

Page 181: 4.15,4.16, 4.19, 4.20, 4.27

MATLAB 实验五

一、实验目的：

1. Be familiar with complex data and character data 2． ractice 2D plotting

二、实验内容：

1. Write three Matlab functions to calculate the hyperbolic sine, cosine, and tangent functions:

s i n h (

),c o s h (

),t a n ()

2

2

x x

x x

x x

x x e e e e e e h e e

-----+-=+

then plot the shapes of hyperbolic sine, cosine, and tangent functions on one figure,

55x -≤≤.

2. Write a program use the function 32()552f x x x x =-+- and plot the line,

and search for the minimum and maximum in 200 steps over the range of 13x -≤≤, mark the minimum and maximum on the line figure.

3. Write a function to calculate the distance between two points 11(,)x y and 22(,)x y , that the points

should be given by ‘input ’ function.

4. Write a function complex_to that accept a complex number var , and returns two output

arguments containing the magnitude mag and angle theta of the complex number. The output angle should be in degrees.

Write another function polar_to_complex that accepts two input arguments containing the magnitude mag and angle theta of the complex number in degrees, and returns the complex number var .

5. Write a program that accepts a series of strings from a user with the input function, sorts the strings into ascending order, and prints them out. Home works: 需提交实验报告

Page 181: 5.8, 5.10, 5.14, 5.16, 5.17, 5.22, 5.24

MATLAB 实验六

一、实验目的：

1. Practice 2D plotting and 3D plotting

2. Learn to use fplot function

3. Be familiar with cell arrays and structure arrays

二、实验内容：

1. Give the 3D plot figure of 0.30.1()sin(3),()cos()t t x t e t y t e t --== use plot3 function,

020

x ≤≤, and grid on, line width is 3.0.

-20

20

x

y

z

t=-10:0.1:10;

x=exp(-0.3*t).*sin(3*t); y=exp(-0.1*t).*cos(t);

plot3(x,y,t,'linewidth',3);grid on xlabel('x') ylabel('y') zlabel('z')

2. Plot the function sin x y e x -=, 02x ≤≤, step 0.1. Create the following plot types: (a) stem

plot; (b) stair plot; (c) bar plot; (d) compass plot. x=0:0.1:2; y=exp(-x).*sin(x); subplot(221);stem(x,y); subplot(222);stairs(x,y); subplot(223);bar(x,y); subplot(224);compass(x,y);

00.51 1.52

0.10.20.30.4-1

1

2

3

00.10.20.3

0.4

3. Plot the function ()f x = over the range 0.110.0x ≤≤ using function fplot, and

grid on. function y=myfun(x) y=1./sqrt(x); fh=@myfun; fplot(fh,[0 20]);grid on

2

4

6

8

10

12

14

16

18

20

00.511.522.533.544.5

5

4. Create a cell arrays:

a={[1 4 3;0 5 8;7 2 9],'anne smith',3+7i;5-2i,-3.14,[];3.14,[],5}

a =

[3x3 double] 'anne smith' [3.0000+ 7.0000i] [5.0000- 2.0000i] [ -3.1400] [] [ 3.1400] [] [ 5]

5. Create a structure arrays and to calculate the mean billing of three patients:

Home works:需提交实验报告

Page 285:6.4, 6.9, 6.11

MATLAB实验七

一、实验目的：

Be familiar with Input/Output functions

二、实验内容：

1. Write a m-file. The m-file creates an array containing 150

?random values, sorts the array into ascending order, opens a user-specified file for writing only, then writes the array to disk in 32-bit floating-point format , and close the file. It then opens the file and read the data back into 510

?array.

?square matrix, then import the array use

2. Edit a file as data4_4.txt that contains 44

uiimport function, and calculate the inverse of the square matrix .

3. Write a program to read a set of integers from an input data file, and locate the largest and

smallest values within the data file. Print out the largest and smallest values, together with the lines on which they were found in one figure.

MATLAB实验八

一、实验目的：

1.Be familiar with handle graphics

2. Learn about GUI

3. Learn to design PID control system with Simulink

二、实验内容：

1. Second Order System GUI Design

% SOSystem_GUI.m

set(gcf,'defaultuicontrolunits','normalized');

set(gcf,'defaultuicontrolfontsize',12);

str=' Second Order System Step response'; set(gcf,'name',str);

h_axes=axes('position',[0.05,0.2,0.6,0.7]); set(h_axes,'xlim',[0,15]);