MATLAB实验(1)

合集下载

MATLAB实验答案(桂电)

MATLAB实验答案(桂电)

实验一 MATLAB入门(1)1.实验目的:(1)了解MATLAB的体系结构与特点,熟悉其集成开发环境。

(2)熟悉MATLAB界面窗口的功能和使用方法。

(3)熟悉MATLAB的帮助系统及使用方法。

(4)了解MATLAB的的数据类型、基本形式和数组的产生方法。

(5)掌握MATLAB基本的数学运算操作。

2.实验原理(1)MATLAB简介MATLAB是美国MathWorks公司开发的高性能的科学与工程计算软件。

它在数值计算、自动控制、信号处理、神经网络、优化计算、小波分析、图像处理等领域有着广泛的用途。

近年来, MATLAB在国内高等院校、科研院所的应用逐渐普及,成为广大科研、工程技术人员必备的工具之一。

MATLAB具有矩阵和数组运算方便、编程效率极高、易学易用、可扩充性强和移植性好等优点,俗称为“草稿纸式的科学计算语言”。

它把工程技术人员从繁琐的程序代码编写工作中解放出来,可以快速地验证自己的模型和算法。

经过几十年的扩充和完善,MATLAB已经发展成为集科学计算、可视化和编程于一体的高性能的科学计算语言和软件开发环境,整套软件由MATLAB开发环境、MATLAB语言、MATLAB数学函数库、MATLAB图形处理系统和MATLAB应用程序接口(API)等五大部分组成。

MATLAB的主要特点包括强大的计算能力(尤其是矩阵计算能力)、方便的绘图功能及仿真能力、极高的编程效率。

另外,MATLAB还附带了大量的专用工具箱,用于解决各种特定领域的问题。

通过学习软件的基本操作及其编程方法,体会和逐步掌握它在矩阵运算、信号处理等方面的功能及其具体应用。

通过本课程实验的学习,要求学生初步掌握MATLAB的使用方法,初步掌握M文件的编写和运行方法,初步将MATLAB运用于数字信号处理中。

循序渐进地培养学生运用所学知识分析和解决问题的能力。

(2)MATLAB的工作界面(Desktop)与操作MATLAB 安装成功后,第一次启动时,主界面如下图(不同版本可能有差异)所示:其中① 是命令窗口(Command Window ),是MATLAB 的主窗口,默认位于MATLAB界面的右侧,用于输入命令、运行命令并显示运行结果。

(完整版)MATLAB)课后实验[1]

(完整版)MATLAB)课后实验[1]

实验一 MATLAB 运算基础1. 先求以下表达式的值,尔后显示 MATLAB 工作空间的使用情况并保存全部变量。

(1)2sin 85 z1 21 e(2) 12z ln( x 1 x ) ,其中22 x2 1 2i5(3)ae e az sin( a 0.3) ln , a 3.0, 2.9, L , 2.9, 32 22t 0 t 1(4) 2z t 1 1 t 242t 2t 1 2 t 3,其中解:M 文件:z1=2*sin(85*pi/180)/(1+exp(2))x=[2 1+2*i;-.45 5];z2=1/2*log(x+sqrt(1+x^2))a=-3.0:0.1:3.0;z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)t=0:0.5:2.5;z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1)4. 完成以下操作:(1) 求[100,999] 之间能被 21 整除的数的个数。

(2) 建立一个字符串向量,删除其中的大写字母。

解:(1) 结果:m=100:999;n=find(mod(m,21)==0);length(n)ans =43(2). 建立一个字符串向量比方:ch='ABC123d4e56Fg9'; 那么要求结果是:ch='ABC123d4e56Fg9';k=find(ch>='A'&ch<='Z');ch(k)=[]ch =123d4e56g9实验二 MATLAB矩阵解析与办理1. 设有分块矩阵 A E R3 3 3 2O S2 3 2 2,其中 E、R、O、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵,试经过数值计算考据 2A E R RS2O S。

MATLAB数值实验一(数据的插值运算及其应用完整版)

