matlab实验题目.doc

合集下载

MATLAB概率习题

数学实验（概率论）题目一．用MATLAB 计算随机变量的分布1．用MA TLAB 计算二项分布在一级品率为0.2的大批产品中，随机地抽取20个产品，求其中有2个一级品的概率。

1．用MA TLAB 计算泊松分布用MATLAB 计算：保险公司售出某种寿险保单2500份.已知此项寿险每单需交保费120元,当被保人一年内死亡时,其家属可以从保险公司获得2万元的赔偿(即保额为2万元).若此类被保人一年内死亡的概率0.002,试求:(1)保险公司的此项寿险亏损的概率;(2)保险公司从此项寿险获利不少于10万元的概率; (3)获利不少于20万元的概率. 3．用MA TLAB 计算均匀分布乘客到车站候车时间ξ()0,6U ，计算()13P ξ<≤。

4．用MA TLAB 计算指数分布用MA TLAB 计算：某元件寿命ξ服从参数为λ（λ=11000-）的指数分布.3个这样的元件使用1000小时后，都没有损坏的概率是多少? 5。

用MATLAB 计算正态分布某厂生产一种设备，其平均寿命为10年，标准差为2年.如该设备的寿命服从正态分布，求寿命不低于9年的设备占整批设备的比例？二．用MATLAB 计算随机变量的期望和方差 1．用MA TLAB 计算数学期望（1）用MATLAB 计算离散型随机变量的期望 1）。

一批产品中有一、二、三等品、等外品及废品5种，相应的概率分别为0.7、0.1、0.1、0.06及0.04,若其产值分别为6元、5.4元、5元、4元及0元.求产值的平均值 2）。

已知随机变量X 的分布列如下：{}kk X p 21== ,,2,1n k =计算.EX （2）用MATLAB 计算连续型随机变量的数学期望假定国际市场上对我国某种商品的年需求量是一个随机变量ξ（单位：吨），服从区间[],a b 上的均匀分布，其概率密度为： 1()0a x bx b aϕ⎧≤≤⎪=-⎨⎪⎩其它计算我国该种商品在国际市场上年销售量的期望.ξE .（3）用MATLAB 计算随机变量函数的数学期望假定国际市场每年对我国某种商品的需求量是随机变量X （单位：吨），服从[20,40]上的均匀分布，已知该商品每售出1吨，可获利3万美元，若销售不出去，则每吨要损失1万美元，如何组织货源，才可使收益最大？ 2．用MA TLAB 计算方差（1）利用MATLAB 计算：设有甲、乙两种股票，今年的价格都是10元，一年后它们的试比较购买这两种股票时的投资风险.。

MATLAB实训实验

2015/2016学年下学期《信号与系统》实验报告班级:学号:学生姓名:指导教师:2016年3月8 日实验一基本函数仿真实验项目: 基本函数仿真实验时间: 2016年 3 月 8 日星期二第 34 节课实验地点: 1501实验室实验目的:1、学习使用MATLAB 软件2、学习MATLAB 中各种函数，并应用函数分析3、对MATALB 的进一步的学习了解，熟练掌握MATALB 的各种操纵，学会使用MATALB 解决复杂的运算并学会用MATALB 解决平时学习4、了解MATALB 的数值运算5、了解MATALB 的基本函数和命令6、学习掌握MATALB 有关命令实验内容: 1、(1) 题目：应用MA TLAB 方法实现单位阶跃信号和矩形脉冲。

(2) 程序清单（源程序）解：对于阶跃函数，MATLAB 中有专门的stairs 绘图命令。

