实验1 Matlab程序设计基础
- 格式:doc
- 大小:146.00 KB
- 文档页数:4
实验项目1 MATLAB熟悉使用及编程基础MATLAB是美国Mathworks公司推出的一套高件能的数值分析和计算软件,它将矩阵运算、数值分析、图形处理、编程处术结合在一起,为用户提供了一个强有力的科学及工程问题分析计算和程序设计的工具。
MATLAB本身也在不断改进和创新,特别是2000年以出的版本6,无论在界面设计、计算方法、编程阶段和工具等方面都有了巨大的突破,全面引入了面向对象编程的概念和方法,使MATLAB真正成为了具有全部高级语言功能和特征的新一代软件开发平台。
MATLAB开发环境是一组工具和组件的集成,这些工具是图形化的用户接口,它们包括 MATLAB桌面、命令窗口、命令历史窗口、编辑调试窗口以及帮助信信息、工作空间、文件和搜索路径等浏览器。
MATLAB集成了丰富的数学函数库,其强大的计算能力覆盖了从基本函数(如求和、正弦、余弦和复数运算等)到特殊函数(如矩阵求逆、矩阵特征值、贝塞尔函数和快速傅里叶变换等)的范围。
MATLAB语言是一种高级编程语言,包括控制流的描述、函数、数据结构、输入输出及面对对象编程,既可以编制适用于快速使用的小程序,也可以编制大型复杂的应用程序。
MATLAB提供了功能强大的图形系统,既可以完成二维和三维数据的可视化、图像处理、动画和图形表达等功能,也可以定制图形的外观,如建立一个完整的图形用户界由的应用程序。
1.实验目的:掌握MA TLAB编程语言和偏微分方程PDE工具箱对电磁场的基本问题进行仿真;2 实验内容:1、MATLAB启动MA TLAB安装到硬盘后,启动方法有:(1)点击Windows桌面上自动生成的快捷方式图标;(2)点击matlab 6\文件夹下快捷方式图标令MATLAB(3)点击matlab/bin/win32文件夹中的Matlab.exe2、m文件编写熟悉基本指令、基本数学运算,借助help、lookfor等命令实现在线帮助。
建立M文件的一般步骤如下:(1)打开文件编辑器:指的是MA TLAB内部编辑/调试器,可以有几种不同的方法打开文件编辑器,最简单的方法是在操作桌面助工具栏上选择蹦(建立新文件)或选择瞪(打开已有的文件),也可以在命令窗口输入命令edit建立新文件或输入命令edit filename,打开名为fi1ename的M文件。
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之和。
MATLAB基础实验指导书漳州师范学院物电系2010年10月目录实验一MATLAB环境的熟悉与基本运算 (2)实验二MATLAB数值运算 (8)实验三MATLAB语言的程序设计 (12)实验四MATLAB的图形绘制 (12)实验五MATLAB应用实例实验六采用SIMULINK的系统仿真 (17)实验七控制系统的频域与时域分析 (25)实验一 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 数组操作函数三、实验内容1、新建一个文件夹(自己的名字命名)2、启动MATLAB6.5,将该文件夹添加到MATLAB路径管理器中。
方法如下:3、保存,关闭对话框4、学习使用help命令,例如在命令窗口输入help eye,然后根据帮助说明,学习使用指令eye(其它不会用的指令,依照此方法类推)5、学习使用clc、clear,观察command window、command history和workspace等窗口的变化结果。
数学与软件科学学院实验报告学期: 2015 至2016 第学期年月日课程名称:MATLAB程序设计及应用专业: 信息与计算科学级班实验编号:1实验项目MATLAB的基本功能熟悉与应用指导教师__姓名学号:实验成绩:一、实验目的及要求(1)熟悉所用计算机系统的基本操作方法;(2) 学习基本MATLAB2014的界面以及基础菜单的应用;(3) 通过执行简单的MATLAB命令,了解该语言方便快捷的特点;(4)通过简单的M文件,了解M文件的编辑、保存和运行过程.二、实验内容(1) 熟悉MATLAB基础菜单的应用:File 菜单Eidit菜单Debug菜单Desktop菜单Help菜单(2) 熟悉几个基本窗口的特征和功能:Command Window、Command History、Workspace、Current Directory等。
(3) Command Window的功能练习:1>Command Window直接计算数学表达式功能,采用课本第二章课后习题(1-4)2>使用help 命令查找函数sin,cos,sqrt,roots,polyval,polyfit,interp1文件,翻译每个函数的功能,用法,并用一个例子给出验证(4) Workspace 窗口的应用:在此窗口下查看变量的内容,修改编辑变量的相关信息。
(5)编写程序完成下列题目:1> 求和S=1+2+3+ (100)2 >求和e 1111!2!10!1...=++++3>求和S 1112310!1...=++++4>设A 234576138⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦, 求A 的逆、特征值和特征向量;验证Ax=λx 三、实验结果(该部分不够填写.请填写附页)结果一:Command Window 的功能练习*课本第二章课后习题(1-4)*1>> a=87(1) a =87>> b=190b =190>> a+bans =277>> a-bans =-103>> a*bans =16530(2) >> A=uint8(87)A =87>> B=uint8(190) B =190>> A+bans =255>> A+Bans =255>> A-Bans =>> A*Bans =2552.>> sind(60) ans =0.8660 >> exp(3) ans =20.0855 >> cos(3*pi/4) ans =-0.70713.>> u=2u =2>> v=3v =3>> 4*u*v/log(v)ans =21.8457>> (exp(u)+v)^2/(v^2-u) ans =15.4189>> sqrt(u-3*v)/(u*v)ans =0.0000 + 0.4410i4.(1)>> (3-5*i)*(4+2*i)ans =22.0000 -14.0000i (2)>> sin(2-8*i)ans =1.3553e+03 + 6.2026e+02i结果二:编写程序完成下列题目:1> 求和S=1+2+3+ (100)2 >求和e 1111!2!10!1...=++++3>求和S 1112310!1...=++++4>设A 234576138⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦, 求A 的逆、特征值和特征向量;验证Ax=λx (1)(2)(3)(4)注:实验成绩等级分为(90-100分)优,(80-89分)良,(70-79分)中,(60-69分)及格,(59分)不及格。
实验一、MATLAB基本操作一、基本操作1、命令窗口的简单使用(1)简单矩阵的输入(2)求[12+2×(7-4)]÷32的算术运算结果2、有关向量、矩阵或数组的一些运算(1)设A=15;B=20;求C=A+B与c=a+b?(2)设A=[1 2 3;4 5 6;7 8 9],B=[9 8 7;6 5 4;3 2 1];求A*B与A.*B?说明*与.*的运算特点A=[1 2 3;4 5 6;7 8 9],B=[9 8 7;6 5 4;3 2 1](3)设a=10,b=20;求i=a/b与j=a\b?a=10,b=20(4)设a=[1 -2 3;4 5 -4;5 -6 7];请设计出程序,分别找出小于0的矩阵元素及其位置。
(5)在MATLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8];看结果如何?如果改成运行A=[1,2;3,4]+i[5,6;7,8],结果又如何?(6)请写出完成下列计算的指令:a=[1 2 3;3 4 2;5 2 3],求a^2=?,a.^2=?(7)有一段指令如下,请思考并说明运行结果及其原因clearX=[1 2;8 9;3 6];X(:)矩阵变为一维矩阵使用三元组方法,创建下列稀疏矩阵2 0 8 00 0 0 10 4 0 06 0 0 0(8)写出下列指令的运行结果>> A = [ 1 2 3 ]; B = [ 4 5 6 ];>> C = 3.^A >> D = A.^B C =3 9 27 D =1 32 7293、 已知⎪⎭⎫ ⎝⎛+⋅=-334sin 234πt e y t 若需要计算t ∈[-1,1],取间隔为0.01,试计算出相对应的y 值。
二、运算基础1、 设有矩阵A 和B ,A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;1617 18 19 20;21 22 23 24 25],B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]; 1) 求它们的乘积C2)将矩阵C的右下角3x2子矩阵赋给D2、完成下列操作1)求[100,999]之间能被61整除的数及其个数(提示:先利用冒号表达式,再利用find和length函数。
实验1 MA TLAB环境及命令窗口的使用一、实验目的1.熟悉MATLAB的桌面环境;2.熟悉MATLAB的各类窗口;3.熟悉常用选项和工具;二、实验步骤1.启动MATLAB2.设置当前路径为d:\dzxx\mat3.在命令窗口输入a=ones(3);b=[1 2 3;4 -5 6;pi 8 9];c=’计算’;d=a+b*j;4.设置数据显示格式,分别在format +,short e,short g,long情况下查看数值形式。
5.将命令窗口的字体设置为宋体,字体大小10号。
6.通过命令窗口左下侧的“start”按钮—“Desktop Tools”—“workspace”,打开工作窗口(1)使用who和whos命令查看变量内容,并记录;(2)使用“clear”命令删除变量b(3)使用“exist”命令查看变量a,b,c,d,记录结果(4)在工作空间选取变量a,双击编辑数组,改为[1 2 3;1 1 1;1 1 1] 7.打开历史窗口,察看输入过的命令,选中b=[1 2 3;4 -5 6;pi 8 9],单击右键,通过“evaluate selection”将其导入命令窗口运行,查看相应结果。
8.选择a=ones(3);b=[1 2 3;4 -5 6;pi 8 9];c='计算'; d=a+b*j四行命令,单击鼠标右键,选择“Creat M-file”命令生成M文件。
9.将M文件保存在“d:\dzxx\mat\test1.m,可以在“当前目录浏览器”窗口看到test1.m,在命令窗口输入“test1”,运行该文件。
10. 在命令窗口输入“save test1”命令,在“当前目录浏览器”窗口看到生成”test1.mat”数据文件。
11. 在命令窗口输入“exit”命令退出MA TLAB12. 重启MATLAB,设置当前路径为d:\dzxx\mat,在命令窗口输入“test1”命令可运行该文件,键入“type test1”可以查看文件内容,记录之。
MATLAB程序设计软件实验报告专业及班级____通信中兴131_______姓名____魏增_______________学号_____6102213869________日期_____2015.6.15_________南昌大学实验报告学生姓名: 魏增 学 号: 6102213869 班级: 中兴131班 实验类型:□ 验证 □ 综合 ■ 设计 □ 创新 实验日期: 实验成绩:实验一 MA TLAB 的基本使用一、 实验目的1.了解MA TALB 程序设计语言的基本特点,熟悉MA TLAB 软件的运行环境;2.掌握变量、函数等有关概念,掌握M 文件的创建、保存、打开的方法,初步具备将一般数学问题转化为对应计算机模型处理的能力;3.掌握二维图形绘制的方法,并能用这些方法实现计算结果的可视化。
二、 MATLAB 的基础知识通过本课程的学习,应基本掌握以下的基础知识: 一. MA TLAB 简介二. MA TLAB 的启动和退出 三. MA TLAB 使用界面简介 四. 帮助信息的获取五. MA TLAB 的数值计算功能六. 程序流程控制 七. M 文件八. 函数文件九. MATLAB 的可视化 三、上机练习1. 仔细预习第二部分内容,关于MA TLAB 的基础知识。
2. 熟悉MA TLAB 环境,将第二部分所有的例子在计算机上练习一遍3、已知矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=123456789,987654321B A 。
求A*B ,A .* B ,比较二者结果是否相同。
并利用MA TLAB 的内部函数求矩阵A 的大小、元素和、长度以及最大值。
解:>> A=[1 2 3;4 5 6;7 8 9]; >> B=[9 8 7;6 5 4;3 2 1];>> A*Bans =30 24 18 84 69 54 138 114 90 >> A.*B ans =9 16 2124 25 2421 16 9 两者结果不同 >> [m,n]=size(A) m =3 n =3 >> b=sum(A) b =12 15 18 >> a=length(A) a = 3 >>max(A)ans =7 8 94、Fibonacci 数组的元素满足Fibonacci 规则:),2,1(,12=+=++k a a a k k k ;且121==a a 。
实验一:MATLAB基本操作一、实验目的1、学习掌握MA TLAB语言的基本操作方法2、掌握命令窗口的使用3、熟悉MA TLAB的数据表示、基本运算和程序控制语句4、熟悉MA TLAB程序设计的基本方法二、实验内容和要求2.b.帮助命令的使用,查找SQRT函数的使用方法按Start-help-search,输入sqrtc.矩阵运算(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结果为:ans =105 115229 251(2)矩阵除法已知A=[1 2 3 ;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];求A/B,A\BA/B=1.0000 1.0000 1.00004.0000 2.5000 2.00007.0000 4.0000 3.0000A\B=1.0e+016 *0.3152 -1.2609 0.9457-0.6304 2.5218 -1.89130.3152 -1.2609 0.9457(3)矩阵的转置及共轭转置已知A=[15+i,2-i,1;6*i,4,9-i];求A.',A'A=15.0000 + 1.0000i 2.0000 - 1.0000i 1.00000 + 6.0000i 4.0000 9.0000 - 1.0000iA.'=15.0000 + 1.0000i 0 + 6.0000i2.0000 - 1.0000i 4.00001.0000 9.0000 - 1.0000iA'=15.0000 - 1.0000i 0 - 6.0000i2.0000 + 1.0000i 4.00001.0000 9.0000 + 1.0000i(4)使用冒号选出指定元素已知A=[1 2 3;4 5 6;7 8 9];求A中第三列前两个元素;A中所有第二行的元素A =1 2 34 5 67 8 9求A中第三列前两个元素:A([1,2],3)ans =36求A中所有第二行的元素:A(2,1:1:end)ans =4 5 6(5)方括号[]用magic函数生成一个4阶魔方矩阵,删除该矩阵的第四列magic(4)ans =16 2 3 135 11 10 89 7 6 124 14 15 1a=ans;a(:,[1,2,3])ans =16 2 35 11 109 7 64 14 153..多项式(1)求多项式p(X)=x^3-2x-4的根(2)已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4]求矩阵A的特征多项式;求矩阵多项式中未知数为20时的值;把矩阵A作为未知数代入到多项式中;(1)a=[1,0,-2,-4]x=roots(a)x =2.0000-1.0000 + 1.0000i-1.0000 - 1.0000i(2)求特征多项式系数:A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4]p=poly(A)p =1.0000 -6.9000 -77.2600 -86.1300 604.5500求矩阵多项式中未知数为20时的值:c=polyval(p,20)c =7.2778e+004把矩阵A作为未知数代入到多项式中:c=polyval(p,A)c =1.0e+003 *0.3801 -0.4545 -1.9951 0.4601-1.9951 0.2093 -1.9951 -2.8880-0.4545 -4.8978 0.6046 0.43530.4353 0.0840 -0.4545 -1.16174、基本程序设计(1)编写命令文件:计算1+2+3+……+n<2000;n=0;s=0;while s<2000n=n+1;s=s+n;endn=n-1n =62(2)编写函数文件:分别用for和while循环结构编写程序,求2的0到n次幂的和。
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运算基础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 +⎡⎤=⎢⎥⎣⎦。
一、实验目的1.熟悉MATLAB 的使用方法及特点;学会建立MATLAB 搜索路径;熟悉MATLAB 工作空间、MATLAB 集成环境、命令窗口;掌握MATLAB 的通用命令、管理命令和函数、管理变量和工作空间的使用方法;2.掌握MATLAB 基本操作及矩阵基础知识,包括:输入矩阵、矩阵的转置、矩阵元素求和、矩阵下际、矩阵连接、矩阵行列删除、矩阵产生和操作、逻辑和关系运算、操作符和特殊字符、基本矩阵和矩阵操作、基本矩阵和阵列;掌握特殊变量和常数;掌握基本数学函数。
二、实验内容1.利用基本矩阵产生 3x3 和15x8 的单位阵,全1 阵,全0 阵,均匀分布的随机阵([-1,1]之间),正态分布随机阵(方差4,均值1)eye(3),eye(15,8)ones(3),ones(15,8)zeros(3),zeros(15,8)1-2*rand(3),1-2*rand(15,8)2*randn(3)+1, 2*randn(15,8)+12.利用diag()函数和rot90()产生下列矩阵:0 0 8 2 0 4a= 0 -7 5 b= 0 5 02 3 0 7 0 8然后求解a 阵的逆矩阵aa 及b 阵的特征值和对应特征向量,并利用reshape 将aa 阵变换成行向量。
X=diag([2 -7 8]),Y=diag([3 5],-1),a=rot90(X+Y)M=diag([2 5 8]),N=diag([7 0 4]),b=M+rot90(N)aa=inv(a) reshape(aa,1,9) d=eig(b) [V,D]=eig(b)aa =-0.1339 0.2143 0.50000.0893 -0.1429 00.1250 0 0ans = -0.1339 0.0893 0.1250 0.2143 -0.1429 0 0.5000 0 0d =-1.082811.08285.0000V = -0.7921 -0.4030 00 0 1.00000.6104 -0.9152 0D = -1.0828 0 00 11.0828 00 0 5.00003产生一均匀分布在(-5,5)随机阵(50x2),精确到小数点后一位x=5-round(100*rand(50,2))/104.编程实现当α∈[-π,π],间隔为1o 时,求解正弦和余弦的值,并利用plot()函数绘制正弦,余弦曲线。
实验1MATLAB 使用方法和程序设计一、实验目的1、掌握MATLAB 软件使用的基本方法。
2、熟悉MATLAB 的数据表示、基本运算和程序控制语句。
3、熟悉MATLAB 绘图命令及基本绘图控制。
4、熟悉MATLAB 程序设计的基本方法。
二、实验内容1.帮助命令使用Help 命令,查找sqrt (开方)函数的使用方法。
2、矩阵运算(1)矩阵乘法已知A=[12;34];B=[55;78];求A^2*B 。
(2)矩阵除法已知A=[123;456;789];B=[100;020;003];求A\B,A/B 。
(3)矩阵的转置及共轭转置已知A=[5+i ,2-i ,1;6*i ,4,9-i];求A.’,A’。
(4)使用冒号选出指定元素已知A=[123;456;789];求A 中第3列前2个元素;A 中所有列第2,3行的元素。
(5)方括号[]用magic 函数生成一个4阶魔术矩阵,删除该矩阵的第四列。
3、多项式(1)求多项式p(x)=x 3-2x-4的根。
(2)求f (x)=(cos x)2的一次导数。
(3)求微分方程的通解,并验证。
4、基本绘图命令(1)绘制余弦曲线cost =y ,]2,0[π∈t 。
(2)在同一坐标系中绘制余弦曲线0.25)-cos(t =y 和正弦曲线)5.0sin(-=t y ,]2,0[π∈t (3)用plot3函数绘制三维螺线:22xdy xy xe dx -+= sin()cos()x t y t z t =⎧⎪=⎨⎪=⎩(0<t <20)5、基本绘图控制4,0[ 区间上的X1=10sint曲线,并要求绘制](1)线形为点划线,颜色为红色,数据点标记为加号。
(2)使用axis坐标轴控制:显示范围、刻度线、比例、网络线。
(3)使用title、text、xlabel、legend等函数进行标注控制:坐标轴名称、标题、相应文本。
6、MATLAB基本程序设计(1)编写命令文件:计算1+2+…+n<2000时的最大n值。
实验一:Matlab程序设计及绘图一、实验目的1、掌握 Matlab 软件使用的基本方法;2、熟悉 Matlab 的编程方法;3、熟悉 Matlab 绘图命令及基本绘图控制。
二、实验内容1.基本绘图命令(1)绘制余弦曲线 y=cos(t),t∈[0,2π](2)在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5), t∈[0,2π]2.基本绘图控制绘制[0,4π]区间上的 x1=10sint 曲线,并要求:(1)线形为点划线、颜色为红色、数据点标记为加号;(2)坐标轴控制:显示范围、刻度线、比例、网络线(3)标注控制:坐标轴名称、标题、相应文本;3.程序设计:找出从2开始的前100个素数,并依次排列显示。
三、实验原理MATLAB环境是一种为数值计算、数据分析和图形显示服务的交互式的环境。
MATLAB有3种窗口,即:命令窗口(The Command Window)、m-文件编辑窗口(The Edit Window)和图形窗口(The Figure Window),而Simulink另外又有Simulink模型编辑窗口。
1.命令窗口(The Command Window)当MATLAB启动后,出现的最大的窗口就是命令窗口。
用户可以在提示符“>>”后面输入交互的命令,这些命令就立即被执行。
在MATLAB中,一连串命令可以放置在一个文件中,不必把它们直接在命令窗口内输入。
在命令窗口中输入该文件名,这一连串命令就被执行了。
因为这样的文件都是以“.m”为后缀,所以称为m-文件。
2.m-文件编辑窗口(The Edit Window)我们可以用m-文件编辑窗口来产生新的m-文件,或者编辑已经存在的m-文件。
在MATLAB主界面上选择菜单“File/New/M-file”就打开了一个新的m-文件编辑窗口;选择菜单“File/Open”就可以打开一个已经存在的m-文件,并且可以在这个窗口中编辑这个m-文件。
实验1 Matlab 程序设计基础
一、实验名称:Matlab 程序设计基础.
二、实验目的:掌握Matlab 语言,会用Matlab 进行一些简单的程序设计.
三、实验要求:
1、对以下问题,编写M 文件:
(1)用起泡法对10个数由小到大排序.即将相邻两个数比较,将小的调到前头.
(2)有一个4×5矩阵,编程求其最大值及其所处的位置.
(3)编程求∑=20
1!n n .
(4)一球从100米高度自由落下,每次落地后反跳回原高度的一半,再落下.求它在第10次落地时,共经过多少米?第10次反弹有多高?
(5)有一函数f(x,y)=x^2+sin(x*y)+2*y ,写一程序,输入自变量的值,输出函数值.
2、在同一平面中的两个窗口分别画出心形线和马鞍面.
要求:
(1)在图形上加格栅、图例和标注.
(2)定制坐标.
(3)以不同角度观察马鞍面.
3、请编写绘制以下图形的Matlab 命令,并展示绘得的图形:
(1)4,12222=+=+y x y x 分别是椭圆14/2
2=+y x 的内切圆和外切圆; (2)指数函数和对数函数x y ln =的图像关于直线x y =对称.
四、实验内容:
1、1)源程序:
function p=paixu(a)
a=[1,2,5,6,8,7,9,4,3,0];
for i=9:-1:1
for j=1:i
if a(j)>a(j+1)
t=a(j);
a(j)=a(j+1);
a(j+1)=t;
end
end
end
a
调试结果:a =
0 1 2 3 4 5 6 7 8 9
2)源程序:
function j=juzhen(a)
a=rand(4,5);
m=a(1);ind=[1 1];
for i=1:size(a,1)
for j=1:size(a,2)
if m<a(i,j)
m=a(i,j);
ind=[i j];
end
end
end
m
ind
a(ind(1),ind(2))
调试结果:m =0.9502 ind =3 5
ans =0.9502 ans =5
3)源程序:function r = fun(n)
n=20
result = 1;
r = 0;
for ii = 1 : n
result = result * ii;
re= r + result;
end
调试结果:ans = 2.5613e+018
4)源程序:function q=xiaoqiu(s)
s=zeros(1,10);
s(1)=100;
for i=1:9
s(i+1)=s(i)/2;
end
L=s(1)+2*sum(s(2:10))
h=s(10)
调试结果:L = 299.6094 h =0.1953
5)源程序:function f=fun(x,y)
f=x^2+sin(x*y)+2*y;
调试结果:在command window界面中敲入fun(1,2),
得到ans = 5.9093
2、源程序:
clear
ax1=subplot(1,2,1);
t=0:.01*pi:2*pi;
f1=polar(t,1+cos(t));
title('心形线'); legend(f1,'r=1+cost')
ax2=subplot(1,2,2);
x1=linspace(-20,20,100);
y1=linspace(-20,20,100);
[xx1,yy1]=meshgrid(x1,y1);
zz1=xx1.*yy1;
f2=surf(xx1,yy1,zz1);
title('马鞍面'); legend(f2,'z=x^2-y^2','Location','NorthWest') xlabel('x ')
ylabel('y ')
zlabel('z')
axis([-20 20 -20 20 -400 400])
setAllowAxesRotate(rotate3d,ax1,false);
grid on
调试结果:
3、1)源程序:
t=linspace(-pi, pi,120);
x1=2*sin(t);
y1=2*cos(t);
x2=sin(t);
y2=cos(t);
x3=2*sin(t);
y3=cos(t);
plot(x1,y1,x2,y2,x3,y3)
调试结果:
2)源程序:x1=linspace(0,5,10);x2=linspace(0,5,10);
y1=exp(x1);
y2=log(x2);
y3=x1;
plot(x1,y1,x2,y2,x1,y3)
axis([0,5,0,5])
调试结果:。