MATLAB数值实验一(数据的插值运算及其应用完整版)

佛山科学技术学院实 验 报 告课程名称 数值分析 实验项目 插值法与数据拟合 专业班级 机械工程 姓 名 余红杰 学 号 10 指导教师 陈剑 成 绩 日 期 月 日一、实验目的1、学会Lagrange 插值、牛顿插值和三次样条插值等基本插值方法;2、讨论插值的Runge 现象3、学会Matlab 提供的插值函数的使用方法,会用这些函数解决实际问题。

二、实验原理1、拉格朗日插值多项式2、牛顿插值多项式3、三次样条插值 三、实验步骤1、用MATLAB 编写独立的拉格朗日插值多项式函数2、用MATLAB 编写独立的牛顿插值多项式函数3、用MATLAB 编写独立的三次样条函数(边界条件为第一、二种情形)4、已知函数在下列各点的值为:根据步骤1,2,3编好的程序,试分别用4次拉格朗日多项式4()L x 、牛顿插值多项式4()P x 以及三次样条函数()S x (自然边界条件)对数据进行插值,并用图给出 {(,),0.20.08,0,1,2,,10i i i x y x i i =+=},4()L x 、4()P x 和()S x 。

5、在区间[-1,1]上分别取10,20n =用两组等距节点对龙格函数21(),(11)125f x x x=-≤≤+作多项式插值,对不同n 值,分别画出插值函数及()f x 的图形。

6、下列数据点的插值可以得到平方根函数的近似,在区间[0,64]上作图。

(1)用这9个点作8次多项式插值8()L x 。

(2)用三次样条(第一边界条件)程序求()S x 。

7、对于给函数21()125f x x =+在区间[-1,1]上取10.2(0,1,,10)i x i i =-+=,试求3次曲线拟合,试画出拟合曲线并打印出方程,与第5题的结果比较。

四、实验过程与结果:1、Lagrange 插值多项式源代码:function ya=lag(x,y,xa) %x 所有已知插值点 %y 插值点对应函数值 %xa 所求点,自变量 %ya 所求点插值估计量 ya=0; mu=1; %初始化%循环方式求L 系数,并求和: for i = 1:length(y) for j = 1:length(x) if i ~= jmu = mu * (xa - x(j) ) / ( x(i) - x(j) ); else continue end endya = ya + y(i) * mu ; mu = 1; end2、Newton 源代码:function ya = newton(x,y,xa) %x 所有已知插值点 %y 插值点对应函数值 %xa 所求点,自变量 %ya 所求点插值估计量 %建立系数零矩阵D 及初始化:D = zeros(length(x)-1);ya = y(1);xi = 1;%求出矩阵D,该矩阵第一行为牛顿插值多项式系数:for i=1:(length(x)-1)D(i,1) = (y(i+1) -y(i))/(x(i+1) -x(i));endfor j=2:(length(x)-1)for i=1:(length(x)-j)D(i,j) = (D(i+1,j-1) - D(i,j-1)) / (x(i+j) - x(i)); endend%xi为单个多项式(x-x(1))(x-x(2))...的值for i=1:(length(x)-1)for j=1:ixi = xi*(xa - x(j));endya = ya + D(1,i)*xi;xi = 1;end3、三次样条插值多项式(1)(第一边界条件)源代码:function y=yt1(x0,y0,f_0,f_n,x) _____________(1)%第一类边界条件下三次样条插值;%xi 所求点;%yi 所求点函数值;%x 已知插值点;%y 已知插值点函数值;%f_0左端点一次导数值;%f_n右端点一次导数值;n = length(x0);z = length(y0);h = zeros(n-1,1);k=zeros(n-2,1);l=zeros(n-2,1);S=2*eye(n);for i=1:n-1h(i)= x0(i+1)-x0(i);endfor i=1:n-2k(i)= h(i+1)/(h(i+1)+h(i));l(i)= 1-k(i);end%对于第一种边界条件:k = [1;k]; _______________________(2)l = [l;1]; _______________________(3)%构建系数矩阵S:for i = 1:n-1S(i,i+1) = k(i);S(i+1,i) = l(i);end%建立均差表:F=zeros(n-1,2);for i = 1:n-1F(i,1) = (y0(i+1)-y0(i))/(x0(i+1)-x0(i));endD = zeros(n-2,1);for i = 1:n-2F(i,2) = (F(i+1,1)-F(i,1))/(x0(i+2)-x0(i));D(i,1) = 6 * F(i,2);end%构建函数D:d0 = 6*(F(1,2)-f_0)/h(1); ___________(4)dn = 6*(f_n-F(n-1,2))/h(n-1); ___________(5)D = [d0;D;dn]; ______________(6)m= S\D;%寻找x所在位置,并求出对应插值:for i = 1:length(x)for j = 1:n-1if (x(i)<=x0(j+1))&(x(i)>=x0(j))y(i) =( m(j)*(x0(j+1)-x(i))^3)/(6*h(j))+...(m(j+1)*(x(i)-x0(j))^3)/(6*h(j))+...(y0(j)-(m(j)*h(j)^2)/6)*(x0(j+1)-x(i))/h(j)+... (y0(j+1)-(m(j+1)*h(j)^2)/6)*(x(i)-x0(j))/h(j) ; break;else continue;endendend(2)(自然边界条件)源代码:仅仅需要对上面部分标注的位置做如下修改:__(1):function y=yt2(x0,y0,x)__(2):k=[0;k]__(3):l=[l;0]__(4)+(5):删除—(6):D=[0:D:0]4、——————————————PS:另建了一个f方程文件,后面有一题也有用到。

