matlab(二)实验1汇总
MATLAB实验报告
MATLAB实验报告实验报告课程名称 MATLAB基础及应用专业班级电子xxxx姓名学号电气与信息学院实验二 MATLAB 数值计算(一)一实验目的:1.掌握数组的创建与运算方法;2. 掌握矩阵的创建与运算方法;3. 掌握数组的运算方法和矩阵运算方法的区别;4.掌握线性方程的求解方法二实验装置:计算机三实验内容:1.数组的创建和运算创建两个含5个元素的一维数组,并求这两个数组的四则运算。
2.矩阵的创建和运算(1)创建两个3×3的矩阵,并求这两个矩阵的四则运算。
(2)创建一个4×4的矩阵,并求这个矩阵的行列式值。
(3)线性方程的求解解方程组=6613753467294x 。
四实验要求:写出实验程序与仿真结果。
1. >> a=linspace(2,10,5) a =2 4 6 8 10>> b=linspace(1,9,5) b =1 3 5 7 9>> a+bans =3 7 11 15 19>> a-bans =1 1 1 1 1>> a.*bans =2 12 30 56 90>> a./bans =2.0000 1.3333 1.2000 1.1429 1.1111 >> a.\b ans =0.5000 0.7500 0.8333 0.8750 0.9000 2.(1) >> a=[1 1 1;2 2 2;3 3 3];>> b=[4 4 4;5 5 5;6 6 6];>> a+bans =5 5 57 7 79 9 9>> a-bans =-3 -3 -3-3 -3 -3-3 -3 -3>> a*bans =15 15 1530 30 3045 45 45>> a/bWarning: Matrix is singular to working precision. ans =NaN NaN NaNNaN NaN NaNNaN NaN NaN>> a\bWarning: Matrix is singular to working precision. ans =NaN NaN NaNNaN NaN NaNInf Inf Inf(2).>> a=[1 1 1 1;2 2 2 2;3 3 3 3;4 4 4 4];>> det(a)ans =(3).>> a=[4 9 2;7 6 4;3 5 7];>> b=[13;6;6];>> x=inv(a)*bx =-0.53081.7109-0.1374实验三 MATLAB数值计算(二)一实验目的:1.掌握多项式的创建与运算方法;2. 掌握基本的数据分析方法;二实验装置:计算机三实验内容:1.多项式创建输入系数矢量,创建多项式x^3-2*x^2+5*x+3。
Matlab第二次实验报告
Matlab第二次实验报告实验目的:1、了解plot函数和subplot函数的基本用法和matlab绘图的基本原理。
2、了解图形的属性设置。
比如画图的颜色,画图采用的线性标识符等。
二:实验基本知识1:1. 单窗口单曲线绘图;2. 单窗口多曲线绘图;3. 单窗口多曲线分图绘图;4. 多窗口绘图;5.可任意设置颜色与线型;6.图形加注功能;7.fplot——绘制函数图函数;8.ezplot——符号函数的简易绘图函数9:subplot—将画图区域分块函数。
实验内容:将高数课本后的18个图画出上机练习程序第1——6图:clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1);x=-50:1:50;y=a.*x.^3;subplot(3,2,1);plot(x,y,'r');title('y=a*x^3');xlabel('X');ylabel('Y');x=0:1:50;y=sqrt(a.*x.^3);subplot(3,2,2);plot(x,y,x,-y,'r');title('y^2=a*x^3');xlabel('X');ylabel('Y');x=-3:0.1:3;y=exp(-x.^2);subplot(3,2,3);plot(x,y,'r');title('y=e^x^2');xlabel('X');ylabel('Y');x=-3:0.1:3;y=8*a^3./(x.^2+4*a^2);subplot(3,2,4);plot(x,y,'r');title('y=8*a^3./(x.^2+4*a^2)'); xlabel('X');ylabel('Y');x=0:0.01:5;y=sqrt(x.^3./(a*2.-x)); subplot(3,2,5);plot(x,y,x,-y,'r');title('y^2*(2a-x)=x^3'); xlabel('X');ylabel('Y');%t=-1:0.01:5;%x=3*a.*t./(1+t.^3);%y=3*a*t.^2./(1+t.^3);subplot(3,2,6);ezplot('x.^3+y.^3-3*3*x.*y'); %plot(x,y);title('x^3+y^3-3axy=0'); xlabel('X');ylabel('Y');输入a=3:显示结果第7——10图:clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1)t=0:0.1:2*pi;x=a*cos(t).^3;y=a*sin(t).^3;subplot(2,2,1);plot(x,y,'b')title('x^1.5+y^1.5=a^1.5'); xlabel('X');ylabel('Y');t=-2*pi:0.1:2*pi;x=a.*(t-sin(t));y=a.*(1-cos(t));subplot(2,2,2);plot(x,y,'b')title('°ÚÏß');xlabel('X');ylabel('Y');t=0:0.01:2*pi;p=a.*(1-cos(t));subplot(2,2,3);plot(p.*cos(t),p.*sin(t),'b'); title('p=a(1-cos£¨t£©)'); xlabel('X');ylabel('Y');t=0:0.01:2*pi;p=a.*t;subplot(2,2,4);plot(p.*cos(t),p.*sin(t),'b'); title('p=at');xlabel('X');ylabel('Y');输入a=1,第11题图:clear,clca=input('ÇëÊäÈëaµÄÖµ£º') figure(1);t=-2*pi:pi/100:2*pi;p=exp(a.*t);plot(p.*cos(t),p.*sin(t)); title('p=e^at');xlabel('X');ylabel('Y');输入a=0.1,显示结果:clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1)t=0.2*pi:pi/100:100*pi;p=a./t;plot(p.*cos(t),p.*sin(t));title('p=at');xlabel('X');ylabel('Y');输入a=1,显示结果:第13——14题图:clear,clca=input('ÇëÊäÈëaµÄÖµ£º');figure(1);t=0:pi/100:2*pi;p=(a^2.*sin(2*t)).^0.5;subplot(1,2,1);plot(p.*cos(t),p.*sin(t),-p.*cos(t),-p.*sin(t),'g'); title('p^2=a^2*sin(2t)');xlabel('X');p=(a^2.*cos(2*t)).^0.5;subplot(1,2,2);plot(p.*cos(t),p.*sin(t),-p.*cos(t),-p.*sin(t),'g'); title('p^2=a^2*cos(2t)');xlabel('X');ylabel('Y');输入a=1,显示结果:第15-16题图:clear;clc;a=input('ÇëÊäÈëa:');figure(1)t=0:pi/200:pi;p=a.*cos(3*t);subplot(1,2,1);plot(p.*cos(t),p.*sin(t),'r--');title('p=a*cos(3t)');xlabel('X');p=a.*sin(3*t);subplot(1,2,2);plot(p.*cos(t),p.*sin(t),'r--'); title('p=a*sin(3t)');xlabel('X');ylabel('Y');输入a=1,显示结果:第17——18题图:clear;clc;a=input('ÇëÊäÈëa:');figure(1);t=-pi:pi/200:pi;subplot(1,2,1);p=a.*sin(2*t);plot(p.*cos(t),p.*sin(t),'k-.'); title('p=a*sin(2t)');xlabel('X');ylabel('Y');p=a.*cos(2*t);subplot(1,2,2);plot(p.*cos(t),p.*sin(t),'k-.');title('p=a*cos(2t)');xlabel('X');ylabel('Y');输入a=1,显示结果:实验心得:这次实验主要是学习matlab软件的二维绘图功能:matlab软件还有丰富的图形修饰功能,如改变线条的形式和颜色;除此之外其还可以以多种形式如在一个图形中显示多个函数图形,一个figure中显示多个坐标系,还有一个程序中就可以用多个figure一起来表示不同的函数……在画图过程中你可能因为角度的区间取值不一样,画出来的图形可能跟书本上给出的图形有很大的差异,但没关系,只要你思路是对的,区间可以自己慢慢改动知道跟书本上给出的图形想接近。
MATLAB实验报告(1-4)
信号与系统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:3.matlab程序流程控制1).for循环:for循环变量=初值:增量:终值循环体End2).while循环结构:while 逻辑表达式循环体End3).If分支:(单分支表达式)if 逻辑表达式程序模块End(多分支结构的语法格式)if 逻辑表达式1程序模块1Else if 逻辑表达式2程序模块2…else 程序模块nEnd4).switch分支结构Switch 表达式Case 常量1程序模块1Case 常量2程序模块2……Otherwise 程序模块nEnd4.典型信号的MATLAB表示1).实指数信号:y=k*exp(a*t)举例:2).正弦信号:y=k*sin(w*t+phi)3).复指数信号:举例:4).抽样信号5).矩形脉冲信号:y=square(t,DUTY) (width默认为1)6).三角波脉冲信号:y=tripuls(t,width,skew)(skew的取值在-1~+1之间,若skew取值为0则对称)周期三角波信号或锯齿波:Y=sawtooth(t,width)5.单位阶跃信号的MATLAB表示6.信号的时移、反折和尺度变换:Xl=fliplr(x)实现信号的反折7.连续时间信号的微分和积分运算1).连续时间信号的微分运算:语句格式:d iff(function,’variable’,n)Function:需要进行求导运算的函数,variable:求导运算的独立变量,n:求导阶数2).连续时间信号的积分运算:语句格式:int(function,’variable’,a,b)Function:被积函数variable:积分变量a:积分下限b:积分上限(a&b默认是不定积分)8.信号的相加与相乘运算9.信号的奇偶分解四、小结这一次实验让我能够教熟悉的使用这个软件,并且能够输入简单的语句并输出相应的结果和波形图,也在一定程度上巩固了c语言的一些语法。
matlab实验总结
matlab实验总结在这次实验中,我使用MATLAB对一些数学问题进行了求解,并学会了如何使用MATLAB编写函数和脚本来解决实际问题。
首先,我了解了MATLAB的基本语法和使用方法。
MATLAB 是一种专业的数学软件,它具有强大的矩阵运算和数据处理的能力。
我学会了如何创建矩阵、向量和标量,并进行基本的数学运算如加减乘除、幂运算等。
我还了解了MATLAB的数据类型和变量命名规则,以及如何使用MATLAB的函数来实现一些特定的功能。
接着,我学习了如何在MATLAB中绘制图形。
MATLAB提供了许多绘图函数,如plot、scatter、bar等,可以用来可视化数据和函数的图像。
我学会了如何设置图形的标题、坐标轴标签、图例以及颜色和线型等属性。
我还学会了如何在同一个图形中绘制多个曲线、设置坐标轴的范围和刻度,并添加文本注释。
在本次实验中,我也学会了如何编写MATLAB的函数和脚本。
函数是一段独立的代码块,用于完成特定的任务,并可以在需要时被调用。
我了解了如何定义输入参数和输出参数,并学会了如何在函数内部使用局部变量和全局变量。
脚本是一系列命令的集合,可以按照一定的顺序依次执行。
我学会了如何在MATLAB编辑器中创建脚本文件,并使用run命令执行脚本。
我还了解了如何在脚本中保存和加载数据,以及如何将数据导出到其他文件格式。
最后,在这次实验中,我也学习了如何使用MATLAB进行数值计算和求解方程。
MATLAB提供了丰富的数值计算函数,如求解线性方程组、解微分方程、拟合曲线等。
我了解了如何使用这些函数来解决实际的数学问题。
我还学会了如何使用MATLAB的符号计算工具箱来进行符号计算,如求解代数方程、简化表达式等。
总的来说,这次实验让我对MATLAB有了更深入的了解,并提高了我的编程和数学计算能力。
通过实践,我不仅复习了数学的一些基本概念和方法,还学会了如何使用MATLAB来加速数学计算和可视化数据。
我相信这些经验和技能将对我的学习和研究工作有很大的帮助。
MATLAB实验报告(四份!!!)
0380160815880001449451490168779983287916534749293454515731] >> double(x) ans = -33.5325 60.6078 -54.7098 27.3192 实验三 MATLAB 绘图
一 实验环境 计算机 MATLAB 软件 二 实验目的 1. 掌握 MATLAB 的基本绘命令。 2. 掌握运用 MATLAB 绘制一维、二维、三维的图形方法。 3. 给图形加以修饰。 三 预备知识 特殊的二维图形函数有:直方图、柄图、阶梯图、饼图、频数累计柱状 图、极坐标图 四 实验内容和步骤
3.创建一个 3x3 矩阵,并求其转置,逆矩阵。 >> A=magic(3) A= 8 3 4 >> A' ans = 8 1 6 >> inv(A) 3 5 7 4 9 2 1 5 9 6 7 2
ans = 0.1472 -0.0611 -0.0194 -0.1444 0.0222 0.1889 0.0639 0.1056 -0.1028
3. 画出 y=的曲线(x(-5,5)) 。在这曲线上加入相同区间里的 y=的曲线,并且采用绿色折 线标识。 >> x=linspace(-5,5,100); >> y=x.^2; >> plot(x,y) >> hold on >> z=x.^(1/3); >> plot(x,z,'g--')
4. 在同一窗口不同坐标系里分别绘出 y1=sinx,y2=cosx,y3=cinh(x),y4=cosh(x)4 个图形。 >> x=linspace(0,2*pi,30); >> subplot(2,2,1);plot(x,sin(x)); >> subplot(2,2,2);plot(x,cos(x)); >> subplot(2,2,3);plot(x,sinh(x)); >> subplot(2,2,4);plot(x,cosh(x));
matlab实验心得总结
matlab实验心得总结在通过完成一系列的Matlab实验后,我对这个强大的数学计算软件有了更深入的认识。
通过这些实验,我不仅学到了如何使用Matlab进行数据处理和分析,还体会到了它在科学研究和工程应用中的广泛使用。
实验一:Matlab基础操作在第一次接触Matlab时,我首先学习了它的基本操作。
Matlab提供了友好的用户界面和丰富的命令工具,使得数据处理变得简单且高效。
在实验中,我学会了如何定义变量、进行基本的数学运算和使用矩阵操作等。
这些基础操作为后续的实验打下了坚实的基础。
实验二:数据可视化数据可视化在科学研究和工程领域中起着重要的作用。
在这个实验中,我学会了如何利用Matlab绘制各种图形,如折线图、散点图和柱状图等。
通过调整图形的样式和颜色,使得数据更加直观和易于理解。
同时,我还学会了如何添加标题、坐标轴标签和图例,使得图形具有更好的可读性。
实验三:模拟与仿真Matlab不仅可以进行数据处理和图形绘制,还可以进行模拟和仿真。
在这个实验中,我学会了如何使用Matlab进行数学模型的建立和仿真。
通过设定合适的参数和方程,我可以模拟出各种现实世界中的物理、生物和工程现象。
这对于科学研究和工程设计具有重要的意义。
实验四:信号处理信号处理是Matlab的一个重要应用领域。
在这个实验中,我学会了如何使用Matlab对信号进行分析和处理。
通过应用不同的滤波器,我可以去除信号中的噪声和干扰,提取出感兴趣的信息。
同时,我还学会了如何进行频域分析,通过傅里叶变换将信号转换到频率域,进一步分析信号的频谱特性。
实验五:数值计算Matlab还提供了强大的数值计算功能。
在这个实验中,我学会了如何使用Matlab进行数值计算和优化。
通过使用不同的数值求解方法,我可以解决复杂的数学方程和优化问题,得到精确的计算结果。
这对于科学研究和工程计算具有重要的价值。
总结起来,通过这些实验,我对Matlab的应用能力有了明显的提升。
MATLAB实验报告
MATLAB实验报告MATLAB实验报告姓名:专业:学号:实验⼀MATLAB环境的熟悉与基本运算⼀、实验⽬的:1.熟悉MATLAB开发环境2.掌握矩阵、变量、表达式的各种基本运算⼆、实验基本知识:1.熟悉MATLAB环境:MATLAB桌⾯和命令窗⼝、命令历史窗⼝、帮助信息浏览器、⼯作空间浏览器⽂件和搜索路径浏览器。
2.掌握MATLAB常⽤命令3.MATLAB变量与运算符变量命名规则如下:(1)变量名可以由英语字母、数字和下划线组成(2)变量名应以英⽂字母开头(3)长度不⼤于31个(4)区分⼤⼩写MATLAB中设置了⼀些特殊的变量与常量,列于下表。
MATLAB运算符,通过下⾯⼏个表来说明MATLAB的各种常⽤运算符表2MATLAB算术运算符表3MATLAB关系运算符表4MATLAB逻辑运算符表5MATLAB特殊运算4.MATLAB的⼀维、⼆维数组的寻访表6⼦数组访问与赋值常⽤的相关指令格式5.MATLAB的基本运算表7两种运算指令形式和实质内涵的异同表6.MATLAB的常⽤函数表8标准数组⽣成函数表9数组操作函数三、实验内容1、学习安装MATLAB软件。
2、学习使⽤help命令,例如在命令窗⼝输⼊helpeye,然后根据帮助说明,学习使⽤指令eye(其它不会⽤的指令,依照此⽅法类推)3、学习使⽤clc、clear,观察commandwindow、commandhistory和workspace等窗⼝的变化结果。
4、初步程序的编写练习,新建M-file,保存(⾃⼰设定⽂件名,例如exerc1、exerc2、exerc3……),学习使⽤MATLAB的基本运算符、数组寻访指令、标准数组⽣成函数和数组操作函数。
注意:每⼀次M-file的修改后,都要存盘。
四、实验结果练习A:(1)helprand,然后随机⽣成⼀个2×6的数组,观察commandwindow、commandhistory和workspace等窗⼝的变化结果。
matlab实验报告
(2)将上述信号左平移或右平移 得到两个新信号y2和y3,利用hold命令在同一图中显示三条曲线并通过线型和颜色加以区分。
(3)利用subplot命令将y,y2和y3以子图的形式显示。
3.
实验结果及分析
1.(1)
t1 =
0.0780
实验结果及分析
(2)
t2 =
0.0780
实验结果及分析
(3)
t3 =
0
实验结果及分析
2.
实验
名称
实验四Matlab在通信系统中的应用
实验目的
1、练习通信原理中各种调制方法及简单通信系统的Matlab编程实现;
2、学习Matlab与Simulink的混合编程。
(2)y变成虚数后,重复上述运算。
(3)创建矩阵 ,并进行如下操作,取出矩阵A的前两行元素,生成一个新矩阵B。生成一个3*3的单位矩阵B,与矩阵A进行加、减、乘、除、相等、不相等运算。
(4)随机生成2个复数x1和x2,计算x3=x1÷x2,x3的实部real_x3、虚部image_x3、模abs_x3、幅角angle_x3,并把x1,x2,x3,real_x3,image_x3,abs_x3,angle_x3保存至文件complex_学号.mat。利用save命令保存数据为mat格式文件后,用clear命令清除工作空间中的变量,再利用load命令从文件读入数据。
2、流程控制语句相关操作。
(1)使用for语句及while语句求解1至100整数中奇数的和。
(2)求连续自然数的和,当和大于等于1000时,求最后一个自然数以及自然数的和。
matlab实验报告(实验1)
开课学院及实验室:机电学院 2012年11 月16日
(2)在E盘新建一个文件夹并把其加入搜索路径。
截屏相关过程结果并保存。
(3)逐步按照以下步骤进行并记录各个步骤的结果:
1)在指令窗中输入clear指令;
5)再次在指令窗中输入clear指令后,输入
8)把cow.m暂时移动到另一个无关的文件夹,如“我的文档”。
(2)保存全部变量为数据文件Mydata.mat;在对应文件夹中找到该文件,截屏后保存结果;
(4)删除全部内存变量;
Editor/Debugger和脚本编写初步
把以下指令编写为M脚本文件,并运行:
t=0:pi/50:4*pi;
y0=exp(-t/3); y=exp(-t/3).*sin(3*t); plot(t,y,'-r',t,y0,':b',t,-y0,':b') grid
什么?此方法也是查询函数注释的有效方法,请牢记。
matlab实验总结
matlab实验总结MATLAB是一种功能强大的数字计算软件,可用于科学计算、数据分析、可视化和算法开发等。
在实验过程中,使用MATLAB进行数据分析和可视化是一个常用的方法,有助于更好地理解实验结果和发现数据中的模式和趋势。
本文将介绍如何使用MATLAB进行实验总结,包括实验设计、数据处理和结果分析等方面。
我们将介绍一些常用的MATLAB工具和函数,并探讨如何在实验结束后进行数据管理和处理,以确保数据的完整性和准确性。
1. 实验设计在进行MATLAB实验总结之前,需要设计实验并确定实验条件。
实验设计需要考虑以下几个方面:- 实验目的和问题:明确实验的目的和问题,以便确定实验的参数和指标。
- 实验条件:确定实验所需的硬件和软件环境,包括计算机型号、MATLAB版本和编程语言等。
- 实验数据收集:确定如何收集实验数据,包括数据采集方式、数据存储和管理等。
- 实验数据处理:确定如何处理实验数据,包括数据清洗、特征提取和数据分析等。
2. 数据处理在MATLAB实验总结中,数据处理是非常重要的一部分。
数据处理的目标是将实验数据转化为易于理解和分析的形式。
以下是一些常用的MATLAB工具和函数,用于数据处理:- 数据清洗:用于清理和转换数据,以确保数据的质量和准确性。
例如,去除缺失值、异常值和重复值等。
- 数据转换:用于将数据转换为MATLAB能够理解和处理的形式。
例如,将数据转换为矩阵、向量或数组等。
- 特征提取:用于提取数据的特征,以便更好地理解数据。
例如,使用特征值和特征向量等。
- 数据分析:用于对数据进行分析和可视化,以更好地理解数据。
例如,使用聚类分析、回归分析等。
3. 结果分析在MATLAB实验总结中,结果分析是非常重要的一部分。
结果分析可以帮助更好地理解实验结果,发现实验数据中的模式和趋势。
以下是一些常用的MATLAB工具和函数,用于结果分析:- 可视化:用于将实验结果可视化,以便更好地理解实验结果。
(完整word)Matlab实验报告
实验一:Matlab操作环境熟悉一、实验目的1.初步了解Matlab操作环境.2.学习使用图形函数计算器命令funtool及其环境。
二、实验内容熟悉Matlab操作环境,认识命令窗口、内存工作区窗口、历史命令窗口;学会使用format 命令调整命令窗口的数据显示格式;学会使用变量和矩阵的输入,并进行简单的计算;学会使用who和whos命令查看内存变量信息;学会使用图形函数计算器funtool,并进行下列计算:1.单函数运算操作。
求下列函数的符号导数(1)y=sin(x);(2) y=(1+x)^3*(2-x);求下列函数的符号积分(1)y=cos(x);(2)y=1/(1+x^2);(3)y=1/sqrt(1—x^2);(4)y=(x1)/(x+1)/(x+2)求反函数(1)y=(x-1)/(2*x+3); (2) y=exp(x);(3) y=log(x+sqrt(1+x^2));代数式的化简(1)(x+1)*(x-1)*(x-2)/(x-3)/(x—4);(2)sin(x)^2+cos(x)^2;(3)x+sin(x)+2*x—3*cos(x)+4*x*sin(x);2.函数与参数的运算操作。
从y=x^2通过参数的选择去观察下列函数的图形变化(1)y1=(x+1)^2(2)y2=(x+2)^2(3) y3=2*x^2 (4) y4=x^2+2 (5) y5=x^4 (6) y6=x^2/2 3.两个函数之间的操作求和(1)sin(x)+cos(x) (2) 1+x+x^2+x^3+x^4+x^5乘积(1)exp(—x)*sin(x) (2) sin(x)*x商(1)sin(x)/cos(x); (2) x/(1+x^2); (3) 1/(x—1)/(x—2); 求复合函数(1)y=exp(u) u=sin(x) (2) y=sqrt(u) u=1+exp(x^2)(3) y=sin(u) u=asin(x) (4) y=sinh(u) u=-x实验二:MATLAB基本操作与用法一、实验目的1.掌握用MATLAB命令窗口进行简单数学运算。
MATLAB程序设计与应用(第二版)课后实验答案
Matlab课后实验题答案实验一 MATLAB运算基础1. 先求下列表达式的值,然后显示MATLAB工作空间的使用情况并保存全部变量。
(1)0 122sin851ze =+(2)21ln(2z x=,其中2120.455ix+⎡⎤=⎢⎥-⎣⎦(3)0.30.330.3sin(0.3)ln, 3.0, 2.9,,2.9,3.0 22a ae e az a a--+=++=--(4)2242011122123t tz t tt t t⎧≤<⎪=-≤<⎨⎪-+≤<⎩,其中t=0:0.5:2.52. 已知:1234413134787,2033657327A 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 和B123453166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(1) 求它们的乘积C 。
(2) 将矩阵C 的右下角3×2子矩阵赋给D 。
(3) 查看MATLAB 工作空间的使用情况。
4. 完成下列操作:(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量,删除其中的大写字母。
解:(1) 结果:(2).建立一个字符串向量 例如: ch='ABC123d4e56Fg9';则要求结果是:实验二 MATLAB 矩阵分析与处理1。
设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦,其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵,试通过数值计算验证22E R RS A OS +⎡⎤=⎢⎥⎣⎦。
MATLAB 实验二 基本操作
实验二 Matlab 基本操作(二)一 实验目的:1. 掌握矩阵方程的构造和运算方法2. 掌握基本Matlab 控制语句3. 学会使用Matlab 绘图二 实验内容1. 求解下列线性方程,并进行解的验证:⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----1323151122231592127x=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-0174 >> a=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13];b=[4;7;-1;0];x=a\b2、进行下列计算。
(1)k=∑=6322i i>>i=2:63;mysum=sum(2.^i)mysum =1.8447e+019(2)求出y=x*sin(x)在0<x<100条件下的每个峰值。
>>y='x.*sin(x)';fplot(y,[0 100]);min=fmin(y,0,100)min =54.99613、绘制下列图形。
(1)sin(1/t), -1<t<1;t=-1:0.02:1;y=sin(1./t);plot(t,y)(2)1-)7(cos 3t>> t=0:0.02:pi.*3;y=1-cos(7*t).^3;plot(t,y)4、已知系统闭环传递函数G (S ),分析系统稳定性及单位脉冲、单位阶跃响应。
22s 43206s 266)S (G s s s s s 23423+++++++=>> a=[1 3 4 2 2];b=[6 26 6 20];roots(a)ans =-1.4734 + 1.0256i-1.4734 - 1.0256i-0.0266 + 0.7873i-0.0266 - 0.7873i因为无右根,故系统稳定。
当单位脉冲输入时:>> [r p k]=residue(b,a);t=0:0.2:60;>> y1=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t)+r(3)*exp(p(3)*t)+r(4)*exp(p(4)*t); >> plot(t,y1)当输入单位阶跃函数时:>> a=[1 3 4 2 2 0];b=[6 26 6 20];[r p k]=residue(b,a);t=0:0.2:100;y2=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t)+r(3)*exp(p(3)*t)+r(4)*exp(p(4)*t)+r(5); plot(t,y2)。
MATLAB第二章实验报告
MATLAB数据及基本操作班级:11电信一班姓名:何得中学号:20111060108实验目的:1.熟悉MATLAB的组体成体系及工作环境;2.了解MATLAB的帮助系统;3.掌握MATLAB的数据类型及变量的基本操作;4.掌握矩阵建立的基本方法及其运算;5.掌握矩阵索引与分析方法;6.熟悉字符串单元及结构数据的建立及基本操作。
实验仪器:MA TLAB软件,电脑实验数据:>> A=[1 2 3;4 5 6;7 8 9]A = 1 2 34 5 67 8 9>> A=[1 2 34 5 67 8 9]A = 1 2 34 5 67 8 9>> x=0:pi/8:2*pi;>> y=cos(x)y = Columns 1 through 101.0000 0.9239 0.7071 0.3827 0.0000 -0.3827 -0.7071 -0.9239 -1.0000 -0.9239Columns 11 through 17-0.7071 -0.3827 -0.0000 0.3827 0.7071 0.9239 1.0000>> B=ones(3,5)B = 1 1 1 1 11 1 1 1 11 1 1 1 1>> zeros(3)ans = 0 0 00 0 00 0 0>> C=rand(3,4)C = 0.9501 0.4860 0.4565 0.44470.2311 0.8913 0.0185 0.61540.6068 0.7621 0.8214 0.7919>> D=randn(2,3)D = -0.4326 0.1253 -1.1465-1.6656 0.2877 1.1909>> eye(3)ans = 1 0 00 1 00 0 1>> M=magic(4)M = 16 2 3 135 11 10 89 7 6 124 14 15 1>> T=toeplitz(1:5,1:6)T = 1 2 3 4 5 62 1 234 53 2 1 2 3 44 3 2 1 2 35 4 3 2 1 2>> H=hilb(5)H = 1.0000 0.5000 0.33330.2500 0.20000.5000 0.3333 0.2500 0.2000 0.16670.3333 0.2500 0.2000 0.1667 0.14290.2500 0.2000 0.1667 0.1429 0.12500.2000 0.1667 0.1429 0.1250 0.1111>> a=3:18a =3 4 5 6 7 8 9 10 1112 13 14 15 16 17 18>> b=1:0.2:6b =Columns 1 through 91.0000 1.2000 1.4000 1.600 1.80002.0000 2.2000 2.4000 2.6000 Columns 10 through 182.80003.0000 3.2000 3.40003.6000 3.80004.0000 4.2000 4.4000 Columns 19 through 264.6000 4.80005.0000 5.20005.4000 5.6000 5.80006.0000>> linspace(1,15,10)ans =Columns 1 through 91.00002.5556 4.1111 5.6667 7.2222 8.7778 10.3333 11.8889 13.4444Column 1015.0000>> y=[0.9 0.3 0.4 0.1 0.5 0.6 0.2 ]; >> [sorted,index]=sort(y)sorted =0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.9000 index =4 7 2 3 5 6 1>> A=rand(4,5)A = 0.9218 0.9355 0.0579 0.13890.27220.7382 0.9169 0.3529 0.2028 0.19880.1763 0.4103 0.8132 0.1987 0.01530.4057 0.8936 0.0099 0.6038 0.7468>> [a,b]=min(A)a =0.1763 0.4103 0.0099 0.13890.0153b = 3 3 4 1 3>> min(A(:))ans = 0.0099>> V=[1 2 3 4]V = 1 2 3 4>> norm(V,2)ans =5.4772>> norm(V,1)ans = 10>> norm(V,inf)ans = 4>> V=[2 5 6 8 4 9]V = 2 5 6 8 4 9>> min(V)ans = 2>> V=[4 9 6;7 8 9;4 5 6]V = 4 9 67 8 94 5 6>> min(V)ans = 4 5 6>> max(V)ans = 7 9 9>> mean(V)ans = 5.0000 7.3333 7.0000>> median(V)ans = 4 8 6>> std(V)ans = 1.7321 2.0817 1.7321>> diff(V)ans =3 -1 3-3 -3 -3>> sort(V)ans =4 5 64 8 67 9 9>> length(V)ans = 3>> norm(V)ans =19.9125>> sum(V)ans = 15 22 21>> prod(V)ans =112 360 324>> cumsum(V)ans = 4 9 611 17 1515 22 21 >> cumprod(V)ans = 4 9 628 72 54112 360 324>> V=[1 2 8]V = 1 2 8>> U=[3; 4;7]U = 347>> dot(V,U)ans =67>> cross(V,U)ans =-18 17 -2>> A=[5 6 4;7 8 9;4 5 6]A = 5 6 47 8 94 5 6>> B=[6 2 7 ;7 5 8;4 6 8]B = 6 2 77 5 84 6 8>> C=A-BC = -1 4 -30 3 10 -1 -2>> C=A*BC = 88 64 115134 108 18583 69 116>> C1=A\BC1 =0.6667 -6.4444 -7.66670.6667 4.8889 7.3333-0.3333 1.2222 0.3333>> C2=B/AC2 = -3.0000 8.3333 -9.3333-1.6667 5.5556 -5.88890.0000 -1.3333 3.3333>> A=[1 2 3;4 5 6;7 8 9]A = 1 2 34 5 67 8 9>> D=A^2D = 30 36 4266 81 96102 126 150>> E=A^0.1E = 0.8466 + 0.2270i 0.3599 + 0.0579i -0.0967 - 0.1015i0.4015 + 0.0216i 0.4525 + 0.0133i 0.4432 - 0.0146i-0.0134 - 0.1740i 0.4848 - 0.0509i 1.0132 + 0.0820i>> A=[1+2*i 3;4 2+i]A =1.0000 + 2.0000i 3.00004.0000 2.0000 + 1.0000i>> A'ans = 1.0000 - 2.0000i 4.00003.0000 2.0000 - 1.0000i>> A=[1 2 3;4 5 6;7 8 9]A = 1 2 34 5 67 8 9>> B=[3 4 6;7 9 8;5 2 01]B = 3 4 67 9 85 2 1>> C=A.*BC = 3 8 1828 45 4835 16 9>> D=A./BD = 0.3333 0.5000 0.50000.5714 0.5556 0.75001.4000 4.0000 9.0000>> E=A.\BE = 3.0000 2.0000 2.00001.7500 1.8000 1.33330.7143 0.2500 0.1111>> A=[1+2*i 3;4 2+i]A = 1.0000 + 2.0000i 3.00004.0000 2.0000 + 1.0000i>> A.^2ans = -3.0000 + 4.0000i 9.000016.0000 3.0000 + 4.0000i>> A.'ans = 1.0000 + 2.0000i 4.00003.0000 2.0000 + 1.0000i>> x=0.1:0.1:0.8x =0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000>> y=sin(2*x).*cos(2*x)y =0.1947 0.3587 0.4660 0.4998 0.4546 0.3377 0.1675 -0.0292>> y=sin(x)y =0.0998 0.1987 0.2955 0.3894 0.4794 0.5646 0.6442 0.7174>> y=cos(x)y =0.9950 0.9801 0.9553 0.9211 0.8776 0.8253 0.7648 0.6967>> y=tan(x)y =0.1003 0.2027 0.3093 0.4228 0.5463 0.6841 0.8423 1.0296>> y=cot(x)y =9.9666 4.9332 3.2327 2.3652 1.8305 1.4617 1.1872 0.9712>> y=sec(x)y = 1.0050 1.0203 1.0468 1.0857 1.1395 1.2116 1.3075 1.4353>> y=csc(x)y = 10.0167 5.0335 3.3839 2.5679 2.0858 1.7710 1.5523 1.3940>> y=asin(x)y = 0.1002 0.2014 0.3047 0.4115 0.5236 0.6435 0.7754 0.9273>> y=acos(x)y = 1.4706 1.3694 1.2661 1.1593 1.0472 0.9273 0.7954 0.6435>> y=atan(x)y =0.0997 0.1974 0.2915 0.3805 0.4636 0.5404 0.6107 0.6747>> y=acot(x)y = 1.4711 1.3734 1.2793 1.1903 1.1071 1.0304 0.9601 0.8961>> y=asec(x)y =Columns 1 through 60 + 2.9932i 0 + 2.2924i0 + 1.8738i 0 + 1.5668i 0 +1.3170i 0 + 1.0986iColumns 7 through 80 + 0.8956i 0 + 0.6931i>> y=sinh(x)y =0.1002 0.2013 0.3045 0.4108 0.5211 0.6367 0.7586 0.8881>> y=cosh(x)y =1.0050 1.0201 1.0453 1.0811 1.1276 1.1855 1.2552 1.3374>> y=tanh(x)y = 0.0997 0.1974 0.2913 0.3799 0.4621 0.5370 0.6044 0.6640>> y=coth(x)y = 10.0333 5.0665 3.4327 2.6319 2.1640 1.8620 1.6546 1.5059>> y=sech(x)y = 0.9950 0.9803 0.9566 0.9250 0.8868 0.8436 0.7967 0.7477>> y=csch(x)y = 9.9834 4.9668 3.2839 2.4346 1.9190 1.5707 1.3182 1.1260>> y=asinh(x)y = 0.0998 0.1987 0.2957 0.3900 0.4812 0.5688 0.6527 0.7327>> y=acosh(x)y =Columns 1 through 60 + 1.4706i 0 + 1.3694i0 + 1.2661i 0 + 1.1593i 0 +1.0472i 0 + 0.9273iColumns 7 through 80 + 0.7954i 0 + 0.6435i>> y=atanh(x)y = 0.1003 0.2027 0.3095 0.4236 0.5493 0.6931 0.8673 1.0986>> y=acoth(x)y = Columns 1 through 60.1003 + 1.5708i 0.2027 + 1.5708i 0.3095 + 1.5708i 0.4236 + 1.5708i 0.5493 + 1.5708i 0.6931 + 1.5708i Columns 7 through 80.8673 + 1.5708i 1.0986 + 1.5708i >> y=asech(x)y = 2.9932 2.2924 1.8738 1.5668 1.3170 1.0986 0.8956 0.6931>> y=acsch(x)y = 2.9982 2.3124 1.9189 1.6472 1.4436 1.2838 1.1545 1.0476>> y=log(x)y =-2.3026 -1.6094 -1.2040 -0.9163 -0.6931 -0.5108 -0.3567 -0.2231>> y=log10(x)y = -1.0000 -0.6990 -0.5229 -0.3979 -0.3010 -0.2218 -0.1549 -0.0969>> y=log2(x)y =-3.3219 -2.3219 -1.7370 -1.3219 -1.0000 -0.7370 -0.5146 -0.3219>> c=complex(a,b)c = 0.1763 + 3.0000i 0.4103 + 3.0000i0.0099 + 4.0000i 0.1389 + 1.0000i 0.0153 + 3.0000i>> zc=conj(x)zc =0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000>> b=imag(x)b = 0 0 0 0 0 00 0>> a=real(x)a = 0.1000 0.2000 0.3000 0.40000.5000 0.6000 0.7000 0.8000>> y=fix(x)y = 0 0 0 0 0 0 0 0>> y=floor(x)y =0 0 0 0 0 0 0 0 >> y=ceil(x)y = 1 1 1 1 1 1 1 1>> y=round(x)y = 0 0 0 0 1 1 1 1>> y=sign(x)y = 1 1 1 1 1 1 1 1>> S=rats(x)S =1/10 1/5 3/10 2/5 1/2 3/5 7/10 4/5>> [N,D]=rat(x)N = 1 1 3 2 1 37 4D = 10 5 10 5 2 510 5>> A=[2 4 6;3 5 7;8 10 9]A = 2 4 63 5 78 10 9>> B=[4 1 5;7 5 0;9 3 6]B = 4 1 57 5 09 3 6>> A>Bans = 0 1 10 0 10 1 1>> A<Bans = 1 0 01 0 01 0 0>> A>=Bans = 0 1 10 1 10 1 1>> A<=Bans = 1 0 01 1 01 0 0>> A==Bans = 0 0 00 1 00 0 0>> A~=Bans =1 1 11 0 11 1 1>> A&Bans = 1 1 11 1 01 1 1>> A|Bans = 1 1 11 1 11 1 1>> ~Aans = 0 0 00 0 00 0 0>> ~Bans = 0 0 00 0 10 0 0>> x=0:pi/100:3*pi;>> y=sin(x);>> plot(x,y);>> y1=(x<pi|x>2*pi).*y;>> figure,plot(x,y1);>> r=(x>pi/3&x<2*pi/3)|(x>7*pi/3&8*pi/3);>> rn=~r;>> y2=r*sin(pi/3)+rn.*y1;>> figure,plot(x,y2)>> x=[0 4 5 9 8 6 74 85 7 4 85 8 4 5 8]x =0 4 5 9 8 6 74 85 7 4 85 8 4 5 8>> all(x)ans =0>> any(x)ans = 1>> find(x)ans =2 3 4 5 6 7 8 9 10 11 12 13 14 15>> isempty(x)ans =>> isglobal(x)Warning: isglobal is obsolete and will be discontinued. Type "help isglobal" for more details.ans =0>> isinf(x)ans =0 0 0 0 0 0 0 0 0 0 0 0 0 0 0>> isnan(x)ans =0 0 0 0 0 0 0 0 0 0 0 0 0 0 0>> isfinite(x)ans =1 1 1 1 1 1 1 1 1 1 1 1 1 1 1>> issparse(x)ans = 0>> isstr(x)ans =0>> not(x)ans =1 0 0 0 0 0 0 0 0 0 0 0 0 0 0>> y=[59 46 78 2 3 98 1 2 3 57 9 5 66 8 900 5]y =59 46 78 2 3 98 1 2 3 57 9 5 66 8 900 5>> y=[59 46 78 2 3 98 1 2 3 57 9 5 66 8 900 5]y = 59 46 78 2 3 98 1 2 3 57 9 5 66 8 900 5>> x=[0 4 5 9 8 6 74 85 7 4 85 8 4 5 8]x = 0 4 5 9 8 6 74 85 7 4 85 8 4 5 8>> A=[4 -40 2; 23 50 17;18 29 14]A = 4 -40 223 50 1718 29 14>> [r,c]=find(A>=20&A<=30)r = 23c = 12>> bitand(12,6)ans = 4>> bitshift(12,1)ans = 24>> bitxor(12,1)ans =13>> bitset(12,1)ans = 13 >> bitget(12,1)ans =0>> A =[ 87 59 56 85 66 21 12 35 62 64 69 88 99 81]A =87 59 56 85 66 21 1235 62 64 69 88 99 81>> A(2)ans =59>> A([1,2,8])ans =87 59 35>> A([9,6,2])ans =62 21 59>> A([end-4:end])ans =64 69 88 99 81>> A([1:5,5:-1:1])ans = 87 59 56 85 66 66 85 56 59 87>> A(3)=-1A = 87 59 -1 85 66 2112 35 62 64 69 88 99 81>> A(15)=-8A =87 59 -1 85 66 2112 35 62 64 69 88 99 81 -8>> A(20)=-1A = Columns 1 through 1687 59 -1 85 66 21 12 35 62 64 69 88 99 81 -8 0Columns 17 through 200 0 0 -1>> A=rand(8)A = 0.9501 0.8214 0.9355 0.13890.4451 0.8381 0.3046 0.37840.2311 0.4447 0.9169 0.2028 0.9318 0.0196 0.1897 0.86000.6068 0.6154 0.4103 0.1987 0.4660 0.6813 0.1934 0.85370.4860 0.7919 0.8936 0.6038 0.4186 0.3795 0.6822 0.59360.8913 0.9218 0.0579 0.2722 0.8462 0.8318 0.3028 0.49660.7621 0.7382 0.3529 0.1988 0.5252 0.5028 0.5417 0.89980.4565 0.1763 0.8132 0.0153 0.2026 0.7095 0.1509 0.82160.0185 0.4057 0.0099 0.7468 0.6721 0.4289 0.6979 0.6449>> A(3,3)ans =0.4103>> A(8)ans =0.0185>> sub2ind(size(A),3,3)ans =19>> [i,j]=ind2sub(size(A),19)i =3j =3>> A(:,4)ans =0.13890.20280.19870.60380.27220.19880.01530.7468>> A(3,:)ans =0.6068 0.6154 0.4103 0.1987 0.4660 0.6813 0.1934 0.8537>> A(end,:)ans =0.0185 0.4057 0.0099 0.7468 0.6721 0.4289 0.6979 0.6449>> A(3:4,5:6)ans =0.4660 0.68130.4186 0.3795>> A(2:2:4,1:2:5)ans =0.2311 0.9169 0.93180.4860 0.8936 0.4186>> A=reshape(A,5,5)A = 1 6 11 16 212 7 12 17 223 8 13 18 234 9 14 19 245 10 15 20 25>> A(:)'ans =Columns 1 through 161 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16Columns 17 through 2517 18 19 20 21 22 23 24 25>> A(:,1:2:5)=[]A = 6 167 178 189 1910 20>> A=[9 8 7;6 5 4;3 2 1];>> D=diag(A)D =951>> D1=diag(A,1)D1 =84>> D2=diag(A,-1)D2 =62>> V=[1 2 3 4];>> diag(V)ans = 1 0 0 00 2 0 00 0 3 00 0 0 4>> diag(1:3,-1)ans = 0 0 0 01 0 0 00 2 0 00 0 3 0>>A=[17,1,0,15;5,7,14,16;4,0,13,0;10,12,19,21]; >> D=diag(4:-1:1);>> D*Aans = 68 4 0 6015 21 42 488 0 26 010 12 19 21>> triu(A)ans =17 1 0 150 7 14 160 0 13 00 0 0 21>> triu(A,1)ans = 0 1 0 150 0 14 160 0 0 00 0 0 0>> B=A'B =17 5 4 101 7 0 120 14 13 1915 16 0 21>> rot90(A)ans =15 16 0 210 14 13 191 7 0 1217 5 4 10>> rot90(A,3)ans =10 4 5 1712 0 7 119 13 14 021 0 16 15>> B=fliplr(A)B = 15 0 1 1716 14 7 50 13 0 421 19 12 10>> B=flipud(A)B = 10 12 19 214 0 13 05 7 14 1617 1 0 15>> A=[4 5 6 9;7 8 9 5;4 5 6 9;7 8 9 5]; >> det(A)ans =0>> rank(A)ans =2>> cond(A)ans =1.0968e+017>> inv(A)Warning: Matrix is singular to working precision.ans = Inf Inf Inf InfInf Inf Inf InfInf Inf Inf InfInf Inf Inf Inf>> A=[4 8 9;8 2 7;6 3 8];>> inv(A)ans = 0.0568 0.4205 -0.43180.2500 0.2500 -0.5000-0.1364 -0.4091 0.6364>> [V,D]=eig(A)V = -0.6438 -0.7444 -0.5061-0.5467 0.6392 -0.5515-0.5354 0.1931 0.6631D = 18.2776 0 00 -5.2030 00 0 0.9254>> [U,S,D]=svd(A)U = -0.6378 0.7375 -0.2220-0.5426 -0.6348 -0.5500-0.5465 -0.2304 0.8051S = 18.9278 0 00 5.2884 00 0 0.8791D = -0.5374 -0.6639 -0.5200-0.4135 0.7449 -0.5236-0.7350 0.0663 0.6749>> [L,U]=lu(A)L = 0.5000 1.0000 01.0000 0 00.7500 0.2143 1.0000U = 8.0000 2.0000 7.00000 7.0000 5.50000 0 1.5714>> [Q,R]=qr(A)Q = -0.3714 0.9114 -0.1774-0.7428 -0.4063 -0.5322-0.5571 -0.0659 0.8278 R =-10.7703 -6.1279 -12.99870 6.2808 4.83140 0 1.3009>> a=1234a = 1234>> class(a)ans =double>> size(a)ans = 1 1>> b='1234'b =1234>> class(b)ans =char>> size(b)ans = 1 4>> 'I am astudent'ans =I am astudent>> a='This is No.2.23 Example!'a =This is No.2.23 Example!>> b=a(1:7)b =This is>> c=a(12:end)c =2.23 Example!>> a='Good';>> b='Noon';>> length(a)==length(b)ans = 1>> c=[a,'',b]c =GoodNoon>> d=[a;b]d =GoodNoon>> size(c)ans = 1 8>> size(d)ans = 2 4>> a='Good Noon';>> b=double(a)b = 71 111 111 100 32 78 111 111 110>> c='再见!'c =再见!>> c='再见!'>> d=double(c)d = 20877 35265 65281>> char(d)ans =再见!>> a='The first string';>> b='The second string';>> c=strcmp(a,b)c =0>> d=strncmp(a,b,4)d =1>> X='A friend in need is a friend indeed'; >> Y='friend';>> a=findstr(Y,X)a = 3 23>> b=strfind(Y,X)b = []>> S=['1 2 3';'2 3 4'];>> A=str2num(S)A = 1 2 32 3 4>> B=str2num('6-8i')B = 6.0000 - 8.0000i>> D=num2str(rand(2,3),6)D =0.950129 0.606843 0.8912990.231139 0.485982 0.762097>> d=189;>> h=dec2hex(d)h =BD>> c=dec2base(d,7)c =360>> b=dec2bin(d)b =10111101>> bin2dec(b)ans = 189>> A={ones(3,3,3),'Welcome';30.34,1:200} A = [3x3x3 double] 'Welcome'[ 30.3400] [1x200 double]>>B=[{ones(3,3,3)},{'Welcome'};{30.34},{1:20 0}]B = [3x3x3 double] 'Welcome'[ 30.3400] [1x200 double]>> C={5}C = [5]>> C(2,3)={7}C = [5] [] [][] [] [7]>> isequal(A,B)ans =1>> A={ones(3,3,3),'Welcome';30.34,1:200}; >> b=A(1,2)b = 'Welcome'>> class(B)ans =cell>> C=A{1,2}C =Welcome>> class Cans =char>> D=A{1,2}(6)D =m>> E=A{2,2}([end:-1:190])E =200 199 198 197 196 195 194 193 192 191 190>> class(E)ans =double>> N=A{3}([1 3 5 7])N =Wloe>> A={ones(3,3,3),'Welcome';30.34,1:200}; >> B=cell(2);>> B(:,1)={char('Good','Morning');1:10}B = [2x7 char ] [][1x10 double] []>> C=[A,B]C = [3x3x3 double] 'Welcome' [2x7 char ] [][ 30.3400] [1x200 double] [1x10 double] []>> D=[A,B;C]D = [3x3x3 double] 'Welcome' [2x7 char ] [][ 30.3400] [1x200 double] [1x10 double] [][3x3x3 double] 'Welcome' [2x7 char ] [][ 30.3400] [1x200 double] [1x10 double] []>> D(4,:)=[]D = [3x3x3 double] 'Welcome' [2x7 char ] [][ 30.3400] [1x200 double] [1x10 double] [][3x3x3 double] 'Welcome' [2x7 char ] []>> E=reshape(D,2,3,2)E(:,:,1) = [3x3x3 double] [3x3x3 double] [1x200 double][ 30.3400] 'Welcome' 'Welcome'E(:,:,2) = [2x7 char ] [2x7 char] [][1x10 double] [] []>>A={randn(3,3,2),'Good',pi;29,4+7*i,zeros(4)} A = [3x3x2 double] 'Good' [ 3.1416][ 29] [4.0000+ 7.0000i] [4x4 double]>> B=cellfun('isreal',A)B = 1 1 11 0 1>> C=cellfun('length',A)C = 3 4 11 1 4>> Member.code='09021';>> ='Liu';>> Member.age=22;>> Member.grade=uint16(3);>> MemberMember =code: '09021'name: 'Liu'age: 22grade: 3>> Member(4).name='Wang';>> Member(4).grade=2;>> Member(2)ans = code: []name: []age: []grade: []>>Memebr=struct('code','09021','name','Liu','a ge',22,'grade',uint16(3))Memebr = code: '09021'name: 'Liu'age: 22grade: 3>>Member=struct('code',{'09021','09034'},'name' ,{'Liu','Wamg'},'age',{22,24},'grade',{2,3}) Member = 1x2 struct array with fields:codenameagegrade>>Member=struct('code',{},'name',{},'age',{},'gr ade',{})Member = 0x0 struct array with fields:codenameagegrade>>Member=struct('code',{'09021','09034'},'na me',{'Liu','Wang'},'age',{22,24},'grade',{2,3},' score',{[78 89;90 68],[91 76;89 97]}) Member = 1x2 struct array with fields:codenameagegradescore>> Member(1).scoreans =78 8990 68>> Member(2).score(2,:)ans =89 97>> Member.codeans =09021ans =09034>> Member.('name')ans =Liuans =Wang>> ='Liu';M.ID=1;>> M(2,2).name='Wang';M(2,2).ID=2; >> M2=setfield(M,{2,1},'name','Zhang'); >> ans =Liuans =[]ans = []ans =Wang>> ans =Liuans =Zhangans = []ans =Wang>> fieldnames(M)ans = 'name''ID'>> M3=orderfields(M)M3 = 2x2 struct array with fields:IDname。
Matlab实验报告_2
实验一 Matlab基础知识一、实验目的:1.熟悉启动和退出Matlab的方法。
2.熟悉Matlab命令窗口的组成。
3.掌握建立矩阵的方法。
4.掌握Matlab各种表达式的书写规则以及常用函数的使用。
二、实验内容:1.求[100,999]之间能被21整除的数的个数。
(rem)2.建立一个字符串向量,删除其中的大写字母。
(find)3.输入矩阵,并找出其中大于或等于5的元素。
(find)4.不采用循环的形式求出和式6312ii=∑的数值解。
(sum)三、实验步骤:●求[100,199]之间能被21整除的数的个数。
(rem)1.开始→程序→Matlab2.输入命令:»m=100:999;»p=rem(m,21);»q=sum(p==0)ans=43●建立一个字符串向量,删除其中的大写字母。
(find)1.输入命令:»k=input('’,’s’);Eie48458DHUEI4778»f=find(k>=’A’&k<=’Z’);f=9 10 11 12 13»k(f)=[ ]K=eie484584778●输入矩阵,并找出其中大于或等于5的元素。
(find)1.输入命令:»h=[4 8 10;3 6 9; 5 7 3];»[i,j]=find(h>=5)i=3 j=11 22 23 21 32 3●不采用循环的形式求出和式的数值解。
(sum)1.输入命令:»w=1:63;»q=sum(2.^w)q=1.8447e+019实验二 Matlab 基本程序一、 实验目的:1. 熟悉Matlab 的环境与工作空间。
2. 熟悉M 文件与M 函数的编写与应用。
3. 熟悉Matlab 的控制语句。
4. 掌握if,switch,for 等语句的使用。
二、 实验内容:1. 根据y=1+1/3+1/5+……+1/(2n-1),编程求:y<5时最大n 值以及对应的y 值。
matlab实验总结
matlab实验总结Matlab实验总结近年来,计算机科学与技术领域的发展迅猛,各种编程语言和工具层出不穷。
而在这些编程语言中,Matlab作为一种高效、灵活且易于学习的工具,被广泛应用于科学计算和工程领域。
本文将对我在Matlab实验中的经验和总结进行分享。
在进行Matlab实验时,我们需要充分了解Matlab的基本语法和常用函数。
Matlab拥有丰富的函数库,可以实现各种复杂的计算和数据处理。
因此,对于初学者来说,学会查找和使用Matlab函数是至关重要的。
通过查阅Matlab官方文档、教材和网络资源,我们可以快速找到我们所需要的函数,并了解它们的使用方法和参数设置。
编写Matlab代码时,我们需要注重代码的可读性和可维护性。
良好的代码结构和注释可以使我们的代码更易于理解和修改。
在编写代码时,我们应该尽量使用有意义的变量名和函数名,避免使用过于简单或过于复杂的命名方式。
此外,我们还可以使用Matlab提供的代码自动缩进和格式化功能,以保持代码的整洁和统一。
在实验过程中,我们还需要注意数据的处理和分析。
Matlab提供了丰富的数据处理和分析函数,可以帮助我们对实验数据进行统计、绘图和模型拟合等操作。
在使用这些函数时,我们需要了解函数的输入参数和输出结果的含义,以确保我们得到准确的结果。
在进行Matlab实验时,我们还需要注意代码的效率和性能。
Matlab 是一种解释性语言,相比于编译型语言,它的执行效率较低。
因此,在处理大量数据或进行复杂计算时,我们应该尽量优化代码,减少不必要的计算和内存占用。
例如,可以使用矩阵运算代替循环,或者使用向量化的方法进行计算。
我们还需要注意Matlab实验的结果分析和可视化。
Matlab提供了丰富的绘图函数,可以帮助我们将实验结果以直观的方式展示出来。
在进行结果分析时,我们应该注重对实验数据的解释和推断,以及对实验结果的合理性和可靠性的评估。
Matlab作为一种强大的科学计算工具,可以帮助我们进行各种复杂的计算和数据处理。
matlab实验总结
matlab实验总结本次实验是一个MATLAB实验的总结,我们通过这次实验学习了MATLAB的一些基本功能和使用方法。
在实验中,我们探索了MATLAB的各种工具箱和函数,包括数据可视化、矩阵运算、符号计算等等。
通过这些实验,我们对MATLAB有了更深入的了解,并且对如何使用MATLAB解决实际问题有了更具体的认识。
首先,在本次实验中,我们学习了MATLAB的数据可视化功能。
通过使用plot函数,我们能够绘制出各种曲线图和散点图,以更直观地呈现数据之间的关系。
此外,我们还学习了如何设置图例、坐标轴和标签,以及如何调整图表的样式和颜色。
通过这些实验,我们掌握了如何将数据进行可视化,并且能够从图表中得到更直观的信息。
其次,我们还学习了MATLAB的矩阵运算功能。
MATLAB提供了一系列强大而方便的矩阵运算函数,包括矩阵乘法、矩阵转置、矩阵求逆等等。
通过这些函数,我们能够高效地进行复杂的矩阵运算,并且能够解决一些实际问题,例如线性方程组的求解和最小二乘法拟合等。
除了数据可视化和矩阵运算,我们还学习了MATLAB的符号计算功能。
通过使用符号计算工具箱,我们可以进行各种数学运算和代数操作,例如符号求导、符号积分、多项式运算等等。
这些功能可以帮助我们更好地理解和分析数学问题,解决复杂的数学方程和数值计算问题。
在实验过程中,我们还学习了如何使用MATLAB编辑器和命令窗口。
MATLAB编辑器是一个功能强大的编辑器,用于编写和调试MATLAB代码。
在编辑器中,我们可以方便地编写脚本和函数,并且可以进行代码的调试和查错。
命令窗口则是一个交互式的界面,我们可以直接在命令窗口中输入MATLAB代码,并且可以立即得到执行结果。
总的来说,这次MATLAB实验为我们提供了一个全面的学习机会,我们通过实际操作和实验练习,掌握了MATLAB的基本功能和使用方法。
通过数据可视化、矩阵运算和符号计算,我们能够更好地理解和解决实际问题。
MATLAB实验报告 (2)
广西大学实验报告纸姓名:徐律1002100405电气工程学院自动化102班2012年12月27日【实验名称】MATLAB的图形绘制【实验要求】1、学习MATLAB图形绘制的基本方法;2、熟悉和了解MATLAB图形绘制程序编辑的基本指令;3、熟悉掌握利用MATLAB图形编辑窗口编辑和修改图形界面,并添加图形的各种标注;4、掌握plot、subplot的指令格式和语法。
【实验基本知识】1、plot(x,y)------绘制由x,y所确定的曲线;2、多组变量绘图:plot(x1, y1, 选项1, x2, y2, 选项2, ……);3、双Y轴绘图:plotyy()函数;4、图形窗口的分割;5、图形编辑窗口的使用。
【实验内容】【1】二维曲线绘图基本指令演示。
本例运作后,再试验plot(t), plot(Y), plot(Y,t) ,以观察产生图形的不同。
t=(0:pi/50:2*pi)';k=0.4:0.1:1;Y=cos(t)*k;plot(t,Y)plot(t,Y); plot(t);plot(Y);plot(Y,t); %图形如下图所示【2】用图形表示连续调制波形Y=sin(t)sin(9t)及其包络线。
t=(0:pi/100:pi)';y1=sin(t)*[1,-1];y2=sin(t).*sin(9*t);t3=pi*(0:9)/9;y3=sin(t3).*sin(9*t3);plot(t,y1,'r:',t,y2,'b',t3,y3,'bo')axis([0,pi,-1,1])【3】在一个图形窗口绘制正弦和余弦曲线,要求给图形加标题“正弦和余弦曲线”,X轴Y轴分别标注为“时间t”和“正弦、余弦”,在图形的某个位置标注“sin(t)”“cos(t)”,并加图例,显示网格,坐标为正方形坐标系。
实验程序:实验结果:Sin(t)Cos(t)【4】绘制向量x=[1 3 0.5 2.5 2]的饼形图,并把3对应的部分分离出来。
MATLAB应用课题二GUI拨号实验
《MATLAB基础与应用》课题二:基于MATLAB的双音多频(DTMF)信号的产生、接收以及可视化拨号系统的实现内容一:研究DTMF信号的产生以及基于Goertzel算法的接收原理,并使用MATLAB 编程实现DTMF拨号系统的信号发生及接收。
双音多频系统,简称为DTMF(Double Tone Multi Frequency),是用按键进行电话拨号的一种制式,由美国AT﹠T贝尔实验室开发。
作为实现电话号码快速可靠传输的一种技术,它具有很强的抗干扰能力和较高的传输速度,可广泛应用于电话通信系统中。
但绝大部分是用作电话的音频拨号。
另外,它也可以在数据通信系统中广泛地用来实现各种数据流和语音等信息的远程传输。
每个DTMF信号是由2个单频率信号叠加而成的,其中一个频率处于低频段(小于1KHz),另一个频率处于高频段(1KHz~2KHz)。
高低频段各有四个频率,四个低频表示四行,四个高频表示四列,高、低频组合可以提供16个字符。
第四列的最高频率暂时不用,因此可以提供12个字符,包括0到9中一个特定的十进制数字,或者字符*和#。
而其它四个字符保留以备扩展使用。
字符分配如表1所示。
表1 双音多频通信系统编码表(第四列备用)DTMF是一个典型的小型数字信号处理系统。
它既有模拟信号的生成和传输部分,需要使用D/A转换;又有在接收端对它进行A/D转换并进行数字信号处理即检测的部分。
而且为了提高系统的检测速度并降低成本,还开发了一种特殊的DFT算法,称为Goertzel算法。
这种算法既可以用硬件(专用芯片)也可以用软件实现。
所以DTMF系统的设计问题是理论与工程相结合的一个典型范例。
1、DTMF信号的产生:每个字符对应的DTMF信号由两个正弦波的数字样本叠加而成。
正弦波的采样率为8kHz,硬件必需每125ms输出一个样本。
将这个叠加信号送到D/A转换器变换成模拟信号,再通过电话线传送到交换机。
2、DTMF信号的接收:在接收端,将接收到的模拟音频信号进行A/D转换,恢复为数字信号,然后检测其频谱以确定所发送的字符。
matlab实验2
实验一认识Matlab一、目的初识Matlab,掌握Matlab软件的安装,了解其操作界面的布局,运行环境,及个菜单、工具、窗口的作用。
二、内容1 安装Matlab软件;2 启动Matlab,认识其操作界面和运行环境。
3 在命令窗口输入一些命令如:距阵求和A=[1 2 3;4 5 6]B=[2 3 4;4 3 5]C=A+B绘图x=[0:0.5:360]*pi/180;plot(x,sin(x),x,cos(x));求解方程组p=[3,7,9,0,-23];x=roots(p)4 查看历史命令窗口、工作空间、当前目录窗口,各有什么变化。
三、思考题1 有哪些方法可以使得输出的数(默认保留4位小数)的位数改变。
2 在历史命令窗口,双击某一命令或变量,会有什么反应?3 通过帮助命令Help,如何获得函数如Sin(x)的功能信息?实验二矩阵的建立和基本运算一、目的掌握建立矩阵的基本方法及进行一些基本运算。
二、内容1 启动Matlab2 输入如下命令如:A=[1 2 3;4 5 6]B=[2 3 4;4 3 5]C1=A+BC2=A-BC3=A*BC4=A*B’b=[4 8 12]’三、实验题目M=[1 2 3;4 5 6;7 8 9],G=[8.57 -12.44 3.59;14.21 22.30 -16.16;7.42 -8.93 11.70] Inv(M) %求逆M’ %转置※以下请写出正确代码(1)从M中提取对角元素赋值给变量DM、再将这几个对角元素转化为对角矩阵赋值给变量AM;(2)将M中的第一列除以2、第二列除以3、第三列除以4后赋值给变量H;(3)查找M中大于5的元素的下标。
(4)生成一个5行5列的单位矩阵。
(5)查找M、G中最大元素与最小元素。
(6) 查找M、G中绝对值最大元素与绝对值最小元素。
(7) 查找M、G中绝对值最大元素与绝对值最小元素的下标。
(8)提取M的第三行赋值给变量MH3、提取M的第二列赋值给变量ML3。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
ezplot(‘f(x,y)’)
f ( x, y) 0 在区间 2 x 2 ,2 y 2 上的图形
ezplot(‘f’,’g’,[a,b])
x f (t ), y g (t ) 在区间 a t b 上的图形
例6 1、作图y=sinx、 y=sinx,x∈ [-, ]
例5
x=linspace(0,2*pi,30); y=sin(x); z=cos(x); plot(x,y,'r',x,z,'go');
title('曲线y1=sinx与y2=cosx');
xlabel('x轴');ylabel('y轴'); legend('y1=sin(x)','y2=cos(x)'); gtext('y1=sin(x)'); gtext('y2=cos(x)');
三、极坐标方程绘图: polar 如果想用利用曲线的极坐标方程作图,可使用 polar命令,其基本形式是: polar(theta,rho) 例7 曲线的极坐标方程为: 3 cos3 , 要作出它的图形,应输入: theta=0:0.1:2*pi; >> rho=3*cos(3*theta); >> polar(theta,rho)
(6)给图形加注释 为了进一步使图形具有说服力和可读性,我们还 经常给图形添加注释,例如图形的名称、坐标轴的名 称、图例及文字说明,针对本例图形,给图形取名为 “first figure” ,x和y坐标轴分别取名“横”、“纵”, 图例设置为“ y=cos(x)-sin(x) ”,并在(pi/2,-1)点加上 文字说明“ final point”,输入以下代码:
y cos3 t 的图形。 2、 在[0, 2 ]上画出星形线 3 x sin t
x 在 [ 2,0.5],[0, 2] 上画出隐函数 e +sin(xy)=0的图形。 3、
subplot(2,2,1);ezplot(‘sin(x)’); subplot(2,2,2);ezplot(‘sin(x)’,[-pi,pi]); subplot(2,2,3); ezplot('cos(t)^3','sin(t)^3',[0,2*pi]) subplot(2,2,4); ezplot('exp(x)+sin(x*y)',[-2,0.5,0,2])
• 2.了解常见曲线的图形特征,会画复合函数、参 数函数、隐函数及极坐标下函数的图形 • 3.学会设置图形的属性 • 【实验要求】 • 熟悉MATLAB绘图命令plot、ezplot、polar、plot3 、mesh等
一、二维曲线绘图plot
如何画出 y=sin(x) 在 [0, 2*pi] 上的图像?
点和线的基本属性
plot(x,y,string) 其中 string 是用 单引号 括起来的字符串,用来 指定图形的属性(点、线的形状和颜色)
>> x=[0:0.2:2*pi]; >> plot(x,cos(x)); >> plot(x,cos(x),’r+:’); >> plot(x,cos(x),’bd-.’); >> plot(x,cos(x),’k*-’);
蓝色、点划线 离散点为菱形
红色、虚线、 离散点用加号
属性可以全部指定,也 可以只指定其中某几个 排列顺序任意
黑色、实线 离散点用星号
图形的基本属性汇总
线型
: -. -实线 点线 点划线 虚线
点标记
. o x + * s d ^ v > < p h 实心点 空心圆 叉号 加号 星号 方形 菱形 向上三角 向下三角 向右三角 向左三角 五角星 六角星
x 是向量, y 是矩阵, ①若 x 的长度与 y 的行数相等,则将 x 与 y 中的各列相对应,绘制多条平面曲线; ②若 x 的长度与 y 的列数相等,则将 x 与 y 中的各行相对应,绘制多条平面曲线。 此时 x 的长度必须等于 y 的行数或列数。
练习2: 指出以下各个绘图命令的输出图形分别是什么, 并上机验证
出函数图形 例:>> x=[0:pi/10:2*pi]; >> y=sin(x); >> plot(x,y,'bo-')
在MATLAB中用图形函数绘图的一般操作步 骤分为7步,如表所示,下面以绘制一个简单三 角函数的图形为例,详细介绍各个步骤。
例1
用图形表示
y cos x sin x ,
>>title(‘first figure’); >> xlabel(‘横’);ylabel(‘纵’); >> legend(‘y=cos(x)-sin(x)’); >> text(pi/2,-1, ‘final point’);
(7)图形输出 大家完成图形的绘制和编辑以后,需要将图 形打印或是存放在其他文档中,最简单的方法是使用 【file】菜单中的【save as】选项,并且大家可以自 行选择需要保存的格式类型。
• 手工作图
找点: x=0, pi/3, pi/2, 2*pi/3, pi, … 计算函数值: y=sin(0), sin(pi/3), sin(pi/2), … 描点:在坐标系中画出这些离散点
用直线或曲线连接这些点,得到函数的大致图形
• MATLAB 作图
给出离散点列: x=[0:pi/10:2*pi] 计算函数值: y=sin(x) 画图: MATLAB 二维绘图命令plot(x,y,string)作
针对这种情况,MATLAB提供了若干特殊图形绘 制函数。接下来主要介绍特殊图形的绘制方法,主 要图形包括:条形图、区域图、饼状图、柱状图、 离散图、罗盘图、羽毛图、轮廓图以及箭头图等。
特殊二维绘图函数
bar –––– 绘制直方图 hist –––– 绘制统计直方图 stairs –––– 绘制阶梯图 stem –––– 绘制火柴杆图 rose –––– 绘制统计扇形图 comet –––– 绘制彗星曲线 area –––– 区域图 pie –––– 饼图
颜色
y m c r g b w k 黄色 红紫色 青色 红色 绿色 蓝色 白色 黑色
例4 在 [0,2 ] 用红线画sin(x),用绿圈画cos(x)。 >> x=linspace(0,2*pi,30); >>y=sin(x); >>z=cos(x); >>plot(x,y,'r',x,z,'go'); x=0:pi/10:2*pi;
基本形式
plot(x,y)
x, y 都是向量,则以x中元素为横坐标,y中 元素为纵坐标作平面曲线。此时 x, y 必须具 有相同长度(维数)。
例2:y = cos(x)在[0, 4π]上的图像.
>> x=[0:0.1:4*pi]; >> y=cos(x); >> plot(x,y);
可设置图名、标注等
>> t=[0 1]; x=[1 2]; y=[x;3 4]; z=[y;5 6]; >> plot(t,x); >> plot(t,y); >> plot(t,y'); >> plot(t,z); >> plot(t,z');
具有两个纵坐标标度的图形
在MATLAB中,如果需要绘制出具有不同纵坐 标标度的两个图形,可以使用plotyy绘图函数。调 用格式为:
二、绘制符号函数二维曲线图: ezplot ezplot即Easy to use function plotter,是一个易用的 一元函数绘图函数 。特别是在绘制参数函数和隐函 数的图像时,ezplot要比plot更方便。因为plot绘制图 形时要指定自变量的范围,而ezplot无需数据准备, 直接绘出图形。
自己动手
• 练习1:在0≤x≤2区间内,绘制曲线 •
y=2e-0.5x cos(4πx)
解答 excise1.m
基本形式
plot(x,y)
x, y 都是向量,则以 x 中元素为横坐标, y 中元素为纵坐标作平面曲线。 此时 x, y 必须具有相同长度。 当x ,y 是同维矩阵时,则以x,y对应列元素为 横、纵坐标分别绘制曲线,曲线条数等于矩阵 的列数。 此时 x, y 必须具有相同的大小。
数学实验
(二)
——MATLAB 图形功能与分形
课堂及考试要求
• 课程安排:共8次课,16课时。 • 成绩分配:满分100分 –平时成绩50分。包括:
上课出勤:20分,点名缺勤3次取消考试资格。
平时作业:30分,实验结果每次上交实验报告。
–考试成绩50分。
实验一 MATLAB绘图
• 【实验目的】
• 1. 掌握利用MATLAB绘制图形的基本方法
x [
, ] 2 2
(1)准备图形数据 图形是数据变换规律和相互联系的形象显示,所以 数据输入是图形创建的基础,我们需要选定数据的 范围,选择对应范围的自变量,计算相应的函数值。 本例中,我们需在命令窗口输入如下代码: >>x = -pi/2:.1:pi/2; >>y=cos(x)-sin(x);
>>plot(x,y,‘k-.*’)
(5)设置坐标范围及网格线 完成基本绘制后,在图形概貌的基础上对图形进 行编辑,增强图形的可读性;一般说来,大家需要设 置坐标轴的刻度范围和分隔线。本例系统默认的坐标 轴范围是x,y∈[-2,2],且没有网格线,将本例的x坐标轴 和y坐标轴的范围分别设置为[-π,π],和[-4,4],且给坐 标加上分隔线,可以输入以下代码: >>axis([-pi,pi,-3,3]) grid on