例如，实现)(t 和矩形脉冲的程序如下：t=-1:2; % 定义时间范围向量t x=(t>=0);subplot(1,2,1),stairs(t,x);axis([-1,2,-0.1,1.2]); grid on % 绘制单位阶跃信号波形 t=-1:0.001:1; % 定义时间范围向量t g=(t>=(-1/2))-(t>=(1/2));subplot(1,2,2),stairs(t,g);axis([-1,1,-0.1,1.2]); grid on % 绘制矩形脉冲波形(3) 运行结果（截图）00.20.40.60.8100.20.40.60.81图1 例1图（4）函数解析Subplot:使用方法：subplot （m,n,p ）或者subplot （m n p ）。

是将多个图画到一个平面上的工具。

其中，m 表示是图排成m 行，n 表示图排成n 列，也就是整个figure 中有n 个图是排成一行的，一共m 行，如果m=2就是表示2行图。

实验一 Matlab基本操作及运算

实验一Matlab基本操作及运算一、实验目的：1．熟悉MATLAB基本操作2．掌握矩阵、变量、表达式的各种基本运算二、实验基本知识：1.熟悉MATLAB环境熟悉MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器文件和搜索路径浏览器。

2.掌握MATLAB常用命令3.MATLAB变量与运算符变量命名规则如下：（1）变量名可以由英语字母、数字和下划线组成（2）变量名应以英文字母开头（3）长度不大于31个（4）区分大小写MATLAB中设置了一些特殊的变量与常量，列于下表。

MATLAB运算符，通过下面几个表来说明MATLAB的各种常用运算符表2 MATLAB算术运算符表3 MATLAB关系运算符表4 MATLAB逻辑运算符表5 MATLAB特殊运算4. MATLAB的一维、二维数组的访问表6 子数组访问与赋值常用的相关指令格式5. MATLAB的基本运算表7 两种运算指令形式和实质内涵的异同表6.MATLAB的常用函数表8 标准数组生成函数表9 数组操作函数7.多项式运算poly——产生特征多项式系数向量roots——求多项式的根p=poly2str(c,‘x’)—（将特征多项式系数向量c转换为以习惯方式显示是多项式）conv, convs——多项式乘运算deconv——多项式除运算polyder(p)——求p的微分polyder(a, b)——求多项式a,b乘积的微分[p,q]=polyder(p1,p2)——求解多项式p1/p2微分的有理分式poly(p,A)——按数组运算规则求多项式p在自变量A的值polym(p,A)——按矩阵运算规则求多项式p在自变量A的值三、实验内容1、新建一个文件夹（自己的名字命名）2、启动MATLAB，将该文件夹添加到MATLAB路径管理器中。

方法如下：3、保存，关闭对话框（要求抓取自己实验的图，插入到自己的实验报告中）4、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye（其它不会用的指令，依照此方法类推）5、使用clc 、clear ，观察command window 、command history 和workspace等窗口的变化结果。

matlab简单编程21个题目及答案

1、设⎥⎦⎤⎢⎣⎡++=)1(sin35.0cos2xxxy，把x=0~2π间分为101点，画出以x为横坐标，y为纵坐标的曲线。

第一题的matlab源程序：①考虑cos（x）为一个整体，然后乘以中括号里面的全部x=0:2*pi/100:2*pi; %x的步长以及范围从0到2*pi y=cos(x).*(0.5+3*sin(x)./(1+x.^2)); %y的表达式plot(x,y)%画出图形图如下：②考虑对整体求解cos，先求x乘以括号中的部分x=0:2*pi/100:2*pi; %x的步长以及范围从0到2*pi y=cos(x.*(0.5+3*sin(x)./(1+x.^2))); %y的表达式plot(x,y) %画出图形图如下：2、产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。

并求该矩阵全体数的平均值和均方差。

第二题的matlab源程序如下：R1=randn(8,6) %产生正态分布随机矩阵R1 =1.0933 -0.7697 1.5442 -0.1924 1.4193 0.21571.1093 0.3714 0.0859 0.8886 0.2916 -1.1658-0.8637 -0.2256 -1.4916 -0.7648 0.1978 -1.14800.0774 1.1174 -0.7423 -1.4023 1.5877 0.1049-1.2141 -1.0891 -1.0616 -1.4224 -0.8045 0.7223-1.1135 0.0326 2.3505 0.4882 0.6966 2.5855-0.0068 0.5525 -0.6156 -0.1774 0.8351 -0.66691.5326 1.1006 0.7481 -0.1961 -0.2437 0.1873aver=(sum(R1(1:end,1:end)))./8 %产生各行的平均值aver =0.0768 0.1363 0.1022 -0.3473 0.4975 0.1044a=std(R1(1:end,1:end)) %产生各行的均方差也就是标准差a =1.0819 0.8093 1.3456 0.8233 0.8079 1.2150aver1=(sum(R1(:)))./48 %全体数的平均值aver1 =0.0950b=std(R1(:)) %全体数的均方差即标准差b =1.01033、设x=rcost+3t,y=rsint+3,分别令r=2,3,4，画出参数t=0~10区间生成的x~y 曲线。

(完整word版)噪声干扰信号的Matlab仿真

雷达对抗实验报告实验题目:噪声干扰信号的Matlab仿真院系：电子与信息工程学院班级：姓名：学号:指导教师：实验时间： 2012 年 6 月噪声调幅、调频、调相信号的Matlab仿真一、实验目的通过实验，加深对噪声调幅、调频、调相信号的理解，加深对噪声调幅、调频、调相信号频谱分析的基本思想与实现方法的认识，并掌握Matlab对随机过程的仿真方法与其基本函数和语法的使用。

二、实验原理实验中要仿真的各种噪声的时域表达式及相应的频谱特性:1.射频噪声干扰窄带高斯过程:称为射频噪声干扰。

其中包络函数服从瑞利分布，相位函数服从［0，2］均匀分布，且与相互独立，载频为常数,且远大于的谱宽。

2.噪声调幅干扰广义平稳随机过程：称为噪声调幅干扰。

其中，调制噪声为零均值，方差为,在区间［—,分布的广义平稳随机过程，服从［0,2]均匀分布，且为与独立的随机变量,为常数。