实验一 matlab的基本操作

实验一 matlab的基本操作

实验一:MATLAB 的基本操作 实验名称: MA TLAB 操作实验日期: 2020 年 4 月 18 日姓名: 班级学号:成绩:一、实验目的1、熟悉MATLAB 的界面和基本操作;2、掌握MATLAB 的基本运算方法;3、掌握MATLAB 中帮助命令的使用方法。

二、实验内容及步骤1、进入Matlab 工作环境,熟悉各窗口的功能。

(1)双击桌面图标,或从“开始”菜单打开MATLAB.exe ,启动MATLAB 。

(2)查看MA TLAB 界面各窗口的布局、了解其功能,并完成各窗口之间的切换。

(3)设置当前工作目录。

在D 盘创建mymatlab 目录,并将其作为当前工作目录。

今后的实验过程中以此目录作为当前工作目录。

2、计算225.389.1753cos 54.5e -÷⎪⎭⎫ ⎝⎛+π的值。

(1)在命令窗口(Command Window )中输入程序:5.54^2+cos(3/5*pi)*sqrt(17.89)/3.5-exp(2)(2)按回车键运行,如果出现“Error ”(出错信息),则应找出原因并改正,再运行。

(3)运行结果: 22.9291 3、输入矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=913652824A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333222111B ,在命令窗口中执行下列表达式,掌握其含义:A(2,3) A(:,2) A(3,:) A(:,1:2:3) A(:,3).*B(:,2) A(:,3)*B(2,:) A*B A.*B A^2 A.^2 B/A B./A(1)在命令窗口键入 A=[4,2,8;2,5,6;3,1,9] ,生成矩阵A ;键入B=[1,1,1;2,2,2;3,3,3] ,生成矩阵B ;(2)记录执行以下表达式的结果:A(2, 3)= 6 A(6)= 1 A(:,2)= [2;5; 1] A(3,:) = [3 1 9] A(:,1:2:3)= [4 8;2 6;3 9] A(:,3).*B(:,2)= [8;12;27]A(:,3)*B(2,:)= [16 16 16;12 12 12;18 18 18]A*B= [32 32 32;30 30 30;32 32 32] A.*B= [4 2 8;4 10 12;9 3 27] A^2= [44 26 116;36 35 100;41 20 111]A.^2= [16 4 64;4 25 36;9 1 81]B/A= [0.5000 0.0769 -0.3846;1.0000 0.1538 -0.7692;1.5000 0.2308 -1.1538]B./A= [0.2500 0.5000 0.1250;1.000 0.4000 0.3333;1.0000 3.0000 0.3333]B.\A= [4.0000 2.0000 8.0000;1.0000 2.5000 3.0000;1.0000 0.3333 3.0000]4、产生一个5阶魔方矩阵,将矩阵的第3行4列元素赋值给变量a ;将由矩阵第2,3,4行第2,5列构成的子矩阵赋值给变量b 。

