实验任务
XXXXXX （1）阅读例子程序，观察输出波形，理解每条语句的含义。
（2）已知有限长序列x(n)=[7,6,5,4,3,2],求DFT和IDFT，要求：画出序列傅立叶变换对应的幅度谱和相位谱；画出原信号与傅立叶逆变换IDFT[X(k)]的图形进行比较。（3）已知周期序列的主值x(n)=[7,6,5,4,3,2],求x(n)周期重复次数为3次时的DFS和IDFS。要求：画出原信号序列的主值和周期序列的图形；画出离散傅立叶变换对应的幅度谱和相位谱。
（4）求x(n)=[7,6,5,4,3,2], 0= %% %题（2） close all;clear all;clc; xn=[7,6,5,4,3,2]; N=length(xn); n=0:N-1;k=0:N-1; Xk=xn*exp(-j*2*pi/N).^(n'*k); x=(Xk*exp(j*2*pi/N).^(n'*k))/N; subplot(2,2,1),stem(n,xn,'k'); title('x(n)');axis([-1,N,0.8*min(xn),1.1*max(xn)]);grid; subplot(2,2,2),stem(n,x,'k'); title('IDFT|X(k)|'); axis([-1,N,0.8*min(x),1.1*max(x)]);grid; subplot(2,2,3),stem(k,abs(Xk),'k'); title('|X(k)|'); axis([-1,N,min(abs(Xk)),1.1*max(abs(Xk))]); subplot(2,2,4),stem(k,angle(Xk),'k'); title('arg|X(k)|'); axis([-1,N,min(angle(Xk)),1.1*max(angle(Xk))]); %% %题（3） close all;clear all;clc; xn=[7,6,5,4,3,2]; N=length(xn); n=0:3*N-1;k=0:3*N-1; xn1=xn(mod(n,N)+1); Xk=xn1*exp(-j*2*pi/N).^(n'*k); subplot(2,2,1),stem(xn,'k'); title('原主值信号x(n)'); subplot(2,2,2),stem(n,xn1,'k'); title('周期序列信号'); axis([-1,4*N,1.1*min(xn1),1.1*max(xn1)]); subplot(2,2,3),stem(k,abs(Xk),'k'); title('|X(k)|'); axis([-1,4*N,1.1*min(abs(Xk)),1.1*max(abs(Xk))]); subplot(2,2,4),stem(k,angle(Xk),'k'); title('arg|X(k)|'); axis([-1,4*N,1.1*min(angle(Xk)),1.1*max(angle(Xk))]); %% %题（4） clear all;clc xn=[7,6,5,4,3,2]; N=length(xn); n=0:N-1; w=linspace(-2*pi,2*pi,50); X=xn*exp(-j*n'*w); subplot(221),stem(n,xn,'k'); ylabel('x(n)'); subplot(223),stem(w,abs(X),'k','.'); axis([-2*pi,2*pi,1.1*min(abs(X)),1.1*max(abs(X))]); ylabel('频度谱'); subplot(224),stem(w,angle(X),'k','.'); axis([-2*pi,2*pi,1.1*min(angle(X)),1.1*max(angle(X))]); ylabel('相位谱'); %% clear all;clc N=100; xn=[7,6,5,4,3,2,zeros(1,N-6)]; n=0:N-1;k=0:N-1; Xk=xn*exp(-j*2*pi/N).^(n'*k); x=(Xk*exp(j*2*pi/N).^(n'*k))/N; subplot(2,2,1),stem(n,xn,'k'); title('x(n)');axis([0,N,1.1*min(xn),1.1*max(xn)]); subplot(2,2,2),stem(n,abs(x),'k'); title('IDFT|X(k)|'); axis([0,N,1.1*min(x),1.1*max(x)]); subplot(2,2,3),stem(k,abs(Xk),'k'); title('|X(k)|'); axis([0,N,1.1*min(abs(Xk)),1.1*max(abs(Xk))]); subplot(2,2,4),stem(k,angle(Xk),'k'); title('arg[X(k)]'); axis([0,N,1.1*min(angle(Xk)),1.1*max(angle(Xk))]); 题（2） 题（3） 题（4） 《Matlab 编程训练》 作业 专 业 学生姓名 班级 学 号 指导教师 完成日期 实训一 MATLAB 语言介绍和数值计算 1．先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存变量。 12 2sin851z e =+ . 2. 已知 1234413134787,2033657327A B --???? ????==???? ????-???? ，求下列表达式的值： (1) A+6*B 和A-B+I （其中I 为单位矩阵） A+6*B: A-B+I: (2)A*B和A.*B A*B程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] c=A*B 结果： A.*B程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] D=A.*B 结果： (3)A^3和A.^3 A^3程序： A=[12 34 -4;34 7 87;3 65 7] E=A^3 结果： A．^3程序： A=[12 34 -4;34 7 87;3 65 7] C=A.^3 (4)A/B及B\A A/B程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] C=A/B 结果： B\A程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] D=B\A 结果： (5)将矩阵C=B\A的右下角2*2子矩阵赋给D, 并（3）保存变量（mat文件）程序： A=[12 34 -4;34 7 87;3 65 7]; B=[1 3 -1;2 0 3;3 -2 7]; C=B*inv(A); D=C(2:3,2:3) 结果： 实验一 1、熟悉MATLAB的窗口结构（命令窗口、历史命令窗口、工作区窗口、当前目录 窗口） 2、掌握命令窗口中基本命令的使用 3、在命令窗口中，给定圆的半径r，求得圆的周长c和面积s，并查看工作区窗 口的变化 4、将r,c,s变量保存到磁盘文件abc.mat中，并删除内存变量r,c,s，查看工作 区窗口的变化 5、将abc.mat文件中变量装入内存，查看工作区窗口的变化 6、将历史命令窗口中的命令再装入命令窗口中使用 7、改变当前目录，查看当前目录窗口的变化 8、掌握命令窗口中 cd,quit,help,date,dir,ls,what,who,clocl,fix(clock),format,save,loa d,clc,clear等命令的使用 9、注意各种MATLAB版本的差别 实验二 1、在命令窗口中，输入长方形的长和宽，求长方形的周长和面积 2、输入三角形的三条边（要满足构成三角形的条件），求三角形的周长和面积 3、掌握MATLAB中各标准函数的使用（sin,cos,sind,fix,mod,…） 4、用fprintf输出各种类型的数据（如fprintf('a=%d\n',123) a=123 >> fprintf('b=%f\n',123.456) b=123.456000 >> fprintf('c=%c\n','A') c=A……) 实验三 1、在编辑窗口中：输入学生成绩，输出该成绩的等级。等级规定如下：[90， 100]为A等，[80，90)为B等，[70，80)为C等，[60，70)为D等，[0，60)为E等。要求用if和 switch两种方法实现。 2、商场购物，100件以下，不优惠，100~199件95折，200~399件90折，400~799 件85折，800~1499件80折，1500件以上，75折。输入所购货物的单价、件数，求实际付款数目。要求用if和 switch两种方法实现（在编辑窗口中实现）。 实验四 1、求两个正整数的最大公约数和最小公倍数（在编辑窗口中实现，命令窗口中 调用）。 2、求100~300内所有素数（在编辑窗口中实现，命令窗口中调用）。 第一章 M A T L A B 概况与基本操作 1.选择题（每题2分，共20分）： (1)最初的MATLAB 核心程序是采用D 语言编写的。 (2)即将于2011年9月发布的MATLAB 新版本的编号为C 。 2011Ra 2011Rb R2011b R2011a (3)在默认设置中，MATLAB 中的注释语句显示的颜色是B 。 A.黑色 B.绿色 C.红色 D.蓝色 (4)如果要以科学计数法显示15位有效数字，使用的命令是A 。 long e long long g long d (5)在命令窗口新建变量a 、b ，如果只查看变量a 的详细信息，使用的命令为A 。 a a (6)如果要清除工作空间的所有变量，使用的命令为C 。 all C.两者都可 D.两者都不可 (7)在创建变量时，如果不想立即在命令窗口中输出结果，可以在命令后加上B 。 A.冒号 B.分号 C.空格 D.逗号 (8)如果要重新执行以前输入的命令，可以使用D 键。 A.下箭头↓ B.右箭头→ C.左箭头← D.上箭头↑ (9)如果要查询函数det 的功能和用法，并显示在命令窗口，应使用命令C 。 D.三者均可 (10)如果要启动Notebook 文档，下列D 操作是可行的。 A.在命令窗口输入notebook 命令 B.在命令窗口输入notebook filename 命令 C.在Word 中启动M-book 文档 D.三者均可 2.填空题（每空1分，共20分）： (1)MATLAB 是matrix 和laboratory 两个单词前三个字母的组合，意为“矩阵实验室”，它的创始人是Cleve Moler 和Jack Little 。 (2)在MATLAB 的默认设置中，关键字显示的字体为蓝色，命令、表达式、计算结果显示的字体为黑色，字符串显示的字体为褐红色，注释显示的字体为绿色，错误信息显示的字体为红色。 (3)在命令窗口中，输出结果显示为各行之间添加空行的命令为format loose ，各行之间不添加空行的命令为format compact 。 (4)在MATLAB 中，各种标点符号的作用是不同的。例如，空格的作用是分隔数组每行各个元素，逗号的作用是分隔数组每行各个元素或函数的各个输入参数，分号的作用是作为不显示命令结果的命令行的结尾或分隔数组各列，冒号的作用是生成一维数组或表示数组全部元素，百分号的作用是引导一行注释，…的作用是连接相邻两行，感叹号的作用是调用操作系统命令。 3.程序设计题（每题10分，共40分） (1)以25m/s 的初速度向正上方投球（g=s 2 ），计算到达最高点的时间tp 以及球从出发点到最高点的距离hp 。 解：根据物理学知识，物体上抛运动的速度与经过的时间之间的关系为0p p v v gt =-，因此所需要的时间为0p p v v t g -= 。而到达最高点时的速度0p v =，因此可根据此公式求出tp ： v0=25;g=;vp=0; tp=(v0-vp)/g tp = 选 题 说 明 本人选做第2、4、5、9、11、12、13、14、16、19、24 题。 作业内容题目2：问题描述：在[0 ， 2π]范围内绘制二维曲线图y=cos(5x)*sin(x) （1）问题分析 这是一个二维绘图问题，先写出x的取值范围，再用plot函数画出y的图像。 （2）软件说明及源代码 >> x = 0:pi/100:2.*pi; y=cos(5*x).*sin(2*x); >> plot(x,y) （3）实验结果 题目4：问题描述：创建符号函数并求解，要求写出步骤和运行结果 (1)创建符号函数f=ax2+bx+c (2)求f=0的解 （1）问题分析 这是符号计算问题，首先要确定符号变量，然后创建符号函数，最后利用subs函数求解特值。 （2）软件说明及源代码 >> syms a b c x f; f=a*x^2+b*x+c; subs(f,0) （3）实验结果 ans = c 题目5：问题描述：求积分 （1）问题分析 这是符号计算的积分求解问题，首先需要确定符号变量，然后利用int函数计算积分。 （2）软件说明及源代码 >> syms x y; y=sqrt(1-2*sin(2*x)); >> int(y,x,0,pi/2) （3）实验结果 ans = ellipticE(-pi/4, 4)*1i - ellipticE(pi/4, 4)*1i - ellipticE(-pi/6, 4)*2i + ellipticE(pi/6, 4)*2i 题目9：问题描述：按水平和竖直方向分别合并下述两个矩阵： （1）问题分析 这是考查矩阵的基本操作，首先定义矩阵，然后合并矩阵。 （2）软件说明及源代码 >> A=[1,0,0;1,1,0;0,0,1]; B=[2,3,4;5,6,7;8,9,10]; >> a=[A,B],b=[A;B] （3）实验结果 a = 1 0 0 2 3 4 1 1 0 5 6 7 0 0 1 8 9 10 b = 1 0 0 1 1 0 0 0 1 2 3 4 5 6 7 8 9 10 题目11：问题描述：计算z＝yx2+3y2x+2y3的和： （1）问题分析 这是符号计算问题，首先确定符号变量，然后构造函数，最后利用diff函数进行求导。 （2）软件说明及源代码 >> syms x y z; >> z=y*x^2+3*y^2*x+2*y^3; >> diff(z,y,1),diff(diff(z,y,1),x,1) （3）实验结果 ans = x^2 + 6*x*y + 6*y^2 Matlab程序设计 班级 姓名 学号 《MATLAB程序设计》作业 1、考虑如下x-y 一组实验数据： x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot的原始数据、一次拟合曲线和三次拟合曲线，给出MATLAB代码和运行结果。 代码如下： x=[1,2,3,4,5,6,7,8,9,10]; y=[1.2,3,4,4,5,4.7,5,5.2,6,7.2]; plot(x,y); title('原始数据'); p=polyfit(x,y,1); q=polyval(p,x); figure,plot(x,q); title('一次拟合'); p=polyfit(x,y,2); q=polyval(p,x); figure,plot(x,q); title('二次拟合'); 运行结果如下： 1 2 3 4 5 6 7 8 9 10 12 3 4 5 6 7 8 原始数据 123 456789 102 2.5 3 3.54 4.5 55.56 6.57一次拟合 123456789 101 2 3 4 5 6 7 二次拟合 2、在[0，3π]区间，绘制y=sin(x)曲线（要求消去负半波，即(π，2π)区间内的函数值置零），求出曲线y 的平均值，以及y 的最大值及其最大值的位置。给出执行代码和运行结果。 代码如下： clear clc x=(0:0.01:3*pi); y=sin(x); plot(x,y); y1=(y>=0).*y; figure,plot(x,y1); 1．求解方程组sin(x-y)=0，cos(x-y)=0写出求解的MATLAB程序？（写出求解的MATLAB的程序，采用solve命令） 答案 >> [x,y]=solve('sin(x-y)=0','cos(x+y)=0','x','y') x = -1/4*pi 1/4*pi y = -1/4*pi 1/4*pi 2．设3x’’-5x’+6x=y,x(0)=3,x’(0)=0,求y（5）的值？（写出求解该微分方程的MATLAB程序，用dsolve命令） %ex1_9.m f='3*D2x-5*x+6*x-y=0';x=dsolve(f,'Dx(0)=0,x(0)=0','y') %ex1_10.m y=solve('-sin(1/3*3^(1/2)*y)*3^(1/2)+y-x=0','x=5','x','y') x = -sin(1/3*3^(1/2)*y)*3^(1/2)+y y = x: [1x1 sym] y: [1x1 sym] 3. 1 3 5 7 A= 9 2 4 6 求B=A^T,C=A^-1，用for循环 8 10 3 5 求解A各元素的和，以及A中的最小元 答案 >> A=[1 3 5 7;9 2 4 6;8 10 3 5],B=A' A = 1 3 5 7 9 2 4 6 8 10 3 5 B = 1 9 8 3 2 10 5 4 3 7 6 5 %ex1_5.m s=0 m=A(1,1) for i=1:3 for j=1:4 if ( A(i,j)<=m ) ; m=A(i,j) end s=s+A(i,j) end end m s m=1 s=63 实验一 MATLAB 语言的控制语句 实验学时： 2 学时 实验目的： 1、了解MATLAB 语言的体系结构及MATLAB 语言的特点； 2、掌握MATLAB 循环语句、条件转移语句等的使用方法，能编写相应程序实现所需功能。 实验内容： 一、使用 for 循环、while 循环计算6302i i k == ∑ 1、for 循环 确定循环变量的起始值、终止值； 编写循环语句组，实现所需功能； 运行程序，并记录结果。 2、while 循环 确定循环的条件表达式； 编写循环语句组，实现所需功能； 运行程序，并记录结果。 二、采用搜寻法寻找矩阵23121416524332372148656166 646463A ??????=?????? 中最小的元素，以及它的位置。 实验报告： 一、【实验目的】 了解matlab基础知识和操作 一、【实验任务】 P16 4. 编写函数文件，计算∑k!，并求出当n=20时表达式的值。 P27 2. 矩阵A=[1,2,3;4,5,6;7,8,9]，B=[4,6,8;5,5,6;3,2,2]，计算A*B，A.*B，并比较两者的区 别。 P27 3. 已知矩阵A=[5,2;9,1]，B=[1,2;9,2]，做简单的关系运算A>B，A= =B，A ans = 4 12 24 2 25 36 21 16 18 P27 第3题 >> A=[5 2;9 1]; >> B=[1 2;9 2]; >> A>B ans = 1 0 0 0 >> A==B ans = 0 1 1 0 >> A> (A==B)&(A MATLAB大作业 作业要求： （1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。 （2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果分析、本课程学习体会、列出主要的参考文献。打印稿不要求MATLAB程序，但电子稿要包含MATLAB 程序。 （3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话则长，无话则短。不要抄袭复制，可以参考网上、文献资料的内容，但要理解，要变成自己的语言，按自己的思路组织内容。 （4）从给出的问题中至少选择一题(多做不限，但必须独立完成，严禁抄袭)。 （5）大作业占过程考核的20%，从完成情况、工作量、作业文档方面评分。 第一类：绘制图形。（B级） 问题一：斐波那契（Fibonacci）螺旋线，也称黄金螺旋线（Golden spiral），是根据 斐波那契数列画出来的螺旋曲线，自然界中存在许多斐波那契螺旋线的图案，是自然界最完美的经典黄金比例。斐波那契螺旋线，以斐波那契数为边的正方形拼成的长方形，然后在正方形里面画一个90度的扇形，连起来的弧线就是斐波那契螺旋线，如图所示。 问题二：绘制谢尔宾斯基三角形（Sierpinskitriangle）是一种分形，由波兰数学家谢 尔宾斯基在1915年提出,它是一种典型的自相似集。其生成过程为：取一个实心的三角形（通常使用等边三角形），沿三边中点的连线，将它分成四个小三角形，然后去掉中间的那一个小三角形。接下来对其余三个小三角形重复上述操作，如图所示。 问题三：其他分形曲线或图形。分形曲线还有很多，教材介绍了科赫曲线，其他还有皮 亚诺曲线、分形树、康托（G. Cantor）三分集、Julia集、曼德布罗集合（Mandelbrot set），等等。这方面的资料很多（如https://www.doczj.com/doc/1811215325.html,/content/16/0103/14/5315_525141100.shtml），请分析构图原理并用MATLAB实现。 问题四：模拟掷骰子游戏：掷1000次骰子，统计骰子各个点出现的次数，将结果以下表的形式显示，并绘制出直方图。 点数 1 2 3 4 5 6 出现次数166 150 164 162 184 174 问题五：利用MATLAB软件绘制一朵鲜花，实现一定的仿真效果。 提示：二维/三维绘图，对花瓣、花蕊、叶片、花杆等的形状和颜色进行详细设置。 第二类：插值与拟合。（B级） 问题一：有人对汽车进行了一次实验，具体过程是，在行驶过程中先加速，然后再保持匀速行驶一段时间，接着再加速，然后再保持匀速，如此交替。注意，整个实验过程中从未 （1）分别使用最近点插值、线性插值、三次埃尔米特插值和三次样条插值进行计算[0,110]时间段50个时间点的速度。 （2）绘制插值图形并标注样本点。 问题二：估算矩形平板各个位置的温度。已知平板长为5m，宽为3m，平板上3×5栅格 点上的温度值为44，25，20，24，30；42，21，20，23，38；25，23，19，27，40。 （1）分别使用最近点插值、线性插值和三次样条插值进行计算。 （2）用杆图标注样本点。 （3）绘制平板温度分布图。 问题三：自行车道的设计。对9条道路上的自行车道宽度以及自行车与过往机动车之间 （1）对数据进行线性拟合。 （2）绘制拟合曲线和样本点。 （3）如果自行车与过往机动车之间安全距离的最小距离是1.8m，试计算相应的自行车道宽度的最小值。 问题四：在水资源工程学中，水库的大小与为了蓄水而拦截的河道中的水流速度密切相 关。对于某些河流来说，这种长时间的历史水流记录很难获得。然而通常容易得到过去若干年间关于降水量的气象资料。鉴于此，推导出流速与降水量之间的关系式往往特别有用。只 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:'); Matlab 大作业 （组内成员：彭超杰、南彦东、江明伟） 一、研究模型 （电车）通过控制油门（保持一定角度）来调节电动机能输出稳定的转速，从而控制车速稳定。 数学依据说明如下： 由图可知存在以下关系：a d a a u w k R i dt di L =++ (w k e d d =) L M M dt dw J -= a m i k M = L a m M i k dt dw J -= k为反电势常数，m k为电动机电磁力矩常数，这里忽略阻尼力矩。d 二、数学模型 再看整个研究对象，示意图以课本为依据，不同点是这里将数控的进给运动，转换为汽车行驶所需要的扭矩。（这里不说明扭矩的具体产生过程，仅仅说明输出车轮旋转的角速度w ） 对照课本不同，() s θ变为()s N ，1 221z z w w =，1w 为电动机的转速，2w 为轮胎的转速，1z 为电动机的光轴齿轮的齿数，2z 为与轮胎相连光轴的 齿轮齿数。 )(*10110w x w k x ==，1 21z z k = ()c a m m d b a m x K K K k s k k JRs JLs K K K k s G i 1231+++= () c a m m d M K K K k s k k JRs JLs R Ls K s G L 1231)(++++-= 同理，忽略电枢绕组的电感L ，简化系统传递函数方框图如下 ()JR K K K k JR s k k s JR K K K k s G c a m m d b a m x i 121++= ()JR K K K k JR s k k s K K K K k s k k Rs R K s G c a m m d c a m m d M L 121121++-=++-= 实验一Matlab使用方法和程序设计 一、实验目的 1、掌握Matlab软件使用的基本方法； 2、熟悉Matlab的数据表示、基本运算和程序控制语句 3、熟悉Matlab绘图命令及基本绘图控制 4、熟悉Matlab程序设计的基本方法 二、实验内容： 1、帮助命令 使用help命令，查找sqrt（开方）函数的使用方法； 解：sqrt Square root Syntax B = sqrt(X) Description B = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results. Remarks See sqrtm for the matrix square root. Examples sqrt((-2:2)') ans = 0 + 1.4142i 0 + 1.0000i 1.0000 1.4142 2、矩阵运算 （1）矩阵的乘法 已知A=[1 2;3 4]; B=[5 5;7 8]; 求A^2*B 解：A=[1 2;3 4 ]; B=[5 5;7 8 ]; A^2*B （2）矩阵除法 已知A=[1 2 3;4 5 6;7 8 9]; B=[1 0 0;0 2 0;0 0 3]; A\B,A/B 解：A=[1 2 3;4 5 6;7 8 9 ]; B=[1 0 0;0 2 0;0 0 3 ]; A\B,A/B （3）矩阵的转置及共轭转置 《数学实验与MATLAB》 ——综合实验报告 实验名称：不同温度下PDLC薄膜的通透性 与驱动电压的具体关系式的研究学院：计算机与通信工程学院 专业班级： 姓名： 学号： 同组同学： 2014年 6月10日 一、问题引入 聚合物分散液晶(PDLC)是将低分子液晶与预聚物Kuer UV65胶相混合，在一定条件下经聚合反应，形成微米级的液晶微滴均匀地分散在高分子网络中，再利用液晶分子的介电各向异性获得具有电光响应特性的材料，它主要工作在散射态和透明态之间并具有一定的灰度。聚合物分散液晶膜是将液晶和聚合物结合得到的一种综合性能优异的膜材料。该膜材料能够通过驱动电压来控制其通透性，可以用来制作PDLC型液晶显示器等，具有较大的应用范围。已知PDLC薄膜在相同光强度及驱动电压下，不用的温度对应于不同的通透性，不同温度下的阀值电压也不相同。为了尽量得到不同通透性的PDLC薄膜，有必要进行温度对PDLC薄膜的特性的影响的研究。现有不同温度下PDLC 薄膜透过率与驱动电压的一系列数据，试得出不同温度下PDLC薄膜通透性与驱动电压的具体关系式，使得可以迅速得出在不同温度下一定通透性对应的驱动电压。 二、问题分析 想要得到不同温度下PDLC薄膜通透性与驱动电压的具体关系式可以运用MATLAB多项式农合找出最佳函数式，而运用MATLAB多项式插值可以得出在不同温度下一定通透性所对应的驱动电压。 三、实验数据 选择10、20、30摄氏度三个不同温度，其他条件一致。 （1）、10摄氏度 实验程序： x=2:2:40; y=[5.2,5.4,5.8,6.4,7.2,8.2,9.4,10.8,12.2,14.0,16.6,22.0, 30.4,39.8,51.3,55.0,57.5,58.8,59.6,60.2]; p3=polyfit(x,y,3); p5=polyfit(x,y,5); p7=polyfit(x,y,7); disp('三次拟合函数'),f3=poly2str(p3,'x') disp('五次拟合函数'),f5=poly2str(p5,'x') disp('七次拟合函数'),f7=poly2str(p7,'x') x1=0:1:40; y3=polyval(p3,x1); y5=polyval(p5,x1); y7=polyval(p7,x1); plot(x,y,'rp',x1,y3,'--',x1,y5,'k-.',x1,y7); legend('拟合点','三次拟合','五次拟合','七次拟合') 实验结果： 自控控制原理 MATLAB 程序设计作业 指导老师:汪晓宁 目录 一、题目 (2) 二、运行结果 (3) 三、程序说明 (8) 四、附录 ............................................ 9 代码 . ............................................. 9 参考文献 .. (17) 一、题目 用 Matlab 创建用户界面,并完成以下功能 a 将产生未综合系统的根轨迹图以及 0.707阻尼比线, 你可以交互地选择交点的运行点。界面能显示运行点的坐标、增益值以及近似为二阶系统估算的超调量、调整时间、峰值时间、阻尼比、无阻尼自然震荡频率以及稳态误差 b 显示未综合系统的阶跃响应 c 输入控制器的参数, 绘制综合后系统的根轨迹图以及显示综合的设计点 (主导极点 , 允许不断改变控制器参数,知道所绘制的根轨迹通过设计点 d 对于综合后的系统, 显示运行点的坐标、增益,近似为二阶系统估算的超调量、调整时间、峰值时间、阻尼比、无阻尼自然震荡频率以及误差系数 e 显示综合后系统的阶跃响应 二、运行结果 输入传递函数分子分母 生成根轨迹图 选择点并得到该点各项参数在下方输出面板输出 获得阶跃响应图 用 rltool(辅助,选择合适的插入零点 输入零点,并得到根轨迹图 选择根轨迹图上的任一点,得到数据,在下方输出面板输出得到阶跃响应图 三、运行说明 第一步, 在请输入分子后的输入框输入传递函数分子的矩阵, 在下一输入框输入传递函数分母并按“生成根轨迹图”按钮获得根轨迹 第二步, 按选择点并显示各参数获得根轨迹图上任一点的各项数据, 数据全部输出在下方输出面板 第三步,按“生成阶跃响应图”按钮可以获得该函数的阶跃响应 第四步,在“请输入插入零点”后的输入框中输入参数,并按“生成综合后根轨迹图” 按钮产生根轨迹 (可以通过点击“根轨迹校正”按钮,调用工具箱拖动零点进行快速查看根轨迹图,选择合适的根轨迹再在输入框中输入零点的值 第五步,按“选择点并显示各参数(综合后系统”选取各点,查阅参数,数据输出在下方输出面板上 第六步,按“生成阶跃响应图(综合后系统”可以得到综合后系统的阶跃响应 最后,点击“退出”结束程序 四、附录 代码: function varargout = Liushuai20122510(varargin % LIUSHUAI20122510 MATLAB code for Liushuai20122510.fig % LIUSHUAI20122510, by itself, creates a new LIUSHUAI20122510 or raises the existing % singleton*. % 一、必答题： 1. MATLAB系统由那些部分组成？ 答：MATLAB系统主要由开发环境、MATLAB语言、MATLAB数学函数库、图形功能和应用程序接口五个部分组成。 2. 如何启动M文件编辑/调试器？ 答：在操作界面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入“edit”命令也可以启动M文件编辑/调试器。 3. 存储在工作空间中的数组能编辑吗？如何操作？ 答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 4. 在MATLAB中有几种获得帮助的途径？ 答：在MATLAB中有多种获得帮助的途径： （1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器； （2）help命令：在命令窗口键入“help” 命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息； （3）lookfor命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数 （4）模糊查询：输入命令的前几个字母，然后按Tab键，就可以列出所有以这几个字母开始的命令和函数。 5. 有几种建立矩阵的方法？各有什么优点？ 答：（1）以直接列出元素的形式输入； （2）通过语句和函数产生； （3）.在m文件中创建矩阵； （4）从外部的数据文件中装入。 6. 命令文件与函数文件的主要区别是什么？ 答：命令文件: M文件中最简单的一种，不需输出输入参数，用M 文件可以控制工作空间的所有数据。运行过程中产生的变量都是全局变量。运行一个命令文件等价于从命令窗口中顺序运行文件里的命令，程序不需要预先定义，只要依次将命令编辑在命令文件中即可。 函数文件：如果M文件的第一个可执行行以function开始，便是函数文件，每一个函数文件定义一个函数。函数文件区别于命令文件之处在于命令文件的变量在文件执行完成后保留在工作空间中，而函数文件内定义的变量只在函数文件内起作用，文件执行完后即被清除。 7. 创建符号变量有几种方法？ M a t l a b程序设计大作 业 公司内部档案编码：[OPPTR-OPPT28-OPPTL98-OPPNN08] Matlab程序设计 课程大作业 题目名称：_________________________________ 班级：_________________________________ 姓名：_________________________________ 学号：_________________________________ 课程教师：温海骏 学期： 2015-2016学年第2学期 完成时间： MATLAB优化应用 §1 线性规划模型 一、线性规划问题： 问题1：生产计划问题 假设某厂计划生产甲、乙两种产品，现库存主要材料有A类3600公斤，B类2000公斤，C类3000公斤。每件甲产品需用材料A类9公斤，B类4公斤，C类3公斤。每件乙产品，需用材料A类4公斤，B类5公斤，C类10公斤。甲单位产品的利润70元，乙单位产品的利润120元。问如何安排生产，才能使该厂所获的利润最大。 问题2：投资问题 某公司有一批资金用于4个工程项目的投资，其投资各项目时所得的净收益(投入资金百分比)如下表：工程项目收益表 由于某种原因，决定用于项目A的投资不大于其他各项投资之和而用于项目B和C的投资要大于项目D的投资。试确定该公司收益最大的投资分配方案。 问题3：运输问题 有A 、B 、C 三个食品加工厂，负责供给甲、乙、丙、丁四个市场。三个厂每天生产食品箱数上限如下表： 四个市场每天的需求量如下表： 从各厂运到各市场的运输费(元/每箱)由下表给出： 求在基本满足供需平衡的约束条件下使总运输费用最小。 §2 多目标规划模型 多目标规划定义为在一组约束下，多个不同的目标函数进行优化设计。 数学模型： 12min ()() ().()0,1,2, ,m j f x f x f x st g x j k ???? ≤= 其中x=(x 1 ,x 2 , … ,x n )为一个n 维向量；f i (x)为目标函数，i=1, 2, … ,m; g j (x)为系统约束, j=1, 2, … ,k 。 《MATLAB》期末大作业 学院土木工程与建筑学院 专业 班级 姓名 指导教师李琳 2018 年 5 月16 日 明 作业内容题目2:问题描述：在[0 2π]范围内绘制二维曲线图y=cos(5x)*sin(x) （1）问题分析 这是一个二维绘图问题，先划定x的范围与间距，再列出y的表达式，利用plot函数绘制二维曲线。 （2）软件说明及源代码 >> x = 0:pi/10:2*pi; >>y = cos(5*x).*sin(x); >>plot(x,y) （3）实验结果 题目4:问题描述：创建符号函数并求解，要求写出步骤和运行结果 (1)创建符号函数f=ax2+bx+c (2)求f=0的解 （1）问题分析 这是一个符号函数显示以及符号函数的求解问题，第一问先定义常量与变量，在写出f表达式，利用pretty函数显示f。第二问利用solve函数求解f=0时的解。 （2）软件说明及源代码 第一问 >> syms a b c x; >> f=a*x^2+b*x+c; >> pretty(f) 第二问 >>syms a b c x; >>f=a*x^2+b*x+c; >> solve(f) （3）实验结果 1、 2、 题目5:问题描述：求积分 （1）问题分析 这是一个利用符号函数求积分的问题，先定义变量x，再列出I1表达式，利用int函数求在范围0到Pi/2上的积分。 （2）软件说明及源代码 >> syms x; >> I1=(1-2*sin(2*x))^0.5; >> int(I1,0,0.5*pi) （3）实验结果 题目6:问题描述：分别随机产生一个6×6的整数矩阵（元素可在［－20,20］之间），求该随机阵的秩，特征值和特征向量。 （1）问题分析 这是一个矩阵运算问题，先利用rand函数产生一个6*6的元素在-20到20 Matlab程序设计作业 姓名: 学号: 专业: ? MATLAB 程序设计》作业 1、考虑如下x-y 一组实验数据: x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y 二[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot 的原始数据、一次拟合曲线和三次拟合曲线，给出 原始曲线 MATLAB 代码和运行结果。 7 6 5 4 3 2 2 3 4 5 6 7 8 9 10 7 6.5 6 5.5 5 4.5 4 3.5 3 2.5 10 一次拟合 三次拟合 x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]; y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2]; figure; plot(x,y) p1=polyfit(x,y,1); y1=polyval(p1,x); figure; plot(x,y1) p2=polyfit(x,y,3); y2=polyval(p2,x); figure; plot(x,y2) 2、在［0, 3n区间，绘制y二Sin(x)曲线(要求消去负半波，即(n 2n)区间内的函数值置零)，求出曲线y 的平均值，以及y 的最大值及其最大值的位置。给出执行代码和运行结果。 x=0:pi/1000:3*pi; y=Sin(x); y1=(y>=0).*y; %消去负半波figure(1); plot(x,y1, 'b' ); a=mean(y1) %求出y1 的平均值 b=max(y1) %求出y1 的最大值b, 以及最大值在矩阵中的位置; d=x(find(y1==b)) >> ex1 a = 0.4243 b = 1 d = 1.5708 7.8540 >> Matlab在汽车振动分析中的应用 XXX （昆明理工大学交通工程学院，昆明650500） 摘要：在明确汽车振动的产生原因及其危害的基础上，对汽车振动进行了理论分析、MATLAB编程计算以及试验研究。结果表明：将MATLAB强大的数据处理和可视化技术应用于汽车振动分析与控制中，既可以验证理论分析的结果，又可以预测汽车响应，具有很高的实用价值。 关键字：MATLAB 汽车振动激励阻尼自由度 1 汽车振动的产生原因及其危害 机械振动是一种特殊形式的运动，激励、质量、弹性和阻尼是振动系统的四大要素。如果把汽车作为一个系统来研究，汽车本身就是一个具有质量、弹簧和阻尼的振动系统。由于汽车内部各部分的固有频率不同，汽车在行驶中常因路面不平、车速和运动方向的变化，车轮、发动机和传动系统的不平衡，以及齿轮的冲击等各种外部和内部的激振作用而极易产生整车和局部的强烈振动。汽车的这种振动使汽车的动力性得不到充分的发挥，使经济性变坏。同时，还要影响汽车的通过性、操纵稳定性和平顺性，使乘员产生不舒服和疲乏的感觉，甚至损坏汽车的零部件和运载的货物，缩短汽车的使用寿命[1]。因此，研究汽车振动的目的主要有两方面：一是降低振动对汽车零部件的损伤、对汽车使用性能的危害；二是试图利用振动为汽车设计服务，利用振动机理设计制造振动机械以减轻劳动强度，提高工作效率是不乏先例的。 2汽车振动的理论分析 2.1建立振动的力学模型 当一个实际振动系统较复杂时，建立的模型越复杂，越接近实际情况，模拟越逼真，但往往使分析困难；相反地，建立模型时，分析越容易，但得到的结果可能不精确，因此，在建立振动系统力学模型中，总是在求得简化表达和逼真模拟二者之间的折衷[2]。振动分析的关键就是：根据研究的内容和要求，把所研究的对象以及外界对它的作用简化为一个既简单又能在动态特性方面与原来的研究对象等效的力学模型。 汽车是由多个系统组成的复杂的振动系统，每个系统都存在振动问题。主要包括发动机、传动系统、制动系统、转向系统、悬架系统、车身和车架系统存在的振动问题。研究汽车这样一个复杂的振动系统，要根据所分析的问题进行简化，具体简化方案有以下几种 [3]： 1)当汽车对称于其纵轴线时，汽车车身只有垂直振动和俯仰振动对平顺性影响最大。此时，将汽车简化成如图1a)所示的四个自由度的平面模型，因轮胎阻尼较小，在此予以忽略。在这个模型中，车身质量m2，主要考虑垂直和俯仰两个自由度前、后车轴质量m1f，m1r。有两个垂直自由度。 2)当汽车前、后轴悬架质量分配达到一定值时，前、后悬架系统的垂直振动几乎是独立的。于是可以将汽车进一步简化为如图1b)所示的车身和车轮两个自由度振动系统。 M2：为簧载(车身)质量，m1为非簧载(车轮)质量。分析平顺性时，只考虑两个质量的垂直自由度。 信号与系统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： 举例2：MATLAB编程作业
matlab 作业
MATLAB第一章作业答案
MATLAB大作业
MATLAB程序设计作业
MATLAB 作业
MATLAB第一次实验
MATLAB大作业
matlab第一次实验报告
Matlab大作业
Matlab实验第一次实验答案
matlab与数学实验大作业
自动控制原理Matlab程序作业(精)
MATLAB作业
Matlab程序设计大作业(终审稿)
MATLAB期末大作业模版
matlab程序设计作业
MATLAB期末作业(MATLAB应用)
MATLAB实验报告(1-4)
相关主题
文本预览