噪声调幅信号的波形图,以及联合概率密度分布函数p（)以及各自的概率密度分布密度p(）存在下列关系：3.噪声调频干扰广义平稳随机过程：称为噪声调频干扰，其中调制噪声为零均值、广义平稳的随机过程，服从［0，2］均匀分布且与独立的随机变量,，噪声调频干扰中的调制噪声和噪声调频干扰信号的波形J（t）如下图示:4.噪声调相干扰广义平稳随机过程：称为噪声调频干扰,其中调制噪声为零均值、广义平稳的随机过程，服从［0,2］均匀分布且与独立的随机变量，，噪声调相干扰的功率谱如下图所示：三、实验内容利用Matlab仿真产生视频噪声：;射频噪声:；噪声调幅干扰：视频噪声，调制度m=0.1~1;噪声调频干扰：视频噪声；噪声调相干扰：视频噪声.等一系列干扰信号并分析特性。

四、实验思路与步骤1.产生一个高斯白噪声，2.利用Matlab自带的fir1函数产生一个低通滤波器，限制高斯白噪声的带宽，由此产生了视频噪声.3.利用产生的视频噪声，分别代入噪声调幅干扰的时域表达式，并且进行100次的积累后求平均值,由此画出噪声调幅干扰频域波形，对其进行快速傅里叶变换后，求出功率谱，由此画出噪声调幅干扰的功率谱波形。

MATLAB实验题答案

1、求以下变量的值，并在MATLAB 中验证。

( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 82、下列运算是否合法，为什么？如合法，结果是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 25 8 11 >> result4=b*d result4 =31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13-5 -8 7>> result6=a.*b result6 =2 8 -34 15 30 >> result7=a./b result7 =0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000>> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 =2 4 816 32 64 3、用MATLAB 求解下面的的方程组。

(1)⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⨯⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----017413231511222315921274321x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2)⎪⎪⎩⎪⎪⎨⎧=-++=--=-++=++56533332821w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6]>> B2=[1;8;3;5] >> x2=inv(A1)*B24、已知⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=1323151122231592127A (1)求矩阵A 的秩(rank) (2)求矩阵A 的行列式(determinant) (3)求矩阵A 的逆(inverse) (4)求矩阵A 的特征值及特征向量(eigenvalue and eigenvector) >> A3=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> r=rank(A3) >> b=inv(A3) >> a=det(A3) >> [V,D]=eig(A3) 5、109910101022222++++==---=∑n n y ,求y=？（运行format long g 命令后，查看y 的值） m1=0; for m=-10:10 m1=m1+2^m; end m1 m1 = 2047.9990234375 6、求分段函数的值。

MATLAB实验

MATLAB实验一：MATLAB语言基本概念实验实验目的：1. 熟悉MATLAB语言及使用环境；2.掌握MATLAB的常用命令；3．掌握MATLAB的工作空间的使用；4.掌握MATLAB的获得帮助的途径。

5.掌握科学计算的有关方法，熟悉MATLAB语言及其在科学计算中的运用；6.掌握MATLAB的命令运行方式和M文件运行方式；7.掌握矩阵在MATLAB中的运用。

实验方案分析及设计：本次实验主要目的是了解MATLAB的使用环境，以及常用的一些命令的使用；了解矩阵在MATLAB实验中的具体运用，以及相关的一些符号命令的使用。

实验器材：电脑一台，MATLAB软件实验步骤：打开MATLAB程序，将实验内容中的题目依次输入MATLAB中，运行得到并记录结果，最后再对所得结果进行验证。

实验内容及要求：1．熟悉MATLAB的菜单和快捷键的功能2．熟悉MATLAB的命令窗口的使用3．熟悉常用指令的使用format clc clear help lookfor who whos 4．熟悉命令历史窗口的使用5. 熟悉MATLAB工作空间的功能将工作空间中的变量保存为M文件，并提取该文件中的变量6.熟悉MATLAB获取帮助的途径将所有plot开头的函数列出来，并详细给出plotfis函数的使用方法1. 输入 A=[7 1 5;2 5 6;3 1 5]，B=[1 1 1; 2 2 2;3 3 3]，在命令窗口中执行下列表达式，掌握其含义：A(2, 3) A(:,2) A(3,:) A(:,1:2:3)A(:,3).*B(:,2) A(:,3)*B(2,:) A*BA.*BA^2 A.^2 B/A B./AA=[7 1 5;2 5 6;3 1 5]7 1 52 5 63 1 5>> B=[1 1 1; 2 2 2;3 3 3]1 1 12 2 23 3 3>> A(2, 3)6>> A(:,2)151>> A(3,:)3 1 5>> A(:,1:2:3)7 52 63 5>> A(:,3).*B(:,2)51215>> A(:,3)*B(2,:)10 10 1012 12 1210 10 10>> A*B24 24 2430 30 3020 20 20>> A.*B7 1 54 10 129 3 15>> A^266 17 6642 33 7038 13 46>> A.^249 1 254 25 369 1 25>> B/A0.1842 0.2105 -0.23680.3684 0.4211 -0.47370.5526 0.6316 -0.7105>> B./A0.1429 1.0000 0.20001.0000 0.4000 0.33331.0000 3.0000 0.60002．输入 C=1:2:20，则 C （i ）表示什么？其中 i=1,2,3, (10)1到19差为2，i 代表公差3. 试用 help 命令理解下面程序各指令的含义：cleart =0:0.001:2*pi;subplot(2,2,1);polar(t, 1+cos(t))subplot(2,2,2);plot(cos(t).^3,sin(t).^3)subplot(2,2,3);polar(t,abs(sin(t).*cos(t)))subplot(2,2,4);polar(t,(cos(2*t)).^0.5)4计算矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535与⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242之和。

(完整word版)Matlab数学实验报告

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 onfor r=0:0.3:3.9x=[0.1];for i=2:150x(i)=r*sin(3.14*x(i-1));endpause(0.5)for i=101:150plot(r,x(i),'k.');endtext(r-0.1,max(x(101:150))+0.05,['\it{r}=',num2str(r)]) end加密迭代后clear;clf;axis([0,4,0,4]);hold onfor r=0:0.005:3.9x=[0.1];for i=2:150x(i)=r*sin(3.14*x(i-1));endpause(0.1)for i=101:150plot(r,x(i),'k.');endend运行后得到Feigenbaum图2.2实验题目二2.2.1实验问题某农夫有一个半径10米的圆形牛栏，长满了草。