实验一 MATLAB运算基础

实验一 MATLAB运算基础

实验一MATLAB运算基础一、实验目的1.熟悉启动和退出MATLAB的方法。

2.熟悉MATLAB命令窗口的组成。

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

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

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

(1)z1=2sin85°1+e2(2)z2=12ln(x+1+x2),其中x=21+2i−0.455(3)z3=e0.3a−e−0.3a2sin(a+0.3)+ln0.3+a2,a=−3.0,−2.9,…,,2.9,3.0(4)z4=t2,t2-1,t2-2t+1,0≤t<11≤t<22≤t<3,其中t=0:0.5:2.5解:(1)z1=2*sin(85*pi/180)/(1+exp(1)*exp(1));(2)x=[2,1+2i;-0.45,5];z2=0.5*log(x+sqrt(1+x*x));(3)a=-3.0:0.1:3.0;z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2); (4)t=0:0.5:2.5;z4=t.^2-(1-2.*t).*(t<3&t>=2)-(t<2&t>=1);运行结果:z1 =0.2375z2 =0.7114 - 0.0253i 0.8968 + 0.3658i0.2139 + 0.9343i 1.1541 - 0.0044iz3=Columns 1 through 270.7388 + 3.1416i 0.7696 + 3.1416i 0.7871 + 3.1416i0.7913 + 3.1416i 0.7822 + 3.1416i 0.7602 + 3.1416i0.7254 + 3.1416i 0.6784 + 3.1416i 0.6196 + 3.1416i0.5496 + 3.1416i 0.4688 + 3.1416i 0.3780 + 3.1416i0.2775 + 3.1416i 0.1680 + 3.1416i 0.0497 + 3.1416i-0.0771 + 3.1416i-0.2124 + 3.1416i -0.3566 + 3.1416i-0.5104 + 3.1416i -0.6752 + 3.1416i -0.8536 + 3.1416i-1.0497 + 3.1416i -1.2701 + 3.1416i -1.5271 + 3.1416i-1.8436 + 3.1416i -2.2727 + 3.1416i -2.9837 + 3.1416iColumns 28 through 61-37.0245 -3.0017 -2.3085 -1.8971-1.5978 -1.3575 -1.1531 -0.9723 -0.8083 -0.6567 -0.5151 -0.3819-0.2561 -0.1374 -0.02550.07920.1766 0.2663 0.3478 0.42060.4841 0.5379 0.5815 0.61450.6366 0.6474 0.6470 0.63510.6119 0.5777 0.5327 0.47740.4126 0.3388z4 =Columns 1 through 60 0.2500 0 1.2500 7.0000 10.2500 使用情况:a 1x61 488 doublet 1x6 48 doublex 2x2 64 double complexz1 1x1 8 doublez2 2x2 64 double complexz3 1x61 976 double complexz4 1x6 48 double2.已知:A=1234−4347873657,B=13−12033−27求下列表达式的值:(1)A+6*B和A−B+I(其中I为单位矩阵)(2)A*B和A.*B(3)A^3和A.^3(4)A/B及A\B(5)[A,B]和[A([1,3],:);B^2]解:A=[12,34,-4;34,7,87;3,65,7];B=[1,3,-1;2,0,3;3,-2,7];I=eye(3);(1)A+6*B;A-B+I;(2)A*B;A.*B(3)A^3;A.^3;(4)A/B;A\B(5)[A,B];[A([1,3],:);B^2];运行结果:(1)ans =18 52 -1046 7 10521 53 49ans =12 31 -332 8 840 67 1(2)ans =68 44 62309 -72 596154 -5 241ans =12 102 468 0 2619 -130 49(3)ans =37226 233824 48604247370 149188 60076678688 454142 118820ans =1728 39304 -6439304 343 65850327 274625 343(4)ans =16.4000 -13.6000 7.600035.8000 -76.2000 50.200067.0000 -134.0000 68.0000ans =-0.0313 0.3029 -0.33240.0442 -0.0323 0.10630.0317 -0.1158 0.1558(5)ans =12 34 -4 1 3 -134 7 87 2 0 33 65 7 3 -2 7ans =12 34 -43 65 74 5 111 0 1920 -5 403.设有矩阵A和BA=12367811121349145101516171819202122232425, B=301617−699423713−411(1)求它们的乘积C。

