北科大Matlab_数学实验报告1~6次(全)
《数学实验》报告
实验名称 Matlab 基础知识
学院
专业班级
姓名
学号
2014年 6月
一、【实验目的】
1.认识熟悉Matlab这一软件，并在此基础上学会基本操作。
2.掌握Matlab基本操作和常用命令。
3.了解Matlab常用函数，运算符和表达式。
4.掌握Matlab工作方式和M文件的相关知识。
5.学会Matlab中矩阵和数组的运算。
二、【实验任务】
P16 第4题
编写函数文件，计算
1!
n
k
k =
∑，并求出当k=20时表达式的值。P27第2题
矩阵A=
123
456
789
??
??
??
??
??
，B=
468
556
322
??
??
??
??
??
，计算A*B,A.*B,并比较两者的区别。
P27第3题
已知矩阵A=
52
91
??
??
??
，B=
12
92
??
??
??
，做简单的关系运算A>B,A==B,A (A==B)&(A P34 第1题 用 111 1 4357 π =-+-+……公式求π的近似值，直到某一项的绝对值小于-6 10为止。 三、【实验程序】 P16 第4题 function sum=jiecheng(n) sum=0; y=1; for k=1:n for i=1:k y=y*i; end sum=sum+y; end sum P27第2题 >>A=[1 2 3;4 5 6;7 8 9] >>B=[4 6 8;5 5 6;3 2 2] >>A*B P27第3题 >> A=[5 2;9 1];B=[1 2;9 2]; >>A>B >>A==B >>A >> (A==B)&(A >> (A==B)&(A>B) P34 第1题 t=1; pi=0; n=1; s=1; while abs(t)>=1e-6 pi=pi+t; n=n+2; s=-s; t=s/n; end pi=4*pi; 四、【实验结果】 P16 第4题 P27第2题 两者的区别：A*B是按正规算法进行矩阵的计算, A.*B是对应元素相乘。P27第3题 P34 第1题 >> pi pi= 3.141590653589692e+000 五、【实验总结】 这次实验是第一次接触Matlab这个软件，所以有些生疏，花的时间也比较多，但功夫不怕有心人，而且当一个程序弄出来后感觉也特别开心，以后再继续努力学习。 了解并掌握matlab的基本绘图 二、【实验任务】 P79页1,3,5题 三、【实验程序】 1. clf; x=0:pi/50:4*pi; y1=exp(x/3).*sin(3*x);y2=exp(x/3);y3=-exp(x/3); plot(x,y1,'b*',x,y2,'r-.',x,y3,'r-.'),grid on legend('y1=exp(x/3).*sin(3*x)','y2=+-exp(x/3)') 3. clf; x1=-pi:pi/50:pi; x2=pi:pi/50:4*pi; x3=1:0.1:8; y1=x1.*cos(x1); y2=x2.*tan(x2.^(-1)).*sin(x2.^3); y3=exp(x3.^(-1)).*sin(x3); subplot(2,2,1),plot(x1,y1,'m.'),grid on,title('y=x*cosx') xlabel('x?á'),ylabel('y?á') gtext('y=x*cosx'),legend('y=x*cosx') subplot(2,2,2),plot(x2,y2,'r*'),grid on,title('y=x*tan(1/x)*sin(x^3)') xlabel('x?á'),ylabel('y?á') gtext('y=x*tan(1/x)*sin(x^3)'),legend('y=x*tan(1/x)*sin(x^3)') subplot(2,2,3),plot(x3,y3,'bp'),grid on,title('y=e(1/x3)*sinx') xlabel('x?á'),ylabel('y?á') gtext('y=e(1/x3)*sinx'),legend('y=e(1/x3)*sinx') 5. t=0:pi/50:20*pi; x=t.*cos(t*pi/6); y=t.*sin(t*pi/6); z=2*t; plot3(x,y,z) 1． 3. 5. 五、【实验总结】 一、 【实验目的】 1. 学会用Matlab 进行三维的曲线绘图； 2. 掌握绘图的基本指令和参数设置 二、 【实验任务】 P79 习题5 绘制圆锥螺线的图像并加标注，圆锥螺线的参数方程为； ??? ? ? ? ??? ===t z t t y t t x 26s i n 6c o s ππ )200(π≤≤t 。 P79 习题9 画三维曲线2)x,y (-2-x -52 2≤≤=y z 与平面z=3的交线。 三、 【实验程序】 习题5： clf; t=0:pi/100:20*pi; x=t.*cos(t.*pi/6); y=t.*sin(t.*pi/6); z=2*t; plot3(x,y,z) title('圆锥螺线') xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴') 习题9： clf; t=-2:0.1:2; [x,y]=meshgrid(t); z1=5-x.^2-y.^2; subplot(1,2,1),mesh(x,y,z1),title('曲面z1=5-x.^2-y.^2') z2=3*ones(size(x)); r0=abs(z1-z2)<=0.05; subplot(1,2,2),plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),'.') title('交线') 四、【实验结果】 习题5： 习题9： 五、【实验总结】 这次三维曲线（曲面）的绘制虽然不算复杂，但还是要注意一些细节，而且要注意弄懂其中的原因，不能硬套书上的，否则很容易不明道理的出错。 1. 学会用Matlab 练习使用矩阵的基本运算； 2. 掌握用Matlab 运用矩阵的特征值、特征向量、特征多项式； 3. 学会用Matlab 解线性方程组； 4. 掌握用Matlab 进行数值方法计算定积分 二、 【实验任务】 P114 习题12 随机输入一个六阶方阵，并求其转置、行列式、秩，以及行最简式。 P114 习题14 求矩阵??? ? ? ?????=2 1 11 2 11 1 2A 的特征多项式、特征值和特征向量。 P115 习题20 求下列线性方程组的通解： (1)??? ??=+=++=++0 -43-203-0 -243213214321x x x x x x x x x x x (2)??? ???? =+=+=+21 -32--13--0--432143214321x x x x x x x x x x x x P167 习题17 用三种方法求下列积分的数值解： (2)dx x x x ?+π 2 cos 1sin P167 习题18 用多种数值方法计算定积分?4 sin -11 π dx x ，并与精确值2进行比较，观察 习题12 >> A=[1 9 5 3 6 5;2 4 6 8 1 0;3 4 6 9 7 2;4 6 7 8 10 4;5 0 7 3 2 1;3 8 6 3 1 9] >> A' >> det(A) >> rank(A) >> rref(A) 习题14： >> B=[2 1 1;1 2 1;1 1 2] >> p=poly(B) >> [V D]=eig(B) 习题20： (1) >> A=[1 1 2 -4;-1 1 3 0;2 -3 4 -1] >> rref(A) (2) >> B=[1 -1 -1 1;1 -1 1 -3;1 -1 -2 3] >> rref(B) >> C=[1 -1 -1 1 0;1 -1 1 -3 1;1 -1 -2 3 -1/2] >> rref(C) 习题17： （2） function y=jifen(x) y=x.*sin(x)./(1+cos(x).^2); h=0.01; x=0:h:pi; y0=1+cos(x).^2; y1=x.*sin(x)./y0; t=length(x); s1=sum(y1(1:(t-1)))*h s2=sum(y1(2:t))*h s3=trapz(x,y1) s4=quad('jifen',0,pi) 习题18： function y=jifen(x) y=1./(1-sin(x)); h=0.01; x=0:h:pi/4; y=1./(1-sin(x)); t=length(x); format long s1=sum(y1(1:(t-1)))*h s2=sum(y1(2:t))*h s3=trapz(x,y) s4=quad('jifen',0,pi/4) format short u1=s1-sqrt(2) u2=s2-sqrt(2) u3=s3-sqrt(2) u4=s4-sqrt(2) 四、【实验结果】 习题12 习题14 习题20 （1） 原方程对应的同解方程组为：??? ? ? ????= ==43424 1256572553x x x x x x ，解得方程基础解系为：?????? ? ?? ? ??????????1256572553，所以方程组的通解为： ?????? ????????4321x x x x =???????? ? ???????????12565725531κ (2) 解对应的齐次方程组?? ?=+=4 3 4 21 2x x x x x ，可得一个基础解系：????? ???????=121 2ε 原方程组对应的同解方程组为：??? ??? ?+=++=2122143421x x x x x ，可找到一个特解为： ???? ???? ? ?????????=125125* η 因此，此方程组的通解为： ????? ?? ? ? ?????????+???? ????????=?????? ????????12512512 1214321k x x x x 习题17： （2） 习题18： 五、【实验总结】 在掌握线性代数相关运算和数值积分的理论基础上进行操作，学会了用Matlab相关指令和编程，并进行计算与误差分析，感觉原来很繁琐的计算用Matlab很方便！ 电子科技大学电子工程学院标准实验报告（实验）课程名称MATLAB与数值分析 学生姓名：李培睿 学号：2013020904026 指导教师：程建 一、实验名称 《MATLAB与数值分析》第一次上机实验 二、实验目的 1. 熟练掌握矩阵的生成、加、减、乘、除、转置、行列式、逆、范数等运算 操作。（用.m文件和Matlab函数编写一个对给定矩阵进行运算操作的程序） 2. 熟练掌握算术符号操作和基本运算操作，包括矩阵合并、向量合并、符号 转换、展开符号表达式、符号因式分解、符号表达式的化简、代数方程的符号解析解、特征多项式、函数的反函数、函数计算器、微积分、常微分方程的符号解、符号函数的画图等。（用.m文件编写进行符号因式分解和函数求反的程序） 3. 掌握Matlab函数的编写规范。 4、掌握Matlab常用的绘图处理操作，包括：基本平面图、图形注释命令、 三维曲线和面的填充、三维等高线等。（用.m文件编写在一个图形窗口上绘制正弦和余弦函数的图形，并给出充分的图形注释） 5. 熟练操作MATLAB软件平台，能利用M文件完成MATLAB的程序设计。 三、实验内容 1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。并以x， y为坐标显示图像 x(n+1) = a*x(n)-b*(y(n)-x(n)^2); y(n+1) = b*x(n)+a*(y(n)-x(n)^2) 2. 编程实现奥运5环图，允许用户输入环的直径。 3. 实现对输入任意长度向量元素的冒泡排序的升序排列。不允许使用sort 函数。 四、实验数据及结果分析 题目一： ①在Editor窗口编写函数代码如下： 实验一MATLAB操作基础 实验目的和要求： 1、熟悉MATLAB的操作环境及基本操作方法。 2、掌握MATLAB的搜索路径及设置方法。 3、熟悉MATLAB帮助信息的查阅方法 实验内容： 1、建立自己的工作目录，再设置自己的工作目录设置到MA TLAB搜索路径下，再试 验用help命令能否查询到自己的工作目录。 2、在MA TLAB的操作环境下验证课本；例1-1至例1-4，总结MATLAB的特点。 例1-1 例1-2 例1-3 例1-4 3、利用帮助功能查询inv、plot、max、round等函数的功能。 4、完成下列操作： （1）在matlab命令窗口输入以下命令： x=0:pi/10:2*pi; y=sin(x); （2）在工作空间窗口选择变量y，再在工作空间窗口选择回绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量y的图形，并分析图形的含义。 5、访问mathworks公司的主页，查询有关MATLAB的产品信息。 主要教学环节的组织： 教师讲授实验目的、开发环境界面、演示实验过程，然后同学上机练习。 思考题： 1、如何启动与退出MA TLAB集成环境？ 启动： （1）在windows桌面，单击任务栏上的开始按钮，选择‘所有程序’菜单项，然后选择MA TLAB程序组中的MA TLABR2008b程序选项，即可启动 MATLAB系统。 （2）在MA TLAB的安装路径中找到MA TLAB系统启动程序matlab.exe，然后运行它。 （3）在桌面上建立快捷方式后。双击快捷方式图标，启动MA TLAB。 退出： （1）在MA TLAB主窗口file菜单中选择exitMATLAB命令。 （2）在MA TLAB命令窗口中输入exit或quit命令。 （3）单击MATLAB主窗口的关闭按钮。 2、简述MATLAB的主要功能。 MATLAB是一种应用于科学计算领域的数学软件，它主要包括数值计算和符 号计算功能、绘图功能、编程语言功能以及应用工具箱的扩展功能。 3、如果一个MATLAB命令包含的字符很多，需要分成多行输入，该如何处理？ 数学软件课程设计 题目非线性方程求解 班级数学081 姓名曹曼伦 实验目的：用二分法与Newton迭代法求解非线性方程的根； 用Matlab函数solve、fzero、fsolve求解非线性方程（组）的解。 编程实现二分法及Newton迭代法； 学会使用Matlab函数solve、fzero、fsolve求解非线性方程（组）的解。 通过实例分别用二分法及迭代法解非线性方程组并观察收敛速度。 实验内容： 比较求exp(x)+10*x-2的根的计算量。(要求误差不超过十的五次方） (1)在区间（0，1)内用二分法； (2)用迭代法x=(2-exp(x))/10，取初值x=0 。 试验程序 (1)二分法: format long syms x s=exp(x)+10*x-2 a=0; b=1; A=subs(s,a) B=subs(s,b) f=A*B %若f<0，则为由根区间 n=0; stop=1.0e-5; while f<0&abs(a-b)>=stop&n<=100; Xk=(a+b)/2; %二分 M= subs(s, Xk); if M* A<0 symbol=1 %若M= subs(s, Xk)为正，则与a二分 b= Xk else symbol=0 % 若M= subs(s, Xk)为负，则与b二分 a= Xk end n=n+1 end Xk n (2)牛顿迭代法； format long syms x s= (2-exp(x))/10; %迭代公式 f=diff(s); x=0; %迭代初值 a=subs(f,x); %判断收敛性(a是否小于1) s=(2-exp(x))/10; stop=1.0e-5; %迭代的精度 n=0; while a<1&abs(s-x)>=stop&n<=100; x=s %迭代 s=(2-exp(x))/10; n=n+1 end 实验结果： （1）二分法： symbol =1 b =0.50000000000000 n =1 symbol =1 b =0.25000000000000 n =2 symbol =1 b =0.12500000000000 n =3 symbol =0 a =0.06250000000000 n =4 symbol =1 b =0.09375000000000 n =5 symbol =0 a =0.07812500000000 n =6 symbol =1 b =0.09054565429688 n =15 symbol =1 b =0.09053039550781 n =16 symbol =0 a =0.09052276611328 n =17 Xk =0.09052276611328 n =17 （2）迭代法 由x =0.10000000000000 n =1 x =0.08948290819244 n =2 x =0.09063913585958 n =3 x =0.09051261667437 n =4 x =0.09052646805264 n =5 试验结果可见用二分法需要算17次，而用迭代法求得同样精度的解仅用5次，但由于迭代法一般只具有局部收敛性，因此通常不用二分法来求得非线性方程的精确解，而只用它求得根的一个近似解，再用收敛速度较快的迭代法求得其精确解。 Matlab第一次实验报告 2012029010010 尹康 1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。并以x，y为坐标显示图像 x(n+1) = a*x(n)-b*(y(n)-x(n)^2); y(n+1) = b*x(n)+a*(y(n)-x(n)^2) 程序代码： n=input('input the number of pionts:'); a=input('input a:'); b=input('input b:'); x=[]; y=[]; x(1)=input('input x1:'); y(1)=input('input y1:'); %输入点数、初始值以及系数for i=2:n x(i)=a*x(i-1)-b*(y(i-1)-x(i-1)^2); y(i)=a*x(i-1)+b*(y(i-1)-x(i-1)^2); %根据已输入的数据进行迭代end figure;plot(x,y,'linewidth',2) axis equal %横纵坐标等比例 text(x(1),y(1),'1st point') %标记初始点 运行结果： 心得体会及改进：在输入某些数据时，所绘曲线可能是一条折线（如：n=5，a=b=x1=1,y1=2）甚至只有一个点（如：n=5,a=b=x1=y1=1),此时可能出现曲线与坐标轴重合或无法看到点的情况，为了更清晰地展现曲线，可以使线宽适当加宽并标记初始点。 2.编程实现奥运5环图，允许用户输入环的直径。 程序代码： 函数circle： %在指定的圆心坐标处，用指定颜色、宽度的线条绘出指定半径、圆心角的弧 function f=circle(r,x,y,color,linw,alp1,alp2) alp=linspace(alp1,alp2); X=r*cos(alp)+x; Y=r*sin(alp)+y; plot(X,Y,color,'linewidth',linw) end 主程序代码： r=input('input r:'); 西京学院数学软件实验任务书 实验二十一实验报告 一、实验名称：Romberg 积分法，Gauss 型积分法，高斯-勒让德积分法，高斯-切比雪夫积分法，高斯-拉盖尔积分法，高斯-埃尔米特积分法。 二、实验目的：进一步熟悉Romberg 积分法，Gauss 型积分法，高斯-勒让德积分法，高斯-切比雪夫积分法，高斯-拉盖尔积分法，高斯-埃尔米特积分法。 三、实验要求：运用Matlab/C/C++/Java/Maple/Mathematica 等其中一种语言完成程序设计。 四、实验原理： 1．Romberg 积分法： 龙贝格积分法是用里查森外推算法来加快复合梯形求积公式的收敛速度，它的算法如下，其中()i m T 是通过一系列逼近原定积分的龙贝格分值. 计算(0)1[()()]2 b a T f a f b -= + 对1,2,3,k n = ，计算下列各步： 21()(1)1 111 1(21)()[()]222k k k k k j b a j b a T T f a ---=---=++∑ 对1,2,,m k = 和,1,2,,1i k k k =-- ，计算111 441 m i i i m m m m T T T --+-=- 随着计算的步骤的增加，()i m T 越来越逼近积分()b a f x dx ?。 2．Gauss 型积分法： 高斯积分公式的思想是用n 个不等距的节点123,,,n x x x x 对被积函数进行插值，然后对插值后的函数进行积分，其积分公式为： 1 1 1 ()()n k k k f x dx A f x -=≈∑? 如果积分区间不是[1,1]-，则需转换到此区间： 11()()222 b a b a b a b a f x dx f t dt ---+= +? ? 其中系数k A 、节点k x 与n 的关系如下表所示： 3．高斯-切比雪夫积分法： 第一类切比雪夫积分形式为： 1 1 ()()n k k k f x dx A f x -=≈∑? 其中k A n π= ，21cos 2k k x n π-= 4．高斯-拉盖尔积分法： 高斯-拉盖尔公式有两种形式： 1 ()()n x k k k e f x dx A f x +∞ -=≈∑? 信号与系统MATLAB第一次实验报告 一、实验目的 1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。 2.学会运用MATLAB表示常用连续时间信号的方法 3.观察并熟悉一些信号的波形和特性。 4.学会运用MATLAB进行连续信号时移、反折和尺度变换。 5.学会运用MATLAB进行连续时间微分、积分运算。 6.学会运用MATLAB进行连续信号相加、相乘运算。 7.学会运用MATLAB进行连续信号的奇偶分解。 二、实验任务 将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。 三、实验内容 1.MATLAB软件基本运算入门。 1). MATLAB软件的数值计算： 算数运算 向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn为结束值。 矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开； 矩阵的不同行之间必须用分号”;”或者ENTER分开。2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。 举例：计算一个函数并绘制出在对应区间上对应的值。 2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名” 2.MATLAB软件简单二维图形绘制 1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y) 2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p 表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p) 3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin]) 4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’) 5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’) 6).输出：grid on 举例1： 数学实验报告 制作成员班级学号 2011年6月12日 培养容器温度变化率模型 一、实验目的 利用matlab软件估测培养容器温度变化率 二、实验问题 现在大棚技术越来越好，能够将温度控制在一定温度范围内。为利用这种优势，实验室现在需要培植某种适于在8.16℃到10.74℃下能够快速长大的甜菜品种。为达到实验所需温度，又尽可能地节约成本，研究所决定使用如下方式控制培养容器的温度：1，每天加热一次或两次，每次约两小时； 2，当温度降至8.16℃时，加热装置开始工作；当温度达到10.74℃时，加热装置停止工作。 已知实验的时间是冬天，实验室为了其它实验的需要已经将实验室的温度大致稳定在0℃。下表记录的是该培养容器某一天的温度 时间(h)温度（℃）时间（h）温度（℃）09.68 1.849.31 0.929.45 2.959.13 3.878.981 4.989.65 4.988.811 5.909.41 5.908.691 6.839.18 7.008.5217.938.92 7.938.3919.048.66 8.978.2219.968.43 9.89加热装置工作20.848.22 10.93加热装置工作22.02加热装置工作10.9510.8222.96加热装置工作12.0310.5023.8810.59 12.9510.2124.9910.35 13.889.9425.9110.18 三、建立数学模型 1，分析：由物理学中的傅利叶传热定律知温度变化率只取决于温度 差，与温度本身无关。因为培养容器最低温度和最高温度分别是：8.16℃和10.74℃；即最低温度差和最高温度差分别是：8.16℃和10.74℃。而且，16.8/74.10≈1.1467，约为1，故可以忽略温度对温度变化率的影响2， 将温度变化率看成是时间的连续函数，为计算简单，不妨将温度变化率定义成单位时间温度变化的多少，即温度对时间连续变化的绝对值（温度是下降的），得到结果后再乘以一系数即可。 四、问题求解和程序设计流程1）温度变化率的估计方法 根据上表的数据，利用matlab 做出温度-时间散点图如下： 下面计算温度变化率与时间的关系。由图选择将数据分三段，然后对每一段数据做如下处理：设某段数据为{(0x ,0y ),(1x ,1y ),(2x , 2y ),…,(n x ,n y )},相邻数据中点的平均温度变化率采取公式： 温度变化率=（左端点的温度-右端点的温度）/区间长度算得即：v( 2 1i i x x ++)=(1+-i i y y )/(i i x x - +1). 每段首尾点的温度变化率采用下面的公式计算：v(0x )=(30y -41y +2y )/(2x -0x )v(n x )=(3n y -41+n y +2+n y )/(n x -2-n x ) MATLAB程序设计语言 实 验 报 告 专业及班级：电子信息工程 姓名：王伟 学号：1107050322 日期 2013年6月20日 实验一 MATLAB 的基本使用 【一】 实验目的 1.了解MATALB 程序设计语言的基本特点，熟悉MATLAB 软件的运行环境； 2.掌握变量、函数等有关概念，掌握M 文件的创建、保存、打开的方法，初步具备将一般数学问题转化为对应计算机模型处理的能力； 3.掌握二维图形绘制的方法，并能用这些方法实现计算结果的可视化。 【二】 MATLAB 的基础知识 通过本课程的学习，应基本掌握以下的基础知识： 一. MATLAB 简介 二. MATLAB 的启动和退出 三. MATLAB 使用界面简介 四. 帮助信息的获取 五. MATLAB 的数值计算功能 六. 程序流程控制 七. M 文件 八. 函数文件 九. MATLAB 的可视化 【三】上机练习 1. 仔细预习第二部分内容，关于MATLAB 的基础知识。 2. 熟悉MATLAB 环境，将第二部分所有的例子在计算机上练习一遍 3. 已知矩阵???? ??????=??????????=123456789,987654321B A 。求A*B ，A .* B ，比较二者结果是否相同。并利用MATLAB 的内部函数求矩阵A 的大小、元素和、长度以 及最大值。 程序代码： >> A=[1 2 3;4 5 6;7 8 9]; >> B=[9 8 7;6 5 4;3 2 1]; >> A*B ans = 30 24 18 84 69 54 138 114 90 >> A.*B ans = 9 16 21 24 25 24 21 16 9 两者结果不同 >> [m,n]=size(A) m = 3 n = 3 >> b=sum(A) b = 12 15 18 >> a=length(A) a = 3 >>max(A) ans = 数字信号处理及MATLAB 实验报告 班级： 学号： 姓名： 4.7.2 例4,2 设x(n)是由两个正弦信号及白噪声的叠加，试用FFT文件对其作频谱分析。程序清单 %产生两个正弦加白噪声 N=256; f1=.1;f2=.2;fs=1; a1=5;a2=3; w=2*pi/fs; x=a1*sin(w*f1*(0:N-1))+a2*sin(w*f2*(0:N-1))+randn(1,N); %应用FFT求频谱 subplot(2,2,1); plot(x(1:N/4)); title('原始信号'); f=-0.5:1/N:0.5-1/N; x=fft(x); y=ifft(x); subplot(2,2,2); plot(f,fftshift(abs(x))); title('频域信号'); subplot(2,2,3); plot(real(x(1:N/4))); title('时域信号'); 例4.3 设x（n）为长度N=6的矩形序列，用MATLAB程序分析FFT取不同长度时x（n）频谱的变化。N=8,32,64,时x（n）的FFT MATLAB实现程序如下。 x=[1,1,1,1,1,1]; N=8; y1=fft(x,N); n=0:N-1; subplot(3,1,1);stem(n,abs(y1),'.k');axis([0,9,0,6]); N=32; y2=fft(x,N); n=0:N-1; subplot(3,1,2);stem(n,abs(y2),'.k');axis([0,40,0,6]); N=64; y3=fft(x,N); subplot(3,1,3);stem(n,abs(y3),'.k');axis([0,80,0,6]); 数学实验报告 姓名: 班级: 学号: 第一次实验任务 过程: a=1+3i; b=2-i; 结果: a+b =3.0000 + 2.0000i a-b =-1.0000 + 4.0000i a*b = 5.0000 + 5.0000i a/b = -0.2000 + 1.4000i 过程: x=-4.5*pi/180; y=7.6*pi/180; 结果: sin(abs(x)+y)/sqrt(cos(abs(x+y))) =0.2098 心得:对于matlab 中的角度计算应转为弧度。 (1)过程: x=0:0.01:2*pi; y1=sin(x); y2=cos(x); y3=exp(x); y4=log(x); plot(x,y1,x,y2,x,y3,x,y4) plot(x,y1,x,y2,x,y3,x,y4) 结果: (2)过程:>> subplot(2,2,1) >> plot(x,y1) >> subplot(2,2,2) >> plot(x,y2) ./,,,,2,311b a b a b a b a i b i a ?-+-=+=计算、设有两个复数 6,7,5.4) cos()sin(2=-=++y x y x y x ,其中、计算的图形。 下分别绘制)同一页面四个坐标系)同一坐标系下(、在( x y e y x y x y x ln ,,cos ,sin 213==== >> subplot(2,2,3) >> plot(x,y3) >> subplot(2.2.4) >> subplot(2,2,4) >> plot(x,y4) 结果: 心得:在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与hold on进行区别,后者为在同一坐标系中划出多条曲线。 5、随机生成一个3x3矩阵A及3x2矩阵B,计算(1)AB,(2)对B中每个元素平方后得到的矩阵C,(3)sinB,(4)A的行列式,(5)判断A是否可逆,若可逆,计算A的逆矩阵,(6)解矩阵方程AX=B,(7)矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。 过程:A=fix(rand(3,3).*10) ; B=fix(rand(3,3).*10); 《数学实验》报告 实验名称 Matlab 基础知识 学院 专业班级 姓名 学号 2014年 6月 一、【实验目的】 1.认识熟悉Matlab这一软件，并在此基础上学会基本操作。 2.掌握Matlab基本操作和常用命令。 3.了解Matlab常用函数，运算符和表达式。 4.掌握Matlab工作方式和M文件的相关知识。 5.学会Matlab中矩阵和数组的运算。 二、【实验任务】 P16 第4题 编写函数文件，计算 1! n k k = ∑，并求出当k=20时表达式的值。P27第2题 矩阵A= 123 456 789 ?? ?? ?? ?? ?? ，B= 468 556 322 ?? ?? ?? ?? ?? ，计算A*B,A.*B,并比较两者的区别。 P27第3题 已知矩阵A= 52 91 ?? ?? ?? ，B= 12 92 ?? ?? ?? ，做简单的关系运算A>B,A==B,A P27第3题 >> A=[5 2;9 1];B=[1 2;9 2]; >>A>B >>A==B >>A> (A==B)&(A 实验报告 2. The Branching statements 一、实验目的： 1.To grasp the use of the branching statements; 2.To grasp the top-down program design technique. 二、实验内容及要求： 1．实验内容： 1).编写 MATLAB 语句计算 y(t)的值 (Write the MATLAB program required to calculate y(t) from the equation) ???<+≥+-=0 530 53)(2 2t t t t t y 已知 t 从-5到 5 每隔0.5取一次值。运用循环和选择语句进行计算。 (for values of t between -5 and 5 in steps of 0.5. Use loops and branches to perform this calculation.) 2).用向量算法解决练习 1, 比较这两个方案的耗时。 （tic ，toc 的命令可以帮助你完成的时间计算，请使用'help'函数）。 Rewrite the program 1 using vectorization and compare the consuming time of these two programs. (tic, toc commands can help you to finish the time calculation, please use the …help ? function). 2.实验要求： 在报告中要体现top-down design technique, 对于 3 要写出完整的设计过程。 三、设计思路： 1．用循环和选择语句进行计算： 1).定义自变量t ：t=-5:0.5:5; 2).用循环语句实现对自变量的遍历。 3).用选择语句实现对自变量的判断，选择。 4).将选择语句置入循环语句中，则实现在遍历中对数据的选择，从而实现程序的功能。 2. 用向量法实现： 1).定义自变量t ：t=-5:0.5:5; 2).用 b=t>=0 语句，将t>=0得数据选择出，再通过向量运算y(b)=-3*t(b).^2 + 5; 得出结果。 3).用取反运算，选择出剩下的数据，在进行向量运算，得出结果。 四、实验程序和结果 1.实验程序 实验程序：创建m 文件：y_t.m Matlab实验报告 ——定积分的近似计算 学生姓名： 学号： 专业：数学与应用数学专业 数学实验报告 实验序号:1001114030 日期：2012年10月20日 班级应一姓名陈璐学号1001114030 实验名称：定积分的近似运算 问题背景描述： 利用牛顿—莱布尼茨公式虽然可以精确地计算定积分的值，但它仅适合于被积分函数的原函数能用初等函数表达出来的情形。如果这点办不到或不容易办到， 这就有必要考虑近似计算的方法。在定积分的很多应用问题中，被积函数甚至没 有解析表达式，可能只是一条实验记录曲线，或者是一组离散的采样值，这时只 能应用近似方法去计算相应的定积分。 实验目的： 本实验将主要研究定积分的三种近似计算算法：矩形法、梯形法、抛物线发。对于定积分的近似数值计算，Matlab有专门函数可用。 实验原理与数学模型： 1.sum（a）：求数组a的和。 2.format long：长格式，即屏幕显示15位有效数字。 3.double（）：若输入的是字符则转化为相应的ASCII码；若输入的是整型数之则转化为 相应的实型数值。 4.quad（）：抛物线法求数值积分。格式：quad（fun，a，b）。此处的fun是函数，并且 为数值形式，所以使用*、/、^等运算时要在其前加上小数点。 5.trapz():梯形法求数值积分。格式：trapz（x，y）。其中x为带有步长的积分区间;y为数 值形式的运算。 6.fprintf（文件地址，格式，写入的变量）：把数据写入指定文件。 7.syms 变量1变量2……：定义变量为符号。 8.sym（'表达式'）：将表达式定义为符号。 9.int（f，v，a，b）：求f关于v积分，积分区间由a到b。 10.subs（f，'x'，a）：将a的值赋给符号表达式f中的x，并计算出值。若简单地使用subs （f），则将f的所有符号变量用可能的数值代入，并计算出值。 实验所用软件及版本：Matlab 7.0.1 MATLAB实验报告 题目：第一次实验报告 学生姓名： 学院： 专业班级： 学号： 年月 MATLAB第一次实验报告 ————入门第一次上机实验刘老师就MATLAB软件进行了 大致的讲解，并讲了如何建立M文件，定义函数数 组矩阵，如何绘图。先就老师讲解及自己学习的情 况做汇报。 一、建立M文件 <1>M文件建立方法： 1. 在MATLAB中，点:File→New →M-file 2. 在编辑窗口中输入程序内容 3. 点File →Save，存盘，M文件名必须与函数名 一致 <2>课上实例 例：定义函数f(x1,x2)=100(x2-x12)2+(1-x1)2 答：建立M文件：fun.m function f=fun(x) f=100*(x(2)-x(1)^2)^2+(1-x(1))^2 如此便可以直接使用函数fun.m 例如计算f(1,2), 只需在MATLAB命令窗口键入命 令： x=[1 2] fun(x) 得f = 100. <3>课下作业 题目：有一函数，写一程序,输入自变量的值,输出函数值. 解答：建立M文件：zuoye1.m function f=zuoye1(x,y) f=x^2+sin(x*y)+2*y 命令行输入x=1,y=1 zuoye1(x,y) 得ans = 3.8415 经验算答案正确，所以程序正确。 二、定义数组、矩阵 <1>说明 逗号或空格用于分隔某一行的元素，分号用于区分不同的行. 除了分号，在输入矩阵时，按Enter 键也表示开始新一行. 输入矩阵时，严格要求所有行有相同的列 <2>课后作业 题目:有一个4x5矩阵,编程求出其最大值及其所处的位置. 解答：a=round(10*rand (4,5)) [temp I]=max(a) [am II]=max(temp) p=[I(II) II] 运行得一随机矩阵 a = 7 7 7 3 7 0 8 2 0 3 8 7 7 1 10 9 4 0 8 0 temp = Matlab 数学实验报告 一、实验目的 通过以下四组实验,熟悉MATLAB的编程技巧,学会运用MATLAB的一些主要功能、命令，通过建立数学模型解决理论或实际问题。了解诸如分岔、混沌等概念、学会建立Malthu模型和Logistic 模型、懂得最小二乘法、线性规划等基本思想。 二、实验内容 2.1实验题目一 2.1.1实验问题 Feigenbaum曾对超越函数y=λsin(πｘ)（λ为非负实数）进行了分岔与混沌的研究，试进行迭代格式x k+1=λsin(πx k),做出相应的Feigenbaum图 2.1.2程序设计 clear;clf; axis([0,4,0,4]); hold on for r=0:0.3:3.9 x=[0.1]; for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.5) for i=101:150 plot(r,x(i),'k.'); end text(r-0.1,max(x(101:150))+0.05,['\it{r}=',num2str(r)]) end 加密迭代后 clear;clf; axis([0,4,0,4]); hold on for r=0:0.005:3.9 x=[0.1]; for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.1) for i=101:150 plot(r,x(i),'k.'); end end 运行后得到Feigenbaum图 2.2实验题目二 2.2.1实验问题 某农夫有一个半径10米的圆形牛栏，长满了草。他要将一头牛拴在牛栏边界的桩栏上，但只让牛吃到一半草，问拴牛鼻子的绳子应为多长？ 2.2.2问题分析 如图所示，E为圆ABD的圆心，AB为拴牛的绳子，圆ABD为草场，区域ABCD为牛能到达的区域。问题要求区域ABCD等于圆ABC 的一半，可以设BC等于x,只要求出∠a和∠b就能求出所求面积。先计算扇形ABCD的面积，2a÷π×πx2=2aπ2，再求AB的面积，用扇形ABE的面积减去三角形ABE的面积即可。 《数学实验》报告 实验名称 matlab拟合与插值学院机械工程学院 专业班级 姓名 学号 2011年 10月 一、【实验目的】 掌握Matlab关于采用最小二乘法拟合曲线的方法。学会使用matlab求实际中得到数据的插值曲线。 二、【实验任务】 P130第8、10、12题 三、【实验程序】 P130第8题： x=[0.10,0.30,0.40,0.55,0.70,0.80,0.95]; y=[15,18,19,21,22.6,23.8,26]; p1=polyfit(x,y,1); p3=polyfit(x,y,3); p5=polyfit(x,y,5); disp('一阶拟合函数'),f1=poly2str(p1,'x') disp('三阶拟合函数'),f3=poly2str(p3,'x') disp('五阶拟合函数'),f5=poly2str(p5,'x') x1=0.1:0.0017:0.95; y1=polyval(p1,x1); y3=polyval(p3,x1); y5=polyval(p5,x1); plot(x,y,'rp',x1,y1,'--',x1,y3,'k-.',x1,y5); legend('拟合点','一次拟合','三次拟合','七次拟合') P130第10题 x=[10,15,20,25,30]; y=[25.2,29.8,31.2,31.7,29.4]; xi=10:.5:30; yi1=interp1(x,y,xi,'*nearest'); yi2=interp1(x,y,xi,'*linear'); yi3=interp1(x,y,xi,'*spline'); yi4=interp1(x,y,xi,'*cubic'); plot(x,y,'ro',xi,yi1,'--',xi,yi2,'-',xi,yi3,'k.-',xi,yi4,'m:') ,grid on 西安交通大学 高等数学 实验报告 班级 组员与学号 2013年 实验名称：学生成绩管理 一、实验目的 二、实验内容 三、详细编程 clear for i=1:10 a{i}=89+i; b{i}=79+i; c{i}=69+i; d{i}=59+i; end c=[d,c]; Name=input('please input name:'); Score=input('please input score:'); n=length(Score); Rank=cell(1,n); S=struct('Name',Name,'Score',Score,'Rank',Rank); for i=1:n switch S(i).Score case 100 S(i).Rank='满分'; case a S(i).Rank='优秀'; case b S(i).Rank='良好'; case c S(i).Rank='及格'; otherwise S(i).Rank='不及格'; end end disp(['学生姓名 ','得分 ','等级']); for i=1:n disp([S(i).Name,blanks(6),num2str(S(i).Score),blanks(6),S(i).Rank]); end s=0; for i=1:n s=S(i).Score+s; end averscore=s/n; t=S(1).Score; for i=1:(n-1) if(S(i).Score 四川师范大学MATLAB语言实验报告1 系级班年月日 实验名称：Intro, Expressions, Commands 姓名学号指导教师成绩1Objective The objective of this lab is to familiarize you with the MATLAB program development environment and to develop your first programs in this environment. 2Using MATLAB 2.1Starting MATLAB Logon to your computer and start MATLAB by double-clicking on the icon on the desktop or by using the Start Programs menu. MATLAB Desktop window will appear on the screen. The desktop consists of several sub-windows. The most important ones are: ●Command Window (on the right side of the Desktop) is used to do calculations, enter variables and run built-in and your own functions. ●Workspace (on the upper left side) consists of the set of variables (arrays) created during the current MATLAB session and stored in memory. ●Command History (on the lower left side) logs commands entered in the Command Window. You can use this window to view previously run statements, and copy and execute selected statements. You can switch between the Launch Pad window and the Workspace window using the menu tabs under the sub-window on the upper left side. Similarly, you can switch between the Command History and Current Directory windows using the menu tabs under the sub-window on the lower left side. 2.2Executing Commands You can type MATLAB commands at the command prompt “>>” on the Command Window. For example, you can type the formula cos(π/6)2sin(3π/8) as >>(cos(pi/6) ^ 2) * (sin(3 * pi/8)) Try this command. After you finish typing, press enter. The command will be interpreted and the result will be displayed on the Command Window. Try the following by observing how the Workspace window changes: >> a = 2; (M ake note of the usage of “;”) >> b = 3; 数学软件课程设计 题目插值方法与数据拟合 班级数学081 姓名曹曼伦 实验目的：用Matlab语言实现Lagrange插值、Newton基本插值算法，会用Matlab中关于数值微分的函数。 用Matlab语言实现最佳平方逼近及最小二乘算法，会用Matlab中的polyfit、isqnonlin函数进行数据拟和。 实验内容： （1）Lagrange插值： function f=Language(x,y,x0) %求已知数据点的拉格朗日插值多项式 %已知数据点的x坐标向量:x %已知数据点的y坐标向量:y %插值点的x坐标:x0 %求得的拉格朗日插值多项式或在x0处的插值:f x=[0.00.40.81.21.6];%input x data(可替换为自己的数据) y=[00.4283920.7421010.9103140.970348];%input y data(可替换为自己的数据) x0=[0.30.5];%input x0data(可替换为自己的数据) syms t l; if(length(x)==length(y)) n=length(x); else disp('x和y的维数不相等！'); return;%检错 end p=sym(0); for(i=1:n) l=sym(y(i)); for(k=1:i-1) l=l*(t-x(k))/(x(i)-x(k)); end; for(k=i+1:n) l=l*(t-x(k))/(x(i)-x(k)); end; p=p+l; end simplify(p);%简化多项式 f=subs(p,'t',x0);%计算插值点的函数值 f=vpa(f,6);%将插值多项式的值化成6位精度的小数 end （2）Newton基本插值算法 function varargout=newtonliu(varargin) clear,clc x=[0.400.550.650.800.901.05]; fx=[0.410750.578150.696750.888111.026521.25382]; newtonchzh(x,fx); function newtonchzh(x,fx) %由此函数可得差分表 n=length(x); fprintf('*****************差分表*****************************\n'); FF=ones(n,n); FF(:,1)=fx'; for i=2:n for j=i:n FF(j,i)=(FF(j,i-1)-FF(j-1,i-1))/(x(j)-x(j-i+1)); end end for i=1:n fprintf('%4.2f',x(i)); for j=1:i fprintf('%10.5f',FF(i,j)); end fprintf('\n'); end 结果 *****************差分表***************************** 0.400.41075 0.550.57815 1.11600 0.650.69675 1.186000.28000 0.800.88811 1.275730.358930.19733 0.90 1.02652 1.384100.433470.212950.03124 1.05 1.25382 1.515330.524930.228670.031430.00029《MATLAB与数值分析》第一次上机实验报告
MATLAB实验报告50059
Matlab数学实验报告一
matlab第一次实验报告
数学实验“几种常见的求积分近似解的方法”实验报告(内含matlab程序)
MATLAB实验报告(1-4)
matlab数学实验报告5
MATLAB实验报告
MATLAB实验报告
浅析Matlab数学实验报告
MATLAB全实验报告
matlab实验报告
matlab实验报告
MATLAB入门实验报告
MATLAB数学实验报告
《数学实验》报告matlab-第五次作业
高等数学实验报告matlab
MATLAB实验报告(8个实验)
Matlab数学实验报告三