他要将一头牛拴在牛栏边界的桩栏上，但只让牛吃到一半草，问拴牛鼻子的绳子应为多长？2.2.2问题分析如图所示，E为圆ABD的圆心，AB为拴牛的绳子，圆ABD为草场，区域ABCD为牛能到达的区域。

问题要求区域ABCD等于圆ABC的一半，可以设BC等于x,只要求出∠a和∠b就能求出所求面积。

matlab微积分实验题目

matlab微积分实验题目
以下是一些可能的MATLAB微积分实验题目：
1.计算函数f(x)=x^3在区间[0,2]上的定积分，并绘制积分曲线的图形。

2.计算函数f(x)=sin(x)在区间[0,π]上的不定积分，并绘制积分曲线的图形。

3.计算函数f(x)=x^2在区间[1,5]上的导数，并绘制导数曲线的图形。

4.计算函数f(x)=sin(x)/x在区间[1,10]上的极限值，并绘制极限曲线的图形。

5.计算函数f(x)=x^3在区间[0,1]上的极值点，并绘制极值点曲线的图形。

6.计算函数f(x)=x^2在区间[-1,1]上的拐点，并绘制拐点曲线的图形。

7.计算函数f(x)=sin(x)/x在区间[0,π]上的零点，并绘制零点曲线的图形。

8.计算函数f(x)=ln(x)在区间[1,10]上的不定积分，并绘制积分曲线的图形。

9.计算函数f(x)=cos(x)/x在区间[0,π]上的导数，并绘制导数曲线的图形。

10.计算函数f(x)=x^3在区间[-2,2]上的二阶导数，并绘制二阶导数曲线的图形。

1。

信号与系统matlab实验习题3 绘制典型信号及其频谱图

绘制典型信号及其频谱图答案在下面四个常用信号及其傅里叶变换式如表1所示。

(1)绘制单边指数信号及其频谱图的MATLAB程序如下：close all;E=1;a=1;t=0:0.01:4;w=-30:0.01:30;f=E*exp(-a*t);F=1./(a+j*w);plot(t,f);xlabel('t');ylabel('f(t)');figure;plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|');figure;max_logF=max(abs(F));plot(w,20*log10(abs(F)/max_logF));xlabel('\omega');ylabel('|F(\omega)| indB');figure;plot(w,angle(F));xlabel('\omega');ylabel('\phi(\omega)');请更改参数，调试此程序，绘制单边指数信号的波形图和频谱图。

观察参数a 对信号波形及其频谱的影响。

注：题目中阴影部分是幅频特性的对数表示形式，单位是(dB)，请查阅相关资料，了解这种表示方法的意义及其典型数值对应的线性增益大小。

(2)绘制矩形脉冲信号、升余弦脉冲信号和三角脉冲信号的波形图和频谱图，观察并对比各信号的频带宽度和旁瓣的大小。

(3)更改参数，调试程序，绘制单边指数信号的波形图和频谱图。

观察参数a对信号波形及其频谱的影响。

答案附上程序代码：close all;E=1;a=1;t=0:0.01:4;w=-30:0.01:30;f=E*exp(-a*t);F=1./(a+j*w);plot(t,f);xlabel('t');ylabel('f(t)');figure;plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|';E=1，a=1,波形图频谱图更改参数E=2，a=1;更改参数a ，对信号波形及其频谱的影响。

matlab简单编程21个题目及答案