matlab 实验教程 实验一

matlab 实验教程 实验一

MATLAB语言及其应用实验教程实验一实验目的:1.熟悉MATLAB的界面,菜单,会使用DEMO,学会使用帮助(help)。

2.学习MATLAB的基本语法实验内容:1.打开MATLAB,点击各个菜单以了解各个子菜单项。

2.更改当前路径,在指定路径下保存所作的实验。

实验具体步骤:如图所示,在D盘根目录下建立“MATLAB实验”文件夹,并在此文件夹下以各自的学号建立子文件夹。

如图所示,点击工具栏上“Current Directory”栏的浏览按钮。

选择刚刚建立的以各自学号命名的文件夹,“Current Directory”栏显示如下路径格式。

注意!以后所有的实验均保存在此目录中。

3.学习打开和关闭命令窗口(command window),工作空间(workspace)和命令历史窗口(command history)。

实验具体步骤:单击命令窗口右上角的“X”标志,关闭命令窗口;在Desktop菜单下选择“Command window”,打开命令窗口;单击工作空间窗口右上角的“X”标志,关闭工作空间;在Desktop菜单下选择“Workspace”,打开工作空间;单击命令历史窗口右上角的“X”标志,关闭命令历史窗口;在Desktop菜单下选择“Command history”,打开命令历史窗口。

4.练习变量的赋值,包括向量赋值,矩阵赋值以及复数的赋值。

实验具体步骤:变量赋值>> a=100a =100>> b=0.1b =0.1000向量赋值>> a=1:1:10a =1 2 3 4 5 6 7 8 9 10>> b=1:2:10b =1 3 5 7 9矩阵赋值>> a=[1 2 3;4 5 6;7 8 9]a =1 2 34 5 67 8 9复数赋值>> c=3+5.2ic =3.0000 + 5.2000i>> z=[1+2i,3+4i;5+6i,7+8i]z =1.0000 +2.0000i3.0000 +4.0000i5.0000 +6.0000i7.0000 +8.0000i5.用变量检查命令who和whos检查工作空间中的变量。

程序设计实验报告(matlab)

程序设计实验报告(matlab)

程序设计实验报告(matlab)实验一: 程序设计基础实验目的:初步掌握机器人编程语言Matlab。

实验内容:运用Matlab进行简单的程序设计。

实验方法:基于Matlab环境下的简单程序设计。

实验结果:成功掌握简单的程序设计和Matlab基本编程语法。

实验二:多项式拟合与插值实验目的:学习多项式拟合和插值的方法,并能进行相关计算。

实验内容:在Matlab环境下进行多项式拟合和插值的计算。

实验方法:结合Matlab的插值工具箱,进行相关的计算。

实验结果:深入理解多项式拟合和插值的实现原理,成功掌握Matlab的插值工具箱。

实验三:最小二乘法实验目的:了解最小二乘法的基本原理和算法,并能够通过Matlab进行计算。

实验内容:利用Matlab进行最小二乘法计算。

实验方法:基于Matlab的线性代数计算库,进行最小二乘法的计算。

实验结果:成功掌握最小二乘法的计算方法,并了解其在实际应用中的作用。

