MATLAB实验报告 彭羊平Word文件下载.docx
《MATLAB实验报告 彭羊平Word文件下载.docx》由会员分享，可在线阅读，更多相关《MATLAB实验报告 彭羊平Word文件下载.docx（29页珍藏版）》请在冰豆网上搜索。
4、程序如下：
A=fix（100*rand（8））+5
B=A>
100;
k=find（any（B'
））
A（k,:
）=0
结果如下:
实验二：
matlab绘图
（1）
1、了解MATLAB图形系统和各种图形函数。
2、熟悉MATLAB的基本图形操作，具备MATLAB画图能力。
3、熟悉各种数学函数，并通过图形函数画出。
1、画出常用三角函数图形。
2、绘制出多峰函数peaks和三角函数的多条曲线。
3、有一位研究生，一年中平均每月的费用为生活费190元、资料费33元、电话费45元、购买衣服42元及其它费用45元。
请以饼图表示出他每月的消费比例，并分离出表示资料费用的的切片。
请给出图中每一块加以标注。
1、试验程序：
x=0:
pi50:
2*pi;
k=[1265176101];
x（k）=[];
figure
（1）;
subplot（2,2,1）,plot（x,sin（x））,gridon;
subplot（2,2,2）,plot（x,cos（x））,gridon;
subplot（2,2,3）,plot（x,tan（x））,gridon;
subplot（2,2,4）,plot（x,cot（x））,gridon;
实验结果：
2、实验程序：
（1）[x,y]=meshgrid（-3:
.15:
3）;
z=peaks（x,y）
x1=x（1,:
plot（x1,z）,gridon
实验结果:
（2）x=0:
y1=sin（x）;
y2=cos（x）;
figure
（1）
plot（x,y1,'
k-'
x,y2,'
k--'
）,gridon
text（pi,0.05,'
\leftarrowsin（\alpha）'
text（pi4-0.05,0.05,'
cos（\alpha）\rightarrow'
title（'
sin（\alpha）andcos（\alpha）'
xlabel（'
\alpha'
x=[19033454245];
explode=[01000];
colormap（40）
y=randn（40）
z=（x-2）.^2+（y-1.2）.^2
plot3（x,y,z）,gridon
网格曲线
实验程序：
x=randn（40）
mesh（x,y,z）
曲面
[x,y]=meshgrid（-3:
.125:
z=x.^2+y.^2+sin（x*y）;
c=ones（size（x））;
surfc（x,y,z,c）,gridon
等高线
contour3（x,y,z,30）
view（-15,25）
实验四：
matlab程序设计
（1）
1、实验目的：
1、初步学会MATLAB编程，对其有所了解。
2.、熟悉冒泡和选择排序的编程方法。
1、实现冒泡排序。
2、实现选择排序。
1、实验程序：
a=[79955288179727258773]
n=length（a）
fori=1:
n-1
flag=-1;
forj=n:
-1:
i+1
if（a（j）<
a（j-1））
temp=a（j）;
a（j）=a（j-1）;
a（j-1）=temp;
flag=0;
b=a;
end
b
disp（['
循环次数为：
'
num2str（i）]）;
functionb=xuanze（a）
n=length（a）;
min=i;
forj=i+1:
n
a（min））
min=j;
if（min~=1）
temp=a（min）;
a（min）=a（i）;
a（i）=temp;
b=a;
实验五：
matlab程序设计
（2）
一、实验目的：
1、进一步熟悉MATLAB编程语言，熟悉程序设计技术。
2、了解MATLAB的程序调试功能，通过编程熟悉。
2、实验内容：
1、企业发放的奖金按个人完成的利润（I）提成，分段提成比例K如下：
K=10%I<
=10万元
5%10<
I<
=20万元
2%20<
=40万元
1%I>
40万元
即如王某完成25万元利润时，个人可得
y=10*10%+10*5%+5*2%（万元）
据此编写程序，求企业职工的奖金。
2、有一分数序列
求前15项的和。
（扩展求其前n项和）
3、有n个人围成一个圈，按序列编号，从第一个人开始报数，数到m时该人退出，并且下一个人从1开始重新报数，求出出圈人的顺序（n>
m,l例如n=20,m=7）。
3、实验结果：
（函数文件）
functiony=jiangjin（x）
y=0;
if（x>
=0&
&
x<
=10）
y=x*0.1;
elseif（10<
x&
=20）
y=10*0.1+（x-10）*0.05;
elseif（20<
=40）
y=10*0.1+10*0.05+（x-20）*0.02;
elseif（x>
40）
y=10*0.1+10*0.05+20*0.02+（x-40）*0.01;
elsey='
error'
;
functiony=sum15（x）
n11=2;
n1=1;
A1=n11.n1;
fori=1:
x
y=y+a1;
A1=（n11+n1）.n11;
3、实验程序：
functiona=yuesefuhuan（m,n）
b=1:
n;
k=0;
%k为数到几了。
i=1;
m1=1;
while（m1<
=n）
if（b（i）~=0）
k=k+1;
end
if（k==m）
k=0;
a（m1）=i;
m1=m1+1;
b（i）=0;
if（i==n）
i=0;
i=i+1;
实验六：
matlab基本应用
（1）
1.、了解MATLAB的基本应用领域：
线性代数、多项式和内插、数据分析与统计、泛函分析、常微分方程求解等。
2、熟悉线性代数方程和常微分方程的求解与编程。
3、了解数据拟合的方法与编程。
1、求解下列线性代数方程。
（1）x1+2x2+3x3=11
（2）3x1+x2+x4=2
2x1+2x2+5x3=126x2+7x3+3x4=4
3x1+5x2+x3=314x2+3x3=7
2x1-x2+2x3+6x4=8
2、通过测量得到一组数据
分别采用y1（t）=c1+c2e-t和y2（t）=d1+d2e-t进行拟合，并画出拟合曲线进行对比。
3、求解下列线性微分方程，并画出状态轨迹。
（1）5x1-5x2-6x3=x1'
（2）x1+2x2-3x3+x4=x1'
3x1-2x2+5x3=x2'
3x1+x3-2x4=x2'
2x1-x2-4x3=x3'
x1-2x2+5x4=x3'
x（0）=[1,-4,5]'
2x1+3x2+x4=x4'
x（0）=[]'
三、实验结果：
a=[1,2,3;
2,2,5;
3,5,1];
b=[]'
x=a\b
a1=[;
];
b1=[]'
x=a1\b1
t=[1:
10]'
y=[4.842,4.362,3.754,3.368,3.169,3.083,3.034,3.016,3.012,3.005]'
A=[ones（size（t））exp（-t）];
A1=[ones（size（t））t.*exp（-t）];
C=A\y;
C1=A1\y;
T=[0:
.1:
11]'
Y=[ones（size（T））exp（-T）]*C;
Y1=[ones（size（T））T.*exp（-T）]*C1;
subplot（2,1,1）,plot（T,Y,'
-'
t,y,'
o'
）,title（'
y1（t）=c1+c2exp（-t）拟合'
）;
subplot（2,1,2）,plot（T,Y1,'
y2（t）=d1+d2t*exp（-t）拟合'
functiondy=diffss（t,y）
dy=[5*y
（1）-5*y
（2）-6*y（3）;
3*y
（1）-2*y
（2）+5*y（3）;
2*y
（1）-y
（2）-4*y（3）;
X0=[1;
-4;
5];
tspan=[0,100];
[T,X]=ode45（'
diffss'
tspan,X0）;
subplot（3,1,1）,plot（T,X（:
1））,title（'
x_{1}'
）,gridon;
subplot（3,1,2）,plot（T,X（:
2））,title（'
x_{2}'
subplot（3,1,3）,plot（T,X（:
3））,title（'
x_{3}'
figure
（2）;
plot（X（:
1）,X（:
2）,X（:
3）,title（'
系统轨迹'
））,gridon;
实验七：
matlab基本应用
（2）
1.、进一步了解MATLAB的基本应用领域，熟悉高阶微分方程求解和编程等。
2、熟悉定积分的求解与编程。
3、了解求信号DFT的方法与编程。
1、求解范德堡方程
。
其中，u>
0为标量。
2、计算下列定积分。
3、微分方程组
当t=0时，x1（0）=1,x2（0）=-0.5,求微分方程上的解，并画出x1-x2的系统轨迹。
functiondy=diffss111（t,y）
globalmu
dy=[y
（2）;
mu*（1-y
（1）*y
（1））*y
（2）-y
（1）];
globalmu;
mu=1;
X0=[2;
0];
tspan=[0,20];
diffss111'
subplot（2,1,1）;
plot（T,X（:
1）,'
T,X（:
2）,'
--'
范德堡方程的解,u=1'
ylabel（'
Soutionx_{1}andx_{2}'
legend（'
'
tspan=[0,3000];
mu=1000;
[T,Y]=ode15s（'
subplot（2,1,2）;
plot（T,Y（:
范德堡方程的解,u=1000'
时间t'
）,ylabel（'
Soutionx_{1}'
2、实验程序：
.001:
2;
y1=exp（-2*x）;
z1=trapz（x,y1）
y2=exp（2*x）;
z2=trapz（x,y2）
x1=-1:
y3=x.*x-3*x+0.5;
z3=trapz（x1,y3）
functiondy=difffss（t,y）
dy=[0.5-y
（1）;
y
（1）-4*y
（2）];
-0.5];
tspan=[0,25];
difffss'
subplot（3,1,3）,plot（T,（X（:
1）-X（:
1）））,title（'
x_{1}-x_{2}'
实验结果如下：
实验八：
数据阵列与结构
1、熟悉MATLAB中数据阵列类型，如多维阵列，结构阵列和单元阵列的产生，结构及其使用。
2、了解阵列类型与结构的简单编程。
1、利用下标建立多维阵列。
2、建立一个单元阵列并存放一些数据。
3、在一个5*5*2.5m3的房间里，按每间隔50cm的栅格测定各点的温度（可随机产生），建立多维阵列表示，并求出高度从0~2.5m每隔50cm的平面上的平均温度。
a（:
:
2）=[235;
790;
785]
A=cell（2,2）
A{1,1}=[126;
094;
234]
A{1,2}='
AnneSmith'
A{2,1}=3+7i;
A{2,2}=-pi:
pi10:
pi;
celldisp（A）
.5:
5;
y=0:
z=0:
2.5;
Temp=fix（20*rand（length（x）,length（y）,length（z）））+9
fork=1:
6
Temaverage（k）=sum（sum（Temp（:
k）））（length（x）*length（y））;
Temaverage
实验结果为：
实验九：
综合练习
1、进一步熟悉MATLAB编程语言，熟悉程序设计技术。
2、了解MATLAB的基本应用领域。
3、进一步掌握图形的绘制方法。
1、产生一个信号，
其中，
，这一信号表示被噪声污染的信号，设计程序求其DFT，并绘图表示。
2、画出各种大小和形状的球、柱体。
t=0:
0.001:
0.6;
x=3*sin（2*pi*20*t）+10*sin（2*pi*200*t+pi4）+10*randn（1,length（t））;
Y=fft（x,512）;
P=Y.*conj（Y）;
f=1000*（0:
255）512;
plot（f,P（1:
256））
球
sphere
axissquare
[x,y,z]=sphere
x=0.5*x
y=0.5*y
z=0.5*z
surf（x,y,z）
surf（x+1,y+1,z+1
柱体：
t=0:
[x,y,z]=cylinder（2+cos（t））;
surf（x,y,z）,axissquare
升级会员 