信号处理仿真(MATLAB)实验系统仿真(MATLAB)实验2

信号处理仿真(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

2． Give the answer of the following questions for the array 1）What is the size of array1?

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

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

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

3． Give the answer of the following commad

1) a=1:2:5; 2) b=[a’ a’ a’]; 3) c=b(1:2:3,1:2:3);

4) d=a+b(2,:) 5) w=[zeros(1,3) ones(3,1)’ 3:5’]

4． Give the answer of the sub-arrays

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

5． Give the answer of the following operations

1) a+b 2) a*d 3) a.*d 4) a*c

5) a.*c 6)a\b 7) a.\b 8)a.^b

6．Solve the following system of simultaneous equations for x 7．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)

end

hold off;

8．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;

9． 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;

10．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;

MATLAB实验二

一、实验目的：

1. Learn to design branch and loop statements program

2. Be familiar with relational and logical operators

3. Practice 2D plotting

二、实验内容：

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

1. Hold command exercise

n x=-pi:pi/20:pi;

n y1=sin(x); y2=cos(x); plot(x,y1, 'b-'); hold on;

n plot(x,y2, 'k--'); hold off;

n legend ('sinx', 'cosx')

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');

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. Compare the 3 approaches follows (Loops and Vectorization)

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

tic;

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

n for ii=1:10000

n square(ii)=ii^2;

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

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

n end

n toc; t1=toc

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

n tic;

n for ii=1:10000

n square(ii)=ii^2;

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

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

n end

n toc; t2=toc

%C. Perform calculation with vectors

n tic;

ii=1:10000

n square(ii)=ii.^2;

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

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

n end

n toc; n t3=toc

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

1. 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)

2. Write a Matlab program to solve the function, 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.

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

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

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

5. Write an m.file to evalue the equation for all values of x between 0.1 and 3,

in steps of 0.1. Do this twice, once with a for loop and once with vectors.

Plot the resulting function using a 4.0 thick dashed blue line.

MATLAB实验三

一、实验目的：

1. Learn to write MATLAB functions

2. Be familiar with complex data and character data

3. Practice 2D plotting

二、实验内容：

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

tangent functions:

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

2. Write a program use the function and plot the line,

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

3. Write a function to calculate the distance between two points and , 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.

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 use plot3 function, , and grid on, linewidth is

3.0.

2. Plot the function , , step 0.1. Create the following plot types: (a) stem plot;

(b) stair plot; (c) bar plot; (d) compass plot.

3. Plot the function over the range using function fplot, and grid on.

4. Create a cell arrays:

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

patients:

MATLAB实验五

一、实验目的：

Be familiar with Input/Output functions

二、实验内容：

1. Write a m-file. The m-file creates an array containing 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 array.

2. Edit a file as data4_4.txt that contains square matrix, then import the

array use 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.