实验六:常微分方程实验目的:了解ODE的基本概念和解法,并通过Matlab进行计算。

实验内容:利用Matlab求解ODE的一阶微分方程组、变系数ODE、高阶ODE等问题。

实验方法:基于Matlab的ODE工具箱,进行ODE求解。

实验结果:深入理解ODE的基本概念和解法,掌握多种ODE求解方法,熟练掌握Matlab的ODE求解工具箱的使用方法。

总结在Matlab环境下进行程序设计实验,使我对Matlab有了更深刻的认识和了解,也使我对计算机科学在实践中的应用有了更加深入的了解。

通过这些实验的学习,我能够灵活应用Matlab进行各种计算和数值分析,同时也能够深入理解相关的数学原理和算法。

这些知识和技能对我未来的学习和工作都将有着重要的帮助。

MATLAB)课后实验答案

MATLAB)课后实验答案

实验一 MATLAB 运算基础1、 先求下列表达式得值,然后显示MATLAB 工作空间得使用情况并保存全部变量。

(1) 0122sin 851z e =+(2) 21ln(2z x =+,其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--L (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩,其中t =0:0、5:2、5 解: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 +⎡⎤=⎢⎥⎣⎦。

解: M 文件如下;5、 下面就是一个线性方程组:1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程得解。

(2) 将方程右边向量元素b 3改为0、53再求解,并比较b 3得变化与解得相对变化。

(3) 计算系数矩阵A 得条件数并分析结论。

解: M 文件如下:实验三 选择结构程序设计1、 求分段函数得值。

2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现,分别输出x=-5、0,-3、0,1、0,2、0,2、5,3、0,5、0时得y 值。

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

MATLAB实验指导书皖西学院信息工程学院实验一MATLAB编程环境及简单命令的执行一、实验目的1.熟悉MATLAB编程环境二、实验环境1.计算机2.MATLAB7.0集成环境三、实验说明1.首先应熟悉MATLAB7.0运行环境,正确操作2.实验学时:2学时四、实验内容和步骤1.实验内容(1)命令窗口的使用。

(2)工作空间窗口的使用。

(3)工作目录、搜索路径的设置。

(4)命令历史记录窗口的使用。

(5)帮助系统的使用。

(6)了解各菜单的功能。

2.实验步骤(1)启动MATLAB,熟悉MATLAB的桌面。

(2)进入MATLAB7.0集成环境。

(3)在命令窗口执行命令完成以下运算,观察workspace的变化,记录运算结果。

1) (365-52⨯2-70)÷3 =63.66672) >>area=pi*2.5^2 19.63503) 已知x=3,y=4,在MATLAB 中求z :()232y x y x z -= = 576 4) 将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。

m1=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡11514412679810115133216 执行以下命令 >>m1( 2 , 3 ) 10 >>m1( 11 ) 6 >>m1( : , 3 ) 3 10 6 15>>m1( 2 : 3 , 1 : 3 ) 5 11 10 9 7 6>>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1) 34 5) 执行命令>>help abs查看函数abs 的用法及用途,计算abs( 3 + 4i ) 5 6) 执行命令>>x=0:0.1:6*pi; >>y=5*sin(x);>>plot(x,y)7) 运行MATLAB 的演示程序,>>demo ,以便对MATLAB 有一个总体了解。

五、思考题1、以下变量名是否合法?为什么?(1)x2 合法(2)3col 不合法变量名只能字母开头(3)_row 不合法同上(4)for 不合法 for为系统预定义关键字2、求以下变量的值,并在MATLAB中验证。