1、设⎥⎦⎤⎢⎣⎡++=)1(sin35.0cos2xxxy，把x=0~2π间分为101点，画出以x为横坐标，y为纵坐标的曲线。

并求该矩阵全体数的平均值和均方差。

MATLAB教程课后实验报告题目及解答[第一至第五章]

MATLAB教程实验报告实验项目名称实验一 Matlab基本操作学生姓名汪德旺专业班级 09数教（1）班学号 0301090131实验成绩日期一. 实验目的和要求1、了解MATLAB 的开发环境。

2、熟悉Matlab的基本操作。

3、掌握建立矩阵的方法。

4、掌握MATLAB各种表达式的书写规则以及常用函数的使用。

5、填写实验报告，实验报告文件取名为report1.doc。

6、于邮件附件形式将实验报告文件report1.doc 发到邮箱*******************，邮件主题为班级学号姓名，如：09数教1班15号张三。

二、实验内容1、先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。

（1）0 122sin851ze =+（2）2212 1ln(0.4552i z x x+⎡⎤==⎢⎥-⎣⎦其中（3）0.30.330.3sin(0.3)ln,22a ae e az a--+=++a=-3.0,-2.9,-2.8,…, 2.8, 2.9,3.0（4）2242,011,12,0:0.5:2.521,23t tz t tt t t⎧≤<⎪=-≤<=⎨⎪-+≤<⎩其中t2.已知：1234413134787,2033657327 A B--⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：（1）A+6*B和A-B+I(其中I为单位矩阵) （2）A*B和A.*B（3）A^3和A.^3（4）A/B和B\A（5）[A,B]和[A([1,3],:);B^2]3、设有矩阵A 和B123453016678910A=,B=17-691112131415023-41617181920970212223242541311⎡⎤⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎣⎦（1）求它们的乘积C 。

（2）将矩阵C 的右下角3*2子矩阵赋给D 。

（3）查看MATLAB 工作空间的使用情况。

MATLAB实验

MATLAB实验实验⼀MATLAB操作基础⼀、实验⽬的1、熟悉MATLAB的操作环境及基本操作⽅法。

2、掌握MATLAB的搜索路径及其设置⽅法。

3、熟悉MATLAB帮助信息的查阅⽅法。

⼆、实验内容1、先建⽴⾃⼰的⼯作⽬录，再将⾃⼰的⼯作⽬录设置到MATLAB搜索路径下，再试验⽤help命令能否查询到⾃⼰的⼯作⽬录。

2、在MATLAB环境下验证例1-1⾄1-4，并完成以下题⽬：（1）绘制右图所⽰图形（23、利⽤MATLAB的帮助功能分别查询inv，plot、max、round等函数的功能及⽤法。

4、在⼯作空间建⽴⼀个变量a ,同时在当前⽬录下建⽴⼀个M⽂件：a.m，试在命令窗⼝输⼊a，观察结果，并解释原因。

三、思考练习1、help命令和lookfor命令有什么区别？2、什么是⼯作空间？假定有变量A与B存在于⼯作空间中，如何⽤命令保存这两个变量？下次重新进⼊MATLAB后，⼜如何装载这两个变量？实验⼆MATLAB数值计算⼀、实验⽬的1、掌握MATLAB变量和数据操作2、掌握MATLAB矩阵及其操作3、掌握MATLAB矩阵运算⼆、实验内容1、求下列表达⽰的值（1）6(10.3424510)w-=+?（2）22tan()b ca eabcxb c aππ++-+=++，其中a=3.5 ,b=5 ,c=-9.8（3）21ln(2tz e t=+，其中213[]50.65it-=-2、已知1540783617A--=??，831253320B-=??-??求下列表达式的值：（1）A+6B和2A B I-+(I为单位矩阵)（2）A*B、A.*B和B*A（3）A/B及B\A（4）[A,B] 和[A([1,3],:) ; B^2]3、建⽴⼀个均值为3，⽅差为1的10*10的正态分布随机矩阵，并将矩阵中⼤于0的元素置1，⼩于0的置0.4、当[34,,,,,,0]A NaN Inf Inf pi eps=--时，求函数all(A),any(A),isnan(A),isinf(A),isfinite(A)的值。

Matlab数学规划方法及实验题目

MATLAB数学规划问题（实验题目及答案在最后）一、线性规划线性规划问题是目标函数和约束条件均为线性函数的问题，MATLAB6.0及更高版本解决的线性规划问题的标准形式为：min n R',f∈xxsub.to：b⋅A≤x⋅Aeq=xbeq≤lb≤xub其中f、x、b、beq、lb、ub为向量，A、Aeq为矩阵。

其它形式的线性规划问题都可经过适当变换化为此标准形式。

在MATLAB6.0版中，线性规划问题（Linear Programming）已用函数linprog取代了MATLAB5.x版中的lp函数。

在6.0和7.0中依然可以使用lp 函数，但在更高版本中，就只能使用linprog函数了。

函数linprog调用格式：x=linprog(f,A,b)x=linprog(f,A,b,Aeq,beq)- 1 -- 1 -x=linprog(f,A,b,Aeq,beq,lb,ub) x=linprog(f,A,b,Aeq,beq,lb,ub,x0) x=linprog(f,A,b,Aeq,beq,lb,ub,x0,options) [x,fval]=linprog(…)[x, fval, exitflag]=linprog(…) [x, fval, exitflag, output]=linprog(…)[x, fval, exitflag, output, lambda]=linprog(…) 说明：x=linprog(f, A, b) %求min f ' *x, sub.to b x A ≤⋅线性规划的最优解。

返回值x 为最优解向量。

x=linprog(f, A, b, Aeq, beq) %含有等式约束beq x Aeq =⋅，若没有不等式约束b x A ≤⋅，则令A=[ ]，b=[ ]。

x = linprog(f, A, b, Aeq, beq, lb, ub) %指定x 的范围ub x lb ≤≤ x=linprog(f, A, b, Aeq, beq, lb, ub, x0) %设置x0为初值点。

MATLAB仿真实例

一、实验目的 (1)二、实验题目 (1)三、实验内容 (1)3.1傅里叶变换与傅里叶反变换 (1)3.2题目一：正弦信号波形及频谱 (2)3.2.1仿真原理及思路 (2)3.2.2程序流程图 (3)3.2.3仿真程序及运行结果 (3)3.2.4实验结果分析 (5)3.3题目二：单极性归零（RZ）波形及其功率谱 (5)3.3.1仿真原理及思路 (5)3.3.2程序流程图 (6)3.3.3仿真程序及运行结果 (6)3.3.4实验结果分析 (8)3.4题目三：升余弦滚降波形的眼图及其功率谱 (8)3.4.1仿真原理及思路 (8)3.4.2程序流程图 (8)3.4.3仿真程序及运行结果 (8)3.4.4实验结果分析： (10)3.5题目四：完成PCM编码及解码的仿真 (11)3.5.1仿真原理及思路 (11)3.5.2程序流程图 (12)3.5.3仿真程序及运行结果 (12)3.5.4实验结果分析 (15)3.6附加题一：最佳基带系统的Pe~Eb\No曲线，升余弦滚降系数a=0.5，取样值的偏差是Ts/4 (16)3.6.1仿真原理及思路 (16)3.6.2程序流程图 (16)3.6.3仿真程序及运行结果 (16)3.6.4实验结果分析 (18)3.7附加题二：试作出Pe～Eb/No曲线。

升余弦滚降系数a＝0.5，取样时间无偏差，但信道是多径信道，C(f)=|1-0.5-j2 ft|,t=T s/2 (18)3.7.1仿真原理及思路 (18)3.7.2程序流程图 (19)3.7.3仿真程序及运行结果 (19)3.7.4实验结果分析 (21)四、实验心得 (21)一、实验目的⏹ 学会MATLAB 软件的最基本运用。

MATLAB 是一种很实用的数学软件，它易学易用。

MATLAB 对于许多的通信仿真类问题来说是很合适的。

⏹ 了解计算机仿真的基本原理及方法，知道怎样通过仿真的方法去研究通信问题。

⏹ 加深对通信原理课程有关内容的理解。

(完整版)利用MATLAB进行时域分析

自动控制原理与系统课程实验报告实验题目：利用MATLAB进行时域分析班级：机电1131班姓名：刘润学号：38号一、实验目的及内容时域分析法是一种直接在时间域中对系统进行分析的方法，具有直观、准确的优点，并且可以提供系统时间响应的全部信息。

在此实验中，主要介绍时域法进行系统分析，包括一阶系统、二阶系统以及高阶系统,以及系统的性能指标。

通过实验，能够快速掌握、并利用MATLAB及控制系统箱对各种复杂控制系统进行时域分析。

二、实验设备三、实验原理典型的二阶系统在不同的阻尼比的情况下，它们的阶跃响应输出特性的差异是很大的。

若阻尼比过小，则系统的振荡加剧，超调量大幅度增加；若阻尼比过大，则系统的响应过慢，又大大增加了调整时间，下面通过此实验课题分析输出响应变化规律：已知二阶振荡环节的传递函数为：G（s）=ωn*ωn/（s*s+2*ζ*ωn*s+ωn*ωn），其中ωn=0.4，ζ从0变化到2，求此系统的单位阶跃响应曲线，并分析当ζ发生变化时，二阶系统的响应有什么样的变化规律。

四、实验步骤编出程序如下图：五、实验结果画出图表如下图：六、结果分析（1）当ξ=0（无阻尼）（零阻尼）时：无阻尼时的阶跃响应为等幅振荡曲线。

如图ξ=0曲线。

（2）当0<ξ<1（欠阻尼）时：对应不同的ξ，可画出一系列阻尼振荡曲线，且ξ越小，振荡的最大振幅愈大。

如图ξ=0.4曲线。

(3)当ξ=1（临界阻尼）时：临界阻尼时的阶跃响应为单调上升曲线。

如图ξ=1曲线。

（4）当ξ>1（过阻尼）时：过阻尼时的阶跃响应也为单调上升曲线。

不过其上升的斜率较临界阻尼更慢。

如图ξ=1.6曲线七、教师评语。

MATLAB程序设计插值实验题目

1、当2,1,1-=x 时，4,3,0)(-=x f ，求x =0时的值。

2、给出x x f ln )(=的数值表用线性插值及二次插值计算54.0ln 的近似值。

3、4()31f x x x =++，已知x =0,1,2处的值，计算x =0.11的近似值。

4.求i =3和5时的值。

5. 给定数据表，采用牛顿插值方法求i =3和5时的值。

6. 已知函数y =在4, 6.25,9x x x ===的近似值。

7. 已知函数y =在4,9,16x x x ===的近似值。

8. 已知函数y =在0,1,4,9x x x x ==== 9. 给定数据表：构造出函数()f x 的差商表，并计算x =0.5时的值.10. 已知函数y =在8,27,64x x x ===的近似值。

11. 已知函数sin()y x =在0,/4,/2x x x ππ===处的函数值，试通过一个二次插值函数求sin(/5)π的近似值。

12. 已知函数cos()y x =在0,/4,/2x x x ππ===处的函数值，试通过一个二次插值函数求cos(/5)π的近似值。

13. 已知函数tan()y x =在0,/4,/3x x x ππ===处的函数值，试通过一个二次插值函数求tan(/5)π的近似值。

14. 已知函数4y x =在1,2,3x x x ===处的函数值，试通过一个二次插值函数求 1.111x =的近似值。

15. 已知函数sin()y x =在0,/4,/2x x x ππ===处函数值，通过一个二次插值函数求1sin(2/5)π-的近似值。

16. 已知函数cos()y x =在0,/4,/2x x x ππ===处函数值，通过一个二次插值函数求1+cos(/5)π的近似值。

17. 已知函数5y x =在1,2,3x x x ===处的函数值，试通过一个二次插值函数求 1.121x =的近似值。

18. 已知函数41/y x =在1,2,3x x x ===处的函数值，试通过一个二次插值函数求 1.211x =的近似值。

matlab实验3

一实验题目：实验三抽样二实验原理：抽样过程两个基本原理：混叠和重建。

由于在MATLAB 中不能产生模拟信号，所以需要做实轴t 轴的仿真。

因此，把仿真时的△t 与所研究的抽样周期s T 明确的放开是很重要的对连续时间正弦信号考虑下面的表达式：0()sin (2)x t f t πθ=+可以按抽样频率1/s s f T =对x(t)抽样来获得离散时间信号/0[]()|()|sin(2/)sst nT t n f s x n x t x t f n f πθ=====+对连续时间信号抽样时，因为频域以抽样频率延拓，故其频谱显现出混叠效应。

在MATLAB 中，只能仿真这一效应。

仿真包括抽样运算，D/A 转换（包括一个重建滤波器）。

为了仿真模拟信号，必须用非常高的采样率—至少是任何模拟信号所允许的具有的最高频率的5倍。

测试混叠。

三实验内容： 1.以不同组合的0f 和 s f 绘出x[n], s f =8kHz, 0f =300Hz,在10ms 长间隔上采样，相位θ任意指定。

a. 使用stem 绘出产生的离散信号b. 使用plot 绘图c. 把正弦频率从100Hz 变至475Hz ，每次增加125Hz,用subplot 指令把四个图放在一个屏上。

d. 把正弦频率从7525Hz 变至7900Hz ，每次增加125Hz,用subplot 指令把四个图放在一个屏上。

e. 把正弦频率从32,100Hz 变至32，475Hz ，每次增加125Hz,用subplot 指令把四个图放在一个屏上。

提前预测显现的频率将会增加还是减少。

2.用M 文件绘制模拟频谱3.a.产生一个模拟频率为0f 的余弦信号作为仿真的模拟信号△0()cos(2),0x t f t t T πθ=+≤≤在长度为T 的时间间隔上，以速率sim f 产生抽样样本。

选择仿真模拟信号长度T 以取得大约900到1000个样本。

b.用plot 指令将绘出的时间信号样本连接起来，确定用时间连续时间标记时间轴。

相关主题

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

前言一、上机实验目的上机实验的目的是提高学生对算法的理解程度，并掌握用实用工具进行数值计算的方法，通过实践环节理解数值分析的应用和研究方法。

二、实验基本内容本课程实验内容分为6个实验。

学生可以在课内机时先完成指导书中给出的程序或验证算法，理解所学的知识，在此基础上再编写其他应用程序。

指导书中的15个实验如下。

1．拉格朗日插值法。

2．最小二乘拟合。

3．数值积分。

4．范数计算和LU分解。

5．牛顿迭代法。

6．龙格库塔法。

三、实验任务与时间安排本课程是一门实践性很强的课程，除了在课内安排的实验外，鼓励同学在课外用相关技术进行编程练习。

具体实践课时安排如下：实验指导书实验1 拉格朗日插值法一、实验目的[1] 了解lagrange插值法的基本原理和方法；[2] 通过实例掌握用MA TLAB求插值的方法；[3] 编程实现lagrange插值二、实验内容1．在matlab中command窗口或新建M文件运行以下程序>> x0=[0.4 0.5 0.6 0.7 0.8];%输入节点的X值>> y0=[-0.916291 -0.693147 -0.510826 -0.356675 -0.223144];%输入节点的Y值>> T=interp1(x0,y0,0.54,'linear');%对插值节点使用线性插值并计算0.54处的函数值，改变红色处的参数可得到不同的插值结果，可选项为’cubic’,’ nearest’,’ spline’等。

>> x=0.2:0.05:1.0;%从0.2开始到1.0每隔0.05取一个待求点>> y=interp1(x0,y0,x,'linear');%用插值法求出待求点对应的函数值y>> plot(x0,y0,’r*’);%用红色星号标出插值点对>> hold on;plot(x,y);%蓝色线为所得的插值函数图形2．用C语言或matlab编写实现lagrange插值，要求：数据输入项（函数参数）为：插值节点及函数值，及待求点x的值输出为待求点x对应的函数值（程序流程图如下所示）三、实验步骤1．开启软件平台——MA TLAB，编程可以选用自己熟悉的软件；2．根据各种数值解法步骤编写M文件3．观察运行结果(数值或图形)；四、实验要求与任务每个同学独立完成编程，并在课后提交程序；同时验证matlab中相应函数的用法。

实验2 最小二乘拟合一、实验目的[1] 了解最小二乘拟合的基本原理和方法；[2] 掌握用MA TLAB作最小二乘多项式拟合和曲线拟合的方法；[3] 通过实例学习如何用拟合方法解决实际问题，注意与插值方法的区别。

[4] 了解各种参数辨识的原理和方法；[5] 通过范例展现由机理分析确定模型结构，拟合方法辨识参数，误差分析等求解实际问题的过程；二、实验内容1．用MA TLAB 中的函数作一元函数的多项式拟合与曲线拟合； 2．用MA TLAB 中的函数作二元函数的最小二乘拟合，作出误差图； 3．针对预测和确定参数的实际问题，建立数学模型，并求解。

三、实验步骤1．开启软件平台——MA TLAB ，开启MA TLAB 编辑窗口； 2．根据各种数值解法步骤编写M 文件 3．观察运行结果(数值或图形)；四、实验要求与任务根据实验内容和步骤，完成以下具体实验旧车价格预测某年美国旧车价格的调查资料如下表，其中x i 表示轿车的使用年数，y i 表示相应的平均价格。

试分析用什么形式的曲线来拟合上述的数据，并预测使用4.5年后轿车的平均价格大致为多少？多项式函数拟合：[a,s]=polyfit(xdata,ydata,n)其中n 表示多项式的最高阶数，xdata ，ydata 为将要拟合的数据，它是用数组的方式输入．输出参数a 为拟合多项式11++++=n n na x a x a y Λ的系数],,,[11+=n n a a a a Λ,s 为误差。

另：多项式在x 处的值y 可用下面程序计算．y=polyval(a,x)，对上面给出的数据做多项式拟合，可取不同的n 观测此时的误差，看取什么样的n 较好。

并在matlab 中画出拟合曲线(取拟合曲线上若干个点对，用plot 命令实现)实验3 数值积分一、实验目的[1] 了解数值积分的基本原理和方法；[2] 掌握用MATLAB求积分的方法；[3] 通过实例学习如何用几种方法求积分。

[4] 了解书上介绍的几种数值积分的不同原理和方法；二、实验内容1．用MA TLAB中的QUAD函数求出作业所示函数的数值积分；2．用MA TLAB中的QUADL函数求出作业中所示函数的数值积分。

三、实验步骤1．开启软件平台——MA TLAB，开启MA TLAB编辑窗口；2．根据各种数值解法步骤编写M文件3．观察运行结果(数值或图形)；四、实验要求与任务根据实验内容和步骤，完成以下具体实验根据提示计算积分dx e x ⎰-12实验提示：1．Q uad函数的使用方法，参数意义Q = QUAD(FUN,A,B)Fun为待求积的函数A,B为积分区间Quad函数是使用误差为10-6的递归自适应Simpson公式计算Fun函数在A，B区间上的积分值。

例子（摘自帮助）：求dx x x ⎰--2035211）首先，定义myfun 函数为待求积函数function y = myfun(x)y = 1./(x.^3-2*x-5); 2）计算数值积分Q = quad(@myfun,0,2);2． Q uad8函数的使用方法，参数意义Q = QUAD8(FUN,A,B) Fun 为待求积的函数 A,B 为积分区间Quad8函数是使用误差为10-3的递归自适应Newton-Cotes 8公式计算Fun 函数在A ，B 区间上的积分值。

它是收敛速度比较高的一种方法，有些版本已经将此函数改为QUADL具体使用方法同quad3． Q UADL 函数的使用方法，参数意义Q = QUADL(FUN,A,B) Fun 为待求积的函数 A,B 为积分区间Quad8函数是使用自适应Lobatto 求积公式计算Fun 函数在A ，B 区间上的积分值。

它也是收敛速度比较高的一种方法，具体使用方法同quad注意，使用上面三个函数的时候A ，B ，Fun 函数的输出必须是以下三种类型的 float:，double ，single实验4 范数计算和LU 分解一、实验目的[1] 了解矩阵范数的定义和常见的几种矩阵范数； [2] 掌握用直接法解线性方程组的常见几种方法；[3] 通过实例学习如何使用直接法中的LU 分解求解线性方程组。

[4] 掌握判断矩阵，线性方程组本身性质的方法——条件数的计算；二、实验内容1．用MA TLAB 中的NORM 函数求矩阵的范数； 2．用MA TLAB 中的Cond 计算矩阵条件数；3．用MA TLAB 中的Chol 函数计算Cholesky 分解（平方根），Lu 函数计算Doolittle 分解。

三、实验步骤1．开启软件平台——MA TLAB ，开启MA TLAB 编辑窗口； 2．根据各种数值解法步骤编写M 文件 3．观察运行结果(数值或图形)；四、实验要求与任务根据实验内容和步骤，完成以下具体实验实验： 1．用H=Hilb （n ）函数生成一个n 阶hilbert 矩阵H ，（n=4） 2．计算生成的矩阵H 的范数 3．计算H 的条件数 4．⎪⎪⎪⎭⎫⎝⎛=6858137576A 用平方根法分解A5．⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=764142321,461561552621C B ，B ，C 是否存在三角分解，若存在则用LU分解实现。

实验5 牛顿迭代法一、实验目的[1] 了解求解非线性方程的解的常见方法； [2] 掌握用二分法，迭代法解线性方程组；[3] 通过实例学习用牛顿法分解求解非线性方程。

[4] 掌握判断迭代法，牛顿法收敛的条件，和收敛阶的计算；二、实验内容1．用牛顿迭代法求方程0634223=-+-x x x 在区间（1，5）之间的根。

三、实验步骤1．用C 语言或matlab 完成编程 2．观察运行结果实验运行结果；四、实验要求与任务根据实验步骤，完成以下具体实验实验： 1．按照以下流程图完成牛顿迭代法x = newton( fname,dfname,x0,e,N )2．调用牛顿迭代法的函数计算实验内容中第一题。

实验6龙格库塔法一、实验目的[1] 了解matlab中提供的求解常微分方程解的常见方法——龙格库塔法的函数ODE45,ODE32的使用方法；二、实验内容求解⎪⎩⎪⎨⎧=-=1)0(2'x x x三、实验步骤1．在command 窗口中输入help ode32 或者help ode45 查看相关的函数调用格式和使用算法的情况2．用函数求解实验内容的原函数，并画相应的图。

四、实验要求与任务根据实验步骤，先观察帮助里面提供的例子，然后再求解实验内容。