最新版MATLAB实验报告 彭羊平Word文档格式.docx
- 文档编号:18138358
- 上传时间:2022-12-13
- 格式:DOCX
- 页数:30
- 大小:549.68KB
最新版MATLAB实验报告 彭羊平Word文档格式.docx
《最新版MATLAB实验报告 彭羊平Word文档格式.docx》由会员分享,可在线阅读,更多相关《最新版MATLAB实验报告 彭羊平Word文档格式.docx(30页珍藏版)》请在冰豆网上搜索。
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:
pi/50:
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(pi/4-0.05,0.05,'
cos(\alpha)\rightarrow'
title('
sin(\alpha)andcos(\alpha)'
xlabel('
\alpha'
x=[19033454245];
explode=[01000];
colormaphsv;
pie(x,explode);
研究生每月消费比例'
);
实验三:
matlab绘图
(2)
1、实验目的:
1、熟练掌握MATLAB二维曲线的绘制方法;
2、掌握网格数据的生成方法;
3、掌握三维图的绘制方法;
4、理解各种特殊图形的绘制方法。
1、画出下列函数的三维曲线和网格曲线。
2、画出下列函数的曲面及等高线图。
3、画出各种大小和形状的球、柱体。
1、实验程序:
三维曲线
x=randn(40)
y=randn(40)
z=(x-2).^2+(y-1.2).^2
plot3(x,y,z),gridon
网格曲线
实验程序:
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、初步学会MATLAB编程,对其有所了解。
2.、熟悉冒泡和选择排序的编程方法。
1、实现冒泡排序。
2、实现选择排序。
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:
pi/10:
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+pi/4)+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
holdon
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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 最新版MATLAB实验报告 彭羊平 最新版 MATLAB 实验 报告