(1)a = 1 : 2 : 5 ; a=1 3 5(2)b = [ a' a' a' ] ;1 1 13 3 35 5 5(3)c = a + b ( 2 , : )。

4 6 83、已知若需要计算t∈[-1,1],取间隔为0.01,试计算出相对应的y值。

六、实验报告要求1.整理实验结果,按顺序填写实验结果2.小结实验心得体会实验二矩阵基本运算(一)一、实验目的1.通过实验,进一步熟悉MATLAB编程环境2.通过实验掌握建立矩阵的几种方法3.通过实验理解常用的矩阵运算二、实验环境1.计算机2.MATLAB7.0集成环境三、实验说明1.首先应熟悉MATLAB7.0运行环境,正确操作2.自主编写程序,必要时参考相关资料3.实验学时:2学时四、实验内容和步骤1.实验内容(1)设A和B是两个同维同大小的矩阵,问:若A=B= 1 2 34 5 67 8 91)A*B和A.*B的值是否相等?不相等A*B= 30 36 4266 81 96102 126 150A.*B= 1 4 916 25 3649 64 812)A./B和B.\A的值是否相等?相等A./B=B.\A =1 1 11 1 11 1 13)A/B和B\A的值是否相等?不相等4)A/B和B\A所代表的数学含义是什么?数学含义是A除以B(2)写出完成下列操作的命令。

1)将矩阵A第2—5行中第1,3,5列元素赋给矩阵B。

B=A([2:5],[1,3,5]) 2)删除矩阵A的第7号元素。

A(7)=[]3)将矩阵A的每个元素值加30。

A+304)求矩阵A的大小和维数。

size(A) length(A)5)将含有12个元素的向量x转换成3*4矩阵。

reshape(A,3,4)6)求一个字符串的ASCII码。

a=’abc’ int8(a)7)求一个ASCII码所对应的字符。

char(a) a是ASCII码(3)下列命令执行后,L1、L2、L3、L4的值分别是多少?A=1:9;B=10-A;L1=A==B;L1 =0 0 0 0 1 0 0 0 0L2=A<=5;L2 =1 1 1 1 1 0 0 0 0L3=A>3&A<7;L3 =0 0 0 1 1 1 0 0 0L4=find(A>3&A<7);L4 =4 5 6(4)完成下列操作1)求[100,999]之间能被61整除的数及其个数(提示:先利用冒号表达式,再利用find和length函数。

)C=100:999;L=find(rem(C,61)==0) C(L)2)建立一个字符串向量,删除其中的大写字母(提示:利用find 函 数和空矩阵。

) b=find(a>='A'&a<='Z') a(b)=[] 2.实验步骤(1)分析实验内容,写出程序大致框架或完整的程序代码。

(2)进入MATLAB7.0集成环境。

(3)编辑程序并进行保存。

(4)运行程序,若有错误,修改错误后再次运行,如此反复进行到不显示出错为止。

(5)检查程序输出结果。

五、思考题用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] inv(a)*b2) ⎪⎪⎩⎪⎪⎨⎧=-++=--=-++=++56533332821w z y x w y x w z y x z y xa=[1,1,1,0;1,2,1,-1;2,-1,0,-3;3,3,5,-6] b=[1;8;3;5] inv(a)*b六、实验报告要求1.整理实验结果,填入相应表格中2.小结实验心得体会实验三矩阵基本运算(二)一、实验目的1.熟悉矩阵和向量的建立方式2.理解矩阵拆分的方法3.通过实验进一步掌握矩阵的基本运算二、实验环境1.计算机2.MATLAB7.0集成环境三、实验说明1.熟练操作MATLAB7.0运行环境 2.自主编写程序,必要时参考相关资料 3.实验前应写出程序大致框架或完整的程序代码 5.实验学时:2学时四、实验内容和步骤1.实验内容(1)已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=7613870451A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=023352138B 求下列表达式的值:1) A+6B 和A 2-B+I (I 为单位矩阵) A+6B=47 23 -10 12 37 26 -15 73 7 A 2-B+I=-18 -217 17 22 533 109 21 867 526 A*B=14 14 16-10 51 21125 328 180A.*B=-8 15 40 35 24-9 122 0B*A=-11 0 -157 228 533 -1 283)A/B和B\AA/B=1.2234 -0.92552.9787-0.9468 2.3511 -0.95744.6170 3.8723 13.8936B\A-0.5106 -8.6170 -1.12770.7340 17.5745 1.8085-0.8830 -21.2128 0.40434)[A,B]和[A([1,3],:);B^2][A,B]-1 5 -4 8 3 -10 7 8 2 5 33 61 7 -3 2 0 [A([1,3],:);B^2]-1 5 -43 61 773 37 117 37 13-20 1 9(2)已知⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---=14.35454.9632053256545410778.01023A ,取出其前三行构成矩阵B ,其前两列构成矩阵C ,其右下角3×2子矩阵构成矩阵D ,B 与C 的乘积构成矩阵E ,分别求E<D,E&D,E|D 、~E|~D 和find(A>=10&A<25) B=A([1:3],[1:4]) C=A([1:4],[1,2]) D=A([2:4],[3,4]) E= 1.0e+003 *0.9141 -0.2239 1.2080 2.7123 1.1330 -0.2103 E<D0 1 0 0 0 1 E&D1 1 0 1 1 1 E|D1 1 1 1 1 1 ~E|~D 0 0 1 0 0 0 find(A>=10&A<25)1 5(1)分析实验内容,写出程序大致框架或完整的程序代码。

