MATLAB实验报告 彭羊平.docx

1、MATLAB实验报告 彭羊平(此文档为word格式，下载后您可任意编辑修改！)班级： 通信班 姓名： 彭羊平 学号： 033 实验 一 ： matlab的基本操作一、 实验目的： 1、了解MATLAB的集成环境，熟悉其基本操作。 2、了解MATLAB的基础知识，包括矩阵的建立、简单操作、逻辑操作和关系运算。 3、熟悉基本的数学函数和逻辑函数。 4、在命令窗口输入命令完成一些简单的功能，为MATLAB程序设计奠定基础。二、 实验内容：1、利用diag等函数产生下列矩阵： 然后利用reshape函数将它们变换成行向量。2、利用rand函数产生（0,1）间均匀分布的10*10随机矩阵A，然后统计A中

2、大于等 于0.6的元素的个数。 3、有一矩阵A，找出矩阵中值大于1的元素，并将它们重新排列成列向量B。4、在一测量矩阵A(6*6)，存在有奇异值（假设大于100的值为奇异值），编程实现删去奇异值所在的行。 三、 实验结果： 1、程序如下： a=fliplr(diag(8,-7,2)+diag(5,3,-1) a=reshape(a,1,9) b=diag(2,5,8)+diag(4,2)+diag(7,-2) b=reshape(b,1,9) 结果如下：2、程序如下： A=rand(10) A=A0.6 sum(sum(A) 结果如下：3、程序如下： A=fix(5*rand(5) a=A1;

3、 K=find(a); for k=1:length(K) B(k)=A(K(k); end B=B 结果如下： 4、程序如下：A=fix(100*rand(8)+5B=A100;k=find(any(B)A(k,:)=0结果如下:班级： 通信班 姓名： 彭羊平 学号： 033 实验 二 ： matlab绘图（1）一、 实验目的：1、了解MATLAB图形系统和各种图形函数。2、熟悉MATLAB的基本图形操作，具备MATLAB画图能力。3、熟悉各种数学函数，并通过图形函数画出。二、 实验内容：1、画出常用三角函数图形。2、绘制出多峰函数peaks和三角函数的多条曲线。3、有一位研究生，一年中平均

4、每月的费用为生活费190元、资料费33元、电话费45 元、购买衣服42元及其它费用45元。请以饼图表示出他每月的消费比例，并分离出表示资料费用的的切片。请给出图中每一块加以标注。三、 实验结果：1、试验程序：x=0:pi50:2*pi;k=1 26 51 76 101;x(k)=;figure(1);subplot(2,2,1),plot(x,sin(x),grid on;subplot(2,2,2),plot(x,cos(x),grid on;subplot(2,2,3),plot(x,tan(x),grid on;subplot(2,2,4),plot(x,cot(x),grid on;实

5、验结果：2、实验程序：（1）x,y=meshgrid(-3:.15:3);z=peaks(x,y)x1=x(1,:)plot(x1,z),grid on 实验结果:（2）x=0:pi50:2*pi;y1=sin(x);y2=cos(x);figure(1)plot(x,y1,k-,x,y2,k-),grid ontext(pi,0.05,leftarrow sin(alpha)text(pi4-0.05,0.05,cos(alpha)rightarrow)title(sin(alpha) and cos(alpha)xlabel(alpha)2、实验程序： x=190 33 45 42 45;

6、explode=0 1 0 0 0;figure(1);colormap (40)y=randn(40)z=(x-2).2+(y-1.2).2plot3(x,y,z),grid on实验结果：网格曲线实验程序：x=randn(40)y=randn(40)z=(x-2).2+(y-1.2).2mesh(x,y,z)实验结果：2、实验程序：曲面x,y=meshgrid(-3:.125:3);z=x.2+y.2+sin(x*y);c=ones(size(x);surfc(x,y,z,c),grid on实验结果：等高线x,y=meshgrid(-3:.125:3);z=x.2+y.2+sin(x*y

7、);contour3(x,y,z,30)view(-15,25)实验结果：班级： 通信班 姓名： 彭羊平 学号： 033 实验 四 ： matlab程序设计（1）1、实验目的：1、初步学会MATLAB编程，对其有所了解。2.、熟悉冒泡和选择排序的编程方法。二、 实验内容： 1、 实现冒泡排序。 2、实现选择排序。三、 实验结果：1、实验程序：a=79 95 52 88 17 97 27 25 87 73n=length(a)for i=1:n-1 flag=-1; for j=n:-1:i+1 if(a(j)a(j-1) temp=a(j); a(j)=a(j-1); a(j-1)=temp;

8、 flag=0; end end b=a; endbdisp(循环次数为：,num2str(i);实验结果： 2、实验程序：function b=xuanze(a)n=length(a);for i=1:n-1 min=i; for j=i+1:n if(a(j)a(min) min=j; end end if(min=1) temp=a(min); a(min)=a(i); a(i)=temp; endendb=a;实验结果：班级： 通信班 姓名： 彭羊平 学号： 033 实验 五： matlab程序设计（2）一、实验目的：1、进一步熟悉MATLAB编程语言，熟悉程序设计技术。2、了解MAT

9、LAB的程序调试功能，通过编程熟悉。2、实验内容：1、企业发放的奖金按个人完成的利润（I）提成，分段提成比例K如下： K= 10% I=10万元 5% 10I=20万元 2% 20I40万元 即如王某完成25万元利润时，个人可得 y=10*10%+10*5%+5*2%(万元) 据此编写程序，求企业职工的奖金。2、有一分数序列 求前15项的和。（扩展求其前n项和） 3、有n个人围成一个圈，按序列编号，从第一个人开始报数，数到m时该人退出，并且下一个人从1开始重新报数，求出出圈人的顺序（nm,l例如n=20,m=7）。3、实验结果： 1、实验程序：（函数文件） function y=jiangji

10、n(x) y=0; if (x=0 & x=10) y=x*0.1; elseif (10x & x=20) y=10*0.1+(x-10)*0.05; elseif (20x & x40) y=10*0.1+10*0.05+20*0.02+(x-40)*0.01; else y=error; end实验结果： 2、实验程序： function y=sum15(x) y=0; n11=2; n1=1; A1=n11.n1; for i=1:x y=y+a1; A1=(n11+n1).n11; end实验结果：3、实验程序：function a=yuesefuhuan(m,n)b=1:n;k=0

11、; %k为数到几了。i=1;m1=1;while (m10为标量。2、计算下列定积分。 3、微分方程组 当t=0时，x1(0)=1,x2(0)=-0.5,求微分方程上的解，并画出 x1-x2的系统轨迹。3、实验结果：1、实验程序：function dy=diffss111(t,y)global mudy=y(2);mu*(1-y(1)*y(1)*y(2)-y(1);global mu;mu=1;X0=2;0;tspan=0,20;T,X=ode45(diffss111,tspan,X0);figure(1);subplot(2,1,1);plot(T,X(:,1),-,T,X(:,2),-),

12、grid on;title(范德堡方程的解, u=1)ylabel(Soution x_1 and x_2)legend(x_1,x_2)tspan=0,3000;mu=1000;T,Y=ode15s(diffss111,tspan,X0);subplot(2,1,2);plot(T,Y(:,1),title(x_1),grid on;title(范德堡方程的解, u=1000),grid on;xlabel(时间t),ylabel(Soution x_1);实验结果： 2、实验程序：x=0:.001:2;y1=exp(-2*x);z1=trapz(x,y1)y2=exp(2*x);z2=tr

13、apz(x,y2)x1=-1:.001:1;y3=x.*x-3*x+0.5;z3=trapz(x1,y3)实验结果：3、实验程序：function dy=difffss(t,y)dy=0.5-y(1);y(1)-4*y(2);X0=1;-0.5;tspan=0,25;T,X=ode45(difffss,tspan,X0);figure(1);subplot(3,1,1),plot(T,X(:,1),title(x_1),grid on;subplot(3,1,2),plot(T,X(:,1),title(x_2),grid on;subplot(3,1,3),plot(T,(X(:,1)-X(

14、:,1),title(x_1-x_2),grid on;实验结果如下：班级： 通信班 姓名： 彭羊平 学号： 033 实验 八： 数据阵列与结构1、实验目的： 1、熟悉MATLAB中数据阵列类型，如多维阵列，结构阵列和单元阵列的产生，结构及其使用。 2、了解阵列类型与结构的简单编程。2、实验内容： 1、利用下标建立多维阵列。 2、建立一个单元阵列并存放一些数据。 3、在一个5*5*2.5m3的房间里，按每间隔50cm的栅格测定各点的温度（可随机产生），建立多维阵列表示，并求出高度从02.5m每隔50cm的平面上的平均温度。 3、实验结果：1、实验程序：a(:,:,2)=2 3 5;7 9 0;

15、 7 8 5实验结果：2、实验程序：A=cell(2,2)A1,1=1 2 6 ;0 9 4 ;2 3 4A1,2=Anne Smith;A2,1=3+7i;A2,2=-pi:pi10:pi;实验结果：celldisp(A)3、实验程序：x=0:.5:5;y=0:.5:5;z=0:.5:2.5;Temp=fix(20*rand(length(x),length(y),length(z)+9for k=1:6Temaverage(k)=sum(sum(Temp(:,:,k)(length(x)*length(y);endTemaverage实验结果为： 班级： 通信班 姓名： 彭羊平 学号： 0

16、33 实验 九： 综合练习 1、实验目的： 1、进一步熟悉MATLAB编程语言，熟悉程序设计技术。 2、了解MATLAB的基本应用领域。 3、进一步掌握图形的绘制方法。2、实验内容： 1、产生一个信号，其中，这一信号表示被噪声污染的信号，设计程序求其DFT，并绘图表示。 2、画出各种大小和形状的球、柱体。3、实验结果：1、实验程序：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)实验结果： 2、实验程序：球sphere axis squarex,y,z=spherex=0.5*xy=0.5*yz=0.5*z surf(x,y,z)surf(x+1,y+1,z+1柱体： t=0:pi10:2*pi;x,y,z=cylinder(2+cos(t);surf(x,y,z),axis square