(2)进入MATLAB7.0集成环境。

(3)编辑程序并进行保存。

(4)运行程序,若有错误,修改错误后再次运行,如此反复进行到不显示出错为止。

(5)检查程序输出结果。

五、思考题109910101022222++++==---=∑ n n y ,求y=?(循环语句做) x=0;for n=-10:1:10 y = 2^n; x = x+y; endx= 2.0480e+003六、实验报告要求1.整理实验结果,填入相应表格中 2.小结实验心得体会实验四分支结构程序设计一、实验目的1.理解MATLAB的编程方法2.掌握if语句和switch语句的使用3. 掌握分支结构程序设计方法4. 通过练习理解MATLAB编程方法和分支结构程序设计方法。

二、实验环境1.计算机2.MATLAB7.0集成环境三、实验说明1.熟练操作MATLAB7.0运行环境2.自主编写程序,必要时参考相关资料3.实验前应写出程序大致框架或完整的程序代码5.实验学时:3学时四、实验内容和步骤1.实验内容(1)从键盘输入1到5000之内的一个数,例如输入693,输出为396,如果数据超出范围,请输出提示“输入数据有误,请输入1~5000之内的数”将它反向输出,如果数据是个位数,则直接输出。

while 1x=input('请输入1~5000之内的一个数 x:')if x>=0&x<10disp(x);elsey = 0;while x~=0x=fix(x);y = y*10 + rem(x,10);x =fix(x/10);endendif y>5000disp('超出范围请重新输入:');elsebreak;endyend(2)输入一个百分制成绩,要求输出成绩等级A,B,C,D,E其中90-100位A,80-89为B,70-79为C,60-69为D,60以下为E1)分别用if语句和switch语句实现2)输入百分制成绩后要判断成绩的合理性,对不合理的成绩应输出出错信息(若成绩出现小数,则只能是“.5”)while 1x=input('请输入成绩:');x=x/10;x=fix(x);switch xcase 10disp('A')case 9disp('A')case 8disp(‘B’)case 7disp('C')case 6disp('D')otherwiseif x>0&x<=10disp('E')elsedisp('Error!')endendend(3)某商场对顾客所购买的商品实行打折销售,标准如下(商品价格用price 来表示):price<200 没有折扣;200≤price<500 3%折扣;500≤price<1000 5%折扣;1000≤price<2500 8%折扣;2500≤price<5000 10%折扣;5000≤price 14%折扣,输入所售商品的价格,求其实际销售价x = input('Please input the price:');if x<200disp(x);else if x<500x = x*0.97;disp(x);else if x<1000x = x*0.95;disp(x);else if x<2500x = x*0.92;disp(x);else if x<5000x = x*0.9;disp(x);elsex = x*0.86;disp(x);endendendendend2.实验步骤(1)分析实验内容,写出程序大致框架或完整的程序代码。

相关文档
最新文档