中国民航大学 matlab实习报告 实验二

合集下载

MATLAB_实习报告

MATLAB_实习报告

实习报告实习题目:MATLAB的设计与绘图专业:学号:一.概述MATLAB 作为一种高级科学计算软件,是进行算法开发、数据可视化、数据分析以及数值计算的交互式应用开发环境。

世界上许许多多的科研工作者都在使用MATLAB 产品来加快他们的科研进程,缩短数据分析和算法开发的时间,研发出更加先进的产品和技术。

二。

实验目的1。

进一步熟悉MATLAB的软件环境和基本使用方法;2.巩固运用MATLAB、矩阵运算、多项式运算、字符串、单元数组、符号计算功能、绘图及句柄图形及它们的命令语句;3.学习和巩固数字图形图像处理在MATLAB中的知识和运用;4.培养我们独立自学、设计和撰写实验报告的能力;5.通过对MATLAB的进一步学习,是我们加深所学内容的认识,理解,掌握,能把所学知识运用到实际工作中;三.实习日记1.熟悉MATLAB的软件环境和基本使用方法(第一周)MATLAB的语言的基本功能和特点:它是一种高级科学计算语言,相对于传统的C、C++ 或者FORTRAN 语言,MATLAB 提供了高效快速解决各种科学计算问题的方法。

它具有数学计算、开发工具、数据的可视化、交互式编辑创建图形、集成的算法开发编程语言和环境、图形用户界面开发环境——GUIDE 、开放性、可扩展性强 、专业应用工具箱等强大的功能。

具有易于学习、使用方便、支持多种操作系统、丰富的内部函数、强大的图形和符号功能等优点。

路径设置:除 MATLAB 默认的搜索路径外,用户可以设置搜索路径。

设置方法为:选择 MATLAB 窗口中的 File | Set Path 命令,进入路径搜索对话框。

MATLAB 的搜索顺序:当在命令窗口中或者一个 M 文件中输入一个元素名 时,MATLAB 按照下面的顺序搜索该元素的意义,以元素 foo 为例:1)查找工作区中是否存在名为 foo 的变量; 2)在当前路径中查找是否存在名 foo 。

m 的文件;3)按照顺序查找搜索路径中是否存在该文件。

matlab_实习报告

matlab_实习报告

matlab_实习报告在大学的学习生涯中,实习是一个非常重要的环节,它能够让我们将理论知识与实际应用相结合,提升自己的专业技能和综合素质。

本次实习,我选择了使用 MATLAB 这个强大的工具进行实践操作,通过一段时间的学习和实践,我收获颇丰。

一、实习目的本次实习的主要目的是深入了解和掌握MATLAB 软件的使用方法,能够运用其解决实际问题,并提高自己的编程能力和逻辑思维能力。

同时,通过实际项目的操作,培养自己的团队协作精神和解决问题的能力,为今后的学习和工作打下坚实的基础。

二、实习单位及岗位介绍我实习的单位是_____,在实习期间,我主要负责利用 MATLAB 进行数据分析和算法实现的相关工作。

三、实习内容及过程(一)基础学习在实习的初期,我首先对 MATLAB 的基本语法和操作进行了系统的学习。

了解了变量的定义、数据类型、矩阵运算、函数的编写等基础知识。

通过大量的练习和实例,我逐渐熟悉了 MATLAB 的编程环境,能够熟练地编写简单的程序。

例如,在学习矩阵运算时,我通过编写程序实现了矩阵的加法、乘法、求逆等操作,深刻理解了矩阵运算在数学和工程中的重要应用。

(二)项目实践在掌握了基础知识后,我开始参与实际的项目。

其中一个项目是对一组数据进行分析和处理,以提取有用的信息。

首先,我使用MATLAB 读取数据文件,并对数据进行预处理,包括去除噪声、缺失值处理等。

然后,运用统计学方法对数据进行分析,计算均值、方差、相关性等统计量。

最后,通过绘图函数将分析结果以直观的图表形式展示出来,以便更好地理解数据的特征和趋势。

在这个过程中,我遇到了很多问题。

例如,数据的格式不一致导致读取错误,算法的复杂度过高导致运行时间过长等。

通过查阅资料、请教同事和不断地调试,我最终解决了这些问题,顺利完成了项目任务。

(三)算法实现除了数据分析,我还参与了算法的实现工作。

在一个图像识别的项目中,需要使用机器学习算法对图像进行分类。

matlab实习报告

matlab实习报告
(三)程序设计与文件操作
1、程序设计:
Matlab既就是一种语言,又就是一种编程环境。Mat1ab提供了很多方便 用户得工具,用于管理变量、输入输出数据以及生成与管理M文件。
用户可在M atlab得命令窗口键入一个命令,也可以由它定义得语言在编辑 器中编写应用程序,Matlab软件对此进行解释后,在Mat1ab环境下对它进行 处理,最后返回结果、
2、文件操作:
(1)指令驱动模式:即在MATLABM命令行窗口下用户输入单行指令时,MA T LAB立即处理这条指令,并显示结果,这就就是MATLAB命令行方式。缺点:命 令行方式程序可读性差,而且不能存储,当处理复杂问题与大量数据时很不方便。
(2)M文件模式:将MATLAB语句构成得程序存储成以m为扩展名得文件,然 后再执行该程序文件,这种工作模式称为程序文件
(二)矩阵及其运算
MATLAB,即“矩阵实验室”,它就是以矩阵为基本运算单元。因此,本书从 最基本得运算单元出发,介绍MATLAB得命令及其用法。
I、矩阵得表示:MATL AB得强大功能之一体现在能直接处理向量
或矩阵。当然首要任务就是输入待处理得向量或矩阵。不管就是任何矩阵(向 量),我们可以直接按行方式输入每个元素:同一行中得元素用逗号(,)或者用空 格符来分隔,且空格个数不限;不同得行用分号(;)分隔。所有元素处于一方括 号([])内;当矩阵就是多维(三维以上),且方括号内得元素就是维数较低得矩 阵时,会有多重得方括号。其中有:符号矩阵、大矩阵得生成、多维数组得创建、 特殊矩阵得生成(主要有零矩阵(zero).单位阵(eye)、全一阵(ones)、均匀 分布随机矩阵(rand).正态分布随机矩阵(r an d n)等)。
(5)解线性方程与函数极值线性方程:解线性方程包括:线性方程组求解、 非线性方程数值求解、常微分方程初值问题得数值解法、函数极值等。

实验二--MATLAB程序设计(含实验报告)

实验二--MATLAB程序设计(含实验报告)

实验二 MATLAB 程序设计一、 实验目的1.掌握利用if 语句实现选择结构的方法。

2.掌握利用switch 语句实现多分支选择结构的方法。

3.掌握利用for 语句实现循环结构的方法。

4.掌握利用while 语句实现循环结构的方法。

5.掌握MATLAB 函数的编写及调试方法。

二、 实验的设备及条件计算机一台(带有MATLAB7.0以上的软件环境)。

M 文件的编写:启动MATLAB 后,点击File|New|M-File ,启动MATLAB 的程序编辑及调试器(Editor/Debugger ),编辑以下程序,点击File|Save 保存程序,注意文件名最好用英文字符。

点击Debug|Run 运行程序,在命令窗口查看运行结果,程序如有错误则改正三、 实验内容1.编写求解方程02=++c bx ax 的根的函数(这个方程不一定为一元二次方程,因c b a 、、的不同取值而定),这里应根据c b a 、、的不同取值分别处理,有输入参数提示,当0~,0,0===c b a 时应提示“为恒不等式!”。

并输入几组典型值加以检验。

(提示:提示输入使用input 函数)2.输入一个百分制成绩,要求输出成绩等级A+、A 、B 、C 、D 、E 。

其中100分为A+,90分~99分为A ,80分~89分为B ,70分~79分为C ,60分~69分为D ,60分以下为E 。

要求:(1)用switch 语句实现。

(2)输入百分制成绩后要判断该成绩的合理性,对不合理的成绩应输出出错信息。

(提示:注意单元矩阵的用法)3.数论中一个有趣的题目:任意一个正整数,若为偶数,则用2除之,若为奇数,则与3相乘再加上1。

重复此过程,最终得到的结果为1。

如:2?13?10?5?16?8?4?2?16?3?10?5?16?8?4?2?1运行下面的程序,按程序提示输入n=1,2,3,5,7等数来验证这一结论。

请为关键的Matlab 语句填写上相关注释,说明其含义或功能。

MATLAB实习报告

MATLAB实习报告

MATLAB实习报告实习报告一、实习单位概述我在公司完成了为期两个月的MATLAB实习。

该公司是一家致力于研发和生产先进电子产品的科技型企业。

公司专注于电脑硬件和软件产品的开发,拥有一支强大的研发团队,同时与多个国内外知名大学与研究机构合作。

二、实习内容在实习期间,我主要负责开发一个基于MATLAB的图像处理算法,用于自动识别和分类图像中的目标物体。

具体而言,我的任务包括以下几个方面:1.数据预处理:对输入图片进行预处理,包括灰度化、降噪、图像增强等操作,以提高后续算法的准确性和鲁棒性。

2.特征提取:使用MATLAB内置的图像处理工具箱,提取目标物体的特征信息,如颜色、纹理、形状等。

3.目标识别:通过编写MATLAB脚本,实现目标物体的自动识别,使用机器学习算法训练分类器,并应用于测试集中的图像。

4.结果评估:对算法的准确性和性能进行评估,包括计算分类准确率、召回率等指标,并分析算法的优缺点和改进空间。

三、实习收获通过这次实习,我收获了许多宝贵的经验和知识。

首先,我深入学习了MATLAB的图像处理工具箱,了解了常用的图像处理算法和技术。

在实践中,我熟悉了MATLAB的基本语法和函数,并学会了如何利用MATLAB进行图像的读取、处理和保存。

其次,我掌握了图像特征提取和目标识别的方法和技巧。

通过对图像进行灰度化、降噪和图像增强等预处理步骤,我能够有效提取目标物体的特征信息,并通过机器学习算法实现目标的自动识别和分类。

最后,我学会了如何评估和改进算法的准确性和性能。

通过对测试集中的图像进行分类,并计算分类准确率和召回率等指标,我能够客观地评估算法的优劣,并提出针对性的改进方法。

四、实习总结通过这次实习,我不仅得到了理论知识的巩固和实践技能的提升,还培养了团队合作和解决问题的能力。

在与团队成员的合作中,我学会了如何与他人相互协作,分享经验与知识,并且通过互相帮助和配合,达到更好的工作效果。

同时,面对实习过程中遇到的问题和困难,我学会了如何利用各种资源和工具,积极主动地解决问题。

matlab实践报告

matlab实践报告

竭诚为您提供优质文档/双击可除matlab实践报告篇一:matlab实习报告mATLAb实习报告姓名:吴涛专业:电子信息工程班级学号:信息(2)班20XX1605010230指导教师:宋艳霞钱云实习时间:20XX,5,13至20XX,5,17目录:一:实习目的 (3)二:实习任务...........................三:实习要求...........................四:实习目的...........................五:实习体会...........................一:实习目的熟悉matlab的软件环境熟悉matlab的基本使用方法二:实习任务1.掌握matlab语言的基本语法规则及基本操作命令的使用。

2.熟悉运用matlab的数组,矩阵运算,数学运算的运算方法。

3.熟悉matlab的字符串,单元数组。

4.熟悉matlab的程序设计。

5.熟悉matlab的符号计算功能。

6.熟悉matlab的绘图及句柄图形。

7.6熟悉matlab的guI设计。

三:实习要求1.每次上机要签到,记录。

2.做的题目要在上机结束后以电子版和woRD文档形式交给指定负责人。

3.认真做一份实习总结报告。

四:实习内容共六题第1题:求a和b的和。

代码:a=[1,2,3]b=[4,5,6]s=a+b结果:s=579第2题:求a和b的差。

代码:a=[4,5,6]b=[1,2(:matlab实践报告),3]h=a-b结果:h=333第3题:求下列方程组的解。

6x1+3x2+4x3=3-2x1+5x2+7x3=-48x2-4x2-3x3=-7代码:A=[6,3,4;-2,5,7;8,-4,-3] b=[3;-4;-7]a=A\b结果:a=0.60007.0000-5.4000第4题:用for语句求三角函数表。

代码:forx=0:0.1:pi/4disp([x,sin(x),cos(x),tan(x)]),en d结果:00100.10000.09980.99500.10030.20000.19870.98010.20270.30000.29550.95530.30930.40000.38940.92110.4228篇二:matlab课程实践报告mATLAb实践》课程设计学生姓名:林淑真学号:110900824专业班级:通信工程四班指导教师:郑晓明二○一二年六月十五日《目录1.设计目的.........................................................22.题目分析.........................................................23.总体设计.........................................................34.具体设计.........................................................35.结果分析.........................................................226.心得体会.........................................................237.参考书目.........................................................238.附录 (24)1、设计目的综合运用mATLAb工具箱实现图像处理的guI程序设计。

matlab实验报告

matlab实验报告
(1)信号形式 ,利用plot命令绘制该信号在 区间内,且采样间隔为 的曲线。添加网格线,标题为“正弦信号曲线”,横坐标显示“x”,纵坐标显示“y”。
(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语言及其应用》实验报告(第二章--Matlab基本程序设计)

国家开放大学《Matlab语言及其应用》实验报告(第二章--Matlab基本程序设计)

国家开放大学《Matlab语言及其应用》实验报告
——Matlab基本程序设计
姓名:学号:
1.M文件
M文件是一个文本文件,可以用任何文本编辑器来建立和编辑,通常使用 Matlab 自带的 M 文件编辑器。

例:编写一个脚本文件将华氏温度转化为摄氏温度。

c=5
9
(f−32)
●新建一个 M 文件 f2cs.m,内容如下:
clear; % 清除当前工作空间中的变量
f=input('Please input Fahrenheit temperature:');
c=5*(f-32)/9;
fprintf('The centigrade temperature is %g\n',c);
●在命令窗口中输入 f2cs,即可执行该 M 文件。

脚本文件运行后,文件中所产生的所有变量都驻留在 Matl ab 当前工作空间,可以用 whos 查看。

2.程序控制结构
例:计算函数y=cos(x+1)+√x2+1 x≥10d的值。

clear
x=input('请输入x的值:');
if x>=10
y=cos(x+1)+sqrt(x*x+1);
end
y
请输入x的值:12
y =12.9490
3.程序调试
例:拼写错误
clear
t=0:0.1:10;
y=son(t);
例:函数错误
调试器的使用:MatLab的调试器和M文件编辑器集成一体。

matlab实习报告

matlab实习报告

Matlab上机实习实习报告
上周我们进行了本学期第二个实习部分-Matlab上机实习。

本次实习一星期,共做了四个部分。

分别是:一,Matlab基本操作;二,数值数组及其运算;三,Matlab图形绘制基础;四,SIMULINK仿真基础。

下面分别介绍个部分实习内容。

一、实习目的
1、掌握Matlab的启动和退出;
2、熟悉Matlab的命令窗口;
3、熟悉其他常用窗口和工具栏;
二、对本软件有了初步认识后,接下来就开始指导书要求做一些训练了
(一)Simulink 图
(二)程序
1、主程序
global kp;
global ki;
global i;
i=1;
result=fminsearch('liyanan3',[2 1])%[2 1]为寻优的初值
2、调用程序
function ss=optm(x)
global kp;
global ki;
global i;
kp=x(1);
ki=x(2);
i=i+1;%记录寻优次数
[tt,xx,yy]=sim('liyanan1',40,[]);
yylong=length(yy);
ss=yy(yylong);
(三)调试结果
1、T=2时的调试结果result =
3.5566 0.2869
Scop2
Scop1
2、T=3时的调试结果result =
2.4292 0.1888 Scop2
Scop1。

MATLAB实训报告_2

MATLAB实训报告_2

M AT L A B语言编程实训报告题目MATLAB在信号与系统中的应用年级2008 专业电子信息工程目录第一章绪论 (2)1.1 本设计课题目的及意义, 重点解决的问题 (2)1.1.1 MATLAB在信号与系统中应用的目的及意义 (2)1.1.2 重点解决的问题 (2)1.2 课题的社会和技术背景 (2)1.3 实现的具体功能 (3)第二章课题的基本概念和原理 (4)2.1 MATLAB的概念 (4)2.2 信号与系统的概念 (4)2.3 信号与系统分析的基本内容与方法 (5)2.4 离散系统的基本概念 (6)2.5连续系统模型及表示 (6)第三章系统设计和实现 (7)3.1采用的软件及开发平台 (7)3.1.1系统设计软件 (7)3.1.2开发平台和编程环境 (7)3.2系统的详细设计 (8)3.2.1建模 (8)3.2.2解决方法 (9)3.3系统设计的亮点 (9)第四章结束语 (10)参考文献 (10)第一章绪论1.1 本设计课题目的及意义, 重点解决的问题1.1.1 MATLAB在信号与系统中应用的目的及意义MATLAB在信号与系统中应用能够让学生熟悉MATLAB软件平台、工具箱、高效的数值计算及符号计算功能。

熟悉MATLAB软件的信号处理编程方法和结果的可视化。

了解数字信号处理的计算机仿真方法。

进一步加深对信号与系统的基本原理、方法及应用的理解。

MATLAB软件具有强大的数值分析和计算结果可视化的功能.运用MATLAB软件,通过一个应用实例,将信号与系统课程的理论与实践教学有机地结合,有效地解决了教学中的难点问题,说明了MATLAB软件在信号与系统课程的实践教学中具有重要的实际意义1.1.2 重点解决的问题matlab的基本使用方法, matlab的预定义函数, 如何用matlab绘图, 如何在matlab中编程, 基于matlab的矩阵运算、符号运算、数值分析等。

1.2 课题的社会和技术背景长期以来, “信号与系统”课程一直采用黑板式的单一的教学方式, 学生仅依靠做习题来巩固和理解教学内容, 对课程中大量的应用性较强的内容不能实际动手设计、调试、分析, 严重影响和制约了教学效果。

MATLAB实验报告2

MATLAB实验报告2

MATLAB实验报告2实验(2)Matlab 基础操作实验:符号运算求解微积分与线性代数1. 求下列函数极限(Find the limits of the following functions )(1) 0sin lim x xx →syms x r ;r=sin(x)/x;limit(r,x,0)ans =1(2) 1lim(1)x x x →∞+syms x r ;r=(1+1/x)^x;limit(r,x,inf)ans =exp(1) (3) lim (1)xx a x →+∞+syms x r a ;r=(1+a/x)^x;limit(r,x,inf)ans =exp(a)(4) 1lim (1)xx x →-∞+syms x r ;r=(1+1/x)^x;limit(r,x,-inf)ans =exp(1) (5)201cos lim x xx →-syms x r ;r=(1-cos(x))/(x^2);limit(r,x,0)ans =实验(2)Matlab 基础操作实验:符号运算求解微积分与线性代数1/2(6) lim nsyms n r ;r=n^(1/n);limit(r,n,inf)ans =1(7)0lim x x π+→ syms x n r ;r=(cos(sqrt(x)))^(n/x);limit(r,x,0,'right')ans =1/exp(n/2)2.求下列函数的导数或偏导数( Find the derivatives of the following functions)(1)()()f x g x ==syms x a b c f g ;f=sqrt(a*x^2+b*x+c);g=sqrt(exp(x^2)+x*sin(x));a1=diff(f,x)a2=diff(g,x)a1 =(b + 2*a*x)/(2*(a*x^2 + b*x + c)^(1/2))a2 =(sin(x) + 2*x*exp(x^2) + x*cos(x))/(2*(exp(x^2) + x*sin(x))^(1/2))(2)3()ln()f x x =syms x f ;f=log(x^3);实验(2)Matlab基础操作实验:符号运算求解微积分与线性代数diff(f,x)ans =3/x(3)Find the 3rd derivative of f(x).()f x=2xxe-syms x f;f=x*exp(-x^2);diff(f,x,3)ans =(24*x^2)/exp(x^2) - 6/exp(x^2) - (8*x^4)/exp(x^2) (4) )Find23222(235)x x y yx-+-?syms x y f;f=x^3-2*x^2*y^2+3*y-5; diff(f,x,2)ans =6*x - 4*y^2(5)Find2322(235)x x y yx y-+-.syms x y f;f=x^3-2*x^2*y^2+3*y-5;fx=diff(f,x);fxy=diff(fx,y)fxy =-8*x*y3. 求下列函数的不定积分或定积分(Find indefinite integrals or definite integrals of the functions)(1)1(sin2cos3)xx x e dxx--++syms x s;s=sin(x)-2*cos(3*x)+1/x+exp(-x); int(s,x)ans =实验(2)Matlab 基础操作实验:符号运算求解微积分与线性代数log(x) - (2*sin(3*x))/3 - cos(x) - 1/exp(x)(2)sin()x x e e dx ?syms x s ;s=exp(x)*sin(exp(x));int(s,x)ans =-cos(exp(x))(3)2syms x s ;s=x^2/(sqrt(x^6)+4);int(s,x)Warning: Explicit integral could not be found.ans =int(x^2/((x^6)^(1/2) + 4), x)(4)cos3cos5x xdx ?syms x s ;s=cos(3*x)*cos(5*x);int(s,x) ans =sin(2*x)/4 + sin(8*x)/16(5) syms x s a ;s=(sqrt(x^2-a^2))/x;int(s,x)ans =实验(2)Matlab 基础操作实验:符号运算求解微积分与线性代数(x^2 - a^2)^(1/2) - log(((-a^2)^(1/2) + (x^2 -a^2)^(1/2))/x)*(-a^2)^(1/2)(6) sin ax e bxdx ?syms x s a b ;s=exp(a*x)*sin(b*x);int(s,x)ans =-(exp(a*x)*(b*cos(b*x) - a*sin(b*x)))/(a^2 + b^2) (7) 0πsyms x s ;s=sqrt(sin(x)-(sin(x))^3);int(s,x,0,pi)Warning: Explicit integral could not be found.ans =int((sin(x) - sin(x)^3)^(1/2), x = 0..pi) (8) 211dx x +∞? syms x s ;s=1/(x^2);int(s,x,1,inf)14. 解下列方程(Solve the equations.)x =syms x ;solve(sqrt(1-x^2)-x,x)ans =2^(1/2)/2实验(2)Matlab基础操作实验:符号运算求解微积分与线性代数(2)1024x y zx y zx y z++=-+=--=-syms x y z;s=solve(x+y+z-10,x-y+z,2*x-y-z+4,x,y,z) s.x,s.y,s.z s =x: [1x1 sym]y: [1x1 sym]z: [1x1 sym]2ans =5ans =3(3)240330sin0x xy zx yzy z++=+-=+=syms x y z;s=solve(x^2+4*x*y+z,x+3*y*z-3,y+sin(z),x,y,z) s.x,s.y,s.z s =x: [1x1 sym] y: [1x1 sym] z: [1x1 sym]实验(2)Matlab 基础操作实验:符号运算求解微积分与线性代数ans =-6.4840528860501102693511549985753ans =-0.071995758403084664887655346911159ans =-43.9102390493223146672631539563435. 求解下列常微分方程(Solve the following ordinary differential equations.)(1)ln()xy y xy y '=-syms x y ;s=dsolve('x*Dy=y*log(x*y)-y','x')s =1/xexp(exp(C8 + log(x)))/x (2)20(1)5dv t dt v ?+==? syms v t ;s=dsolve('Dv+2*t=0','v(1)=5','t')s =6 - t^2(3)()0y a b y aby '''-++=syms a b x y ;s=dsolve('D2y-(a+b)*Dy+a*b*y=0','x')s =实验(2)Matlab 基础操作实验:符号运算求解微积分与线性代数C13*exp(a*x) + C14*exp(b*x)6. 用MATLAB 验证(Use MALTAB to prove the following identities)(1) 22sin cos 1x y +=syms x ;s=solve((sin(x))^2+cos((x))^2-1,x)s =C_求解此方程所得的解是全体实数因此该等式在实数域上成立。

Matlab实验报告二

Matlab实验报告二

Matlab 程序设计实验报告一.实验目的1. 掌握MATLAB 数据对象的特点和运算规则。

2. 掌握MATLAB 中建立矩阵的方法和矩阵的处理方法。

二.实验内容23100.7780414565532503269.5454 3.14A -⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥-⎣⎦,0.434328.9421B ⎡⎤=⎢⎥-⎣⎦ 1. 将矩阵A 的某个元素用机器零代替>> A=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14]23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 5.0000 0 32.00006.0000 -9.5400 54.0000 3.1400>> A(3,2)ans =5>> A(3,2)=epsA =23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 0.0000 0 32.00006.0000 -9.5400 54.0000 3.1400>>2.将A转化为2x8的矩阵,转化成1x16、16x1的向量>> reshape(A,2,8)ans =23.0000 32.0000 10.0000 0.0000 -0.7780 0 0 32.000041.0000 6.0000 -45.0000 -9.5400 65.0000 54.0000 5.0000 3.1400 >> reshape(A,1,16)ans =Columns 1 through 923.0000 41.0000 32.0000 6.0000 10.0000 -45.0000 0.0000 -9.5400 -0.7780Columns 10 through 16>> reshape(A,16,1)ans =23.000041.000032.00006.000010.0000-45.00000.0000-9.5400-0.778065.000054.00005.000032.00003.140065.0000 0 54.0000 0 5.0000 32.0000 3.14003.自定义一个字符串,求其ASCII>> zifu='abcd 12345ABCD'zifu =abcd 12345ABCD>> abs(zifu)ans =97 98 99 100 32 49 50 51 52 53 65 66 67 68>>4.产生和A一样大小的0矩阵和1矩阵>> A=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14]A =23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 5.0000 0 32.00006.0000 -9.5400 54.0000 3.1400>> ndims(A)ans =2>> size(A)ans =4 4>> zeros(size(A))ans =0 0 0 00 0 0 00 0 0 00 0 0 0>> ones(size(A))ans =1 1 1 11 1 1 11 1 1 11 1 1 1>>5.从A提取主对角线元素,构成对角阵>> diag(diag(A))ans =23.0000 0 0 00 -45.0000 0 00 0 0 00 0 0 3.1400>>6.产生均值为3,方差为1的500个正态分布的随机序列,写入表达式>> y=3+sqrt(1)*rand(10,50)y =Columns 1 through 93.9501 3.6154 3.0579 3.0153 3.8381 3.1934 3.4966 3.7271 3.79483.2311 3.7919 3.3529 3.7468 3.0196 3.6822 3.8998 3.3093 3.95683.6068 3.9218 3.8132 3.4451 3.6813 3.3028 3.8216 3.8385 3.52263.4860 3.7382 3.0099 3.9318 3.3795 3.5417 3.6449 3.5681 3.88013.8913 3.1763 3.1389 3.4660 3.8318 3.1509 3.8180 3.3704 3.17303.7621 3.4057 3.2028 3.4186 3.5028 3.6979 3.6602 3.7027 3.97973.4565 3.9355 3.1987 3.8462 3.7095 3.3784 3.3420 3.5466 3.27143.0185 3.9169 3.6038 3.5252 3.4289 3.8600 3.2897 3.4449 3.25233.8214 3.4103 3.2722 3.2026 3.3046 3.8537 3.3412 3.6946 3.87573.4447 3.8936 3.1988 3.6721 3.1897 3.5936 3.5341 3.6213 3.7373Columns 10 through 183.1365 3.5828 3.2091 3.4154 3.2140 3.6833 3.4514 3.6085 3.08413.0118 3.4235 3.3798 3.3050 3.6435 3.2126 3.0439 3.0158 3.45443.8939 3.5155 3.7833 3.8744 3.3200 3.8392 3.0272 3.0164 3.44183.1991 3.3340 3.6808 3.0150 3.9601 3.6288 3.3127 3.1901 3.35333.2987 3.4329 3.4611 3.7680 3.7266 3.1338 3.0129 3.5869 3.15363.6614 3.2259 3.5678 3.9708 3.4120 3.2071 3.3840 3.0576 3.67563.2844 3.5798 3.7942 3.9901 3.7446 3.6072 3.6831 3.3676 3.69923.4692 3.7604 3.0592 3.7889 3.2679 3.6299 3.0928 3.6315 3.72753.0648 3.5298 3.6029 3.4387 3.4399 3.3705 3.0353 3.7176 3.47843.9883 3.6405 3.0503 3.4983 3.9334 3.5751 3.6124 3.6927 3.5548Columns 19 through 273.1210 3.2319 3.4398 3.9342 3.1370 3.4225 3.2974 3.3759 3.19393.4508 3.2393 3.3400 3.2644 3.8188 3.8560 3.0492 3.0099 3.90483.7159 3.0498 3.3142 3.1603 3.4302 3.4902 3.6932 3.4199 3.56923.8928 3.0784 3.3651 3.8729 3.8903 3.8159 3.6501 3.7537 3.63183.2731 3.6408 3.3932 3.2379 3.7349 3.4608 3.9830 3.7939 3.23443.2548 3.1909 3.5915 3.6458 3.6873 3.4574 3.5527 3.9200 3.54883.8656 3.8439 3.1197 3.9669 3.3461 3.4507 3.4001 3.8447 3.93163.2324 3.1739 3.0381 3.6649 3.1660 3.4122 3.1988 3.3678 3.33523.8049 3.1708 3.4586 3.8704 3.1556 3.9016 3.6252 3.6208 3.65553.9084 3.9943 3.8699 3.0099 3.1911 3.0056 3.7334 3.7313 3.3919Columns 28 through 363.6273 3.7165 3.1146 3.3603 3.4319 3.7327 3.3567 3.6700 3.90903.6991 3.5113 3.6649 3.5485 3.6343 3.4222 3.4983 3.2009 3.59623.3972 3.7764 3.3654 3.2618 3.8030 3.9614 3.4344 3.2731 3.32903.4136 3.4893 3.1400 3.5973 3.0839 3.0721 3.5625 3.6262 3.47823.6552 3.1859 3.5668 3.0493 3.9455 3.5534 3.6166 3.5369 3.59723.8376 3.7006 3.8230 3.5711 3.9159 3.2920 3.1133 3.0595 3.16143.3716 3.9827 3.6739 3.7009 3.6020 3.8580 3.8983 3.0890 3.82953.4253 3.8066 3.9994 3.9623 3.2536 3.3358 3.7546 3.2713 3.95613.5947 3.7036 3.9616 3.7505 3.8735 3.6802 3.7911 3.4091 3.59553.5657 3.4850 3.0589 3.7400 3.5134 3.0534 3.8150 3.4740 3.0287Columns 37 through 453.8121 3.9566 3.1879 3.2460 3.2859 3.7241 3.8289 3.3663 3.22193.6101 3.1472 3.4906 3.5874 3.3941 3.2816 3.1663 3.3025 3.70373.7015 3.8699 3.4093 3.5061 3.5030 3.2618 3.3939 3.8518 3.52213.0922 3.7694 3.4635 3.4648 3.7220 3.7085 3.5208 3.7595 3.93293.4249 3.4442 3.6109 3.5414 3.3062 3.7839 3.7181 3.9498 3.71343.3756 3.6206 3.0712 3.9423 3.1122 3.9862 3.5692 3.5579 3.22803.1662 3.9517 3.3143 3.3418 3.4433 3.4733 3.4608 3.0142 3.44963.8332 3.6400 3.6084 3.4018 3.4668 3.9028 3.4453 3.5962 3.17223.8386 3.2473 3.1750 3.3077 3.0147 3.4511 3.0877 3.8162 3.96883.4516 3.3527 3.6210 3.4116 3.6641 3.8045 3.4435 3.9771 3.3557Columns 46 through 503.0490 3.2093 3.1998 3.5979 3.99953.7553 3.4551 3.0495 3.9492 3.21203.8948 3.0811 3.5667 3.2888 3.49843.2861 3.8511 3.1219 3.8888 3.29053.2512 3.5620 3.5221 3.1016 3.67283.9327 3.3193 3.1171 3.0653 3.95803.1310 3.3749 3.7699 3.2343 3.76663.9408 3.8678 3.3751 3.9331 3.66613.7019 3.3722 3.8234 3.0631 3.13093.8477 3.0737 3.0466 3.2642 3.09547.求A和B的上三角矩阵、下三角矩阵、逆矩阵、行列式的值、秩、范数、条件数、迹A =23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 5.0000 0 32.00006.0000 -9.5400 54.0000 3.1400>> B=[0.43,43,2;-8.9,4,21]B =0.4300 43.0000 2.0000-8.9000 4.0000 21.0000triu(A)ans =23.0000 10.0000 -0.7780 00 -45.0000 65.0000 5.00000 0 0 32.00000 0 0 3.1400tril(A)ans =23.0000 0 0 041.0000 -45.0000 0 032.0000 5.0000 0 06.0000 -9.5400 54.0000 3.1400>> inv(A)ans =0.0307 0.0090 -0.0004 -0.01040.0297 -0.0210 0.0008 0.02570.0039 -0.0044 -0.0017 0.0238-0.0353 -0.0057 0.0315 0.0064>> Pinv(B)Warning: Could not find an exact (case-sensitive) match for 'Pinv'. C:\MATLAB701\toolbox\matlab\matfun\pinv.m is a case-insensitive match and will be used instead. You can improve the performance of your code by using exact name matches and we therefore recommend that you update your usage accordingly. Alternatively, you can disable this warning using warning('off','MATLAB:dispatcher:InexactMatch').ans =0.0022 -0.01750.0234 -0.0017-0.0035 0.0405>> det(A)ans =1.8965e+006>> rank(A)ans =4>> norm(A,1)ans =119.7780>> cond(A,1)ans =11.9307>> trace(A)ans =-18.8600>>8.A=[34,NaN,Inf,-Inf,-pi,eps,0],对A求all、any、isnan、isinf、isfinite函数的值A=[34,NaN,Inf,-Inf,-pi,eps,0]A =34.0000 NaN Inf -Inf -3.1416 0.0000 0>> all(A)ans =>> any(A)ans =1>> isnan(A)ans =0 1 0 0 0 0 0>> isinf(A)ans =0 0 1 1 0 0 0>> isfinite(A)ans =1 0 0 0 1 1 1>>9.建立结构矩阵,并索引其中的元素a.x1=10; a.x2='abc'; a(2).x1=12; a(2).x2='wang';10.建立单元矩阵,并索引其中的元素。

matlab实习报告二5篇

matlab实习报告二5篇

matlab实习报告二5篇第一篇:matlab实习报告二MATLAB实习报告(2)实验二 MATLAB矩阵分析与处理王夏一、实验目的1、掌握生成特殊矩阵的方法。

2、掌握矩阵分析的方法。

3、用矩阵求逆发解线性方程组。

二、实验内容1、设有分块矩阵A=[E3×3 R3×2;O2×3 S2×2],其中E、R、O、S分别为单位矩阵、随机矩阵、零矩阵和对角阵,试通过数值计算验证A²=[E R+RS;O S²]。

程序清单:E=eye(3);R=rand(3,2);O=zeros(2,3);S=diag([4,5]);A=[E R;O S] ;A2=A^2;C=[E R+R*S;O S^2];length(find(A2==C))==25 运行结果:ans =12、产生5阶希尔伯特矩阵H和5阶帕斯卡矩阵P,且求其行列式的值Hh和Hp以及他们的条件数Th和Tp,判断哪个矩阵的性能更好,为什么?程序清单:format rat H=hilb(5);format short P=pascal(5);Hh=det(H);Hp=det(P);Th=cond(A);Tp=cond(P);运行结果:Hh =3.7493e-012 Hp =1 Th =5.5228 Tp =8.5175e+003 实验收获:会建立希尔伯特矩阵和帕斯卡矩阵,知道怎么求矩阵行列式的值以及条件数。

希尔伯特矩阵的性能更好,条件数越接近1的矩阵性能越好。

3、建立一个5×5矩阵,求它的行列式值、迹、秩和范数。

程序清单:A=[1:5;1:0.1:1.5;2 5 7 3 9;2:6;3:0.4:4.6]Ha=det(A);Ra=rank(A);Ta=trace(A);Na=norm(A);运行结果:Ha =1.4421e-031 Ra = 3 Ta =18.7000 Na =19.49664、已知向量A,求A的特征值及特征向量,并分析其数学意义。

matlab实习报告3

matlab实习报告3

matlab实验报告实验二MATLAB的基本数值运算一.实验目的MATLAB具有非常强大的数值计算能力,对各种常量(包括一般常量和固定常量)、各种变量(包括数值变量、字符变量、单元形变量和结构变量)熟悉其用法,向量及其运算,矩阵及其运算,数组及其运算,多项式及其各种运算,线性方程组的求解,数值统计的基本函数及其应用,简单插值函数,简单优化函数,微分方程数值解的基本函数等。

二.实验要求理解常量、变量、向量、矩阵、多项式等概念,掌握向量在MATLAB中的表示,熟练掌握矩阵及其运算,数组及其运算,多项式运算,线性方程组求解,了解数值统计的基本函数,简单插值函数,简单优化函数,微分方程数值解的基本函数等。

三.实验内容实验一第二节:1.简单的矩阵建立与矩阵元素:>> a=[5 6 8 7;7 1 5 4;2 3 8 3]a =5 6 8 77 1 5 42 3 8 3>> b=cos(a)b =0.2837 0.9602 -0.1455 0.75390.7539 0.5403 0.2837 -0.6536-0.4161 -0.9900 -0.1455 -0.9900>> c=a+0.3*(1+b/2)c =5.34256.4440 8.27827.41317.4131 1.3810 5.3425 4.20202.23763.1515 8.2782 3.1515>> d=ones(5)+eye(5)d =2 1 1 1 11 2 1 1 11 12 1 11 1 12 11 1 1 1 2>> x=[-3.2 sqrt(5) (2+5+1)*2/4]x =-3.2000 2.2361 4.0000>> x(6)=abs(x(2))x =-3.2000 2.2361 4.0000 0 0 2.2361>> e=[a,b]e =5.00006.0000 8.00007.0000 0.2837 0.9602 -0.1455 0.75397.0000 1.0000 5.0000 4.0000 0.7539 0.5403 0.2837 -0.65362.00003.0000 8.0000 3.0000 -0.4161 -0.9900 -0.1455 -0.9900 >> f=e(:,[3,6])f =8.0000 0.96025.0000 0.54038.0000 -0.99002.常量变量与表达式:>> g='I can do it'g =I can do it>> A='Hello'A =Hello>> B=2*AB =144 202 216 216 2223.了解命令窗口的常用命令及标点符号的使用;4.学习并掌握内存变量的编辑、保存和清除;5.算术运算:(1)加减法运算:>> h=[2 5 ;7 9];g=ones(2);k=h+gk =3 68 10>> l=k+2l =5 810 12(2)乘除法运算:>> m=k*lm =75 96140 184>> n=pi*mn =235.6194 301.5929439.8230 578.0530>> C=k.*lC =15 4880 120>> D=l.*kD =15 4880 120数组的乘法满足交换率6.乘幂运算:>> a=[5 6 8 7;7 1 5 4;2 3 8 3]a =5 6 8 77 1 5 42 3 8 3 >> a.^2ans =25 36 64 4949 1 25 164 9 64 9 7.转置运算:>> a=[1 2 3;4 2 6;7 4 9]a =1 2 34 2 67 4 9>> c=a'c =1 4 72 2 43 6 98.关系运算与逻辑运算:(1)关系运算:>> a=[1 2 3;4 2 6;7 4 9];b=[5 2 1;4 2 6;3 4 8];c=a>bc =0 0 10 0 01 0 1(2)逻辑运算:>> a=[1 2 3;4 2 6;7 4 9];b=[5 2 1;4 2 6;3 4 8];c=a&bc =1 1 11 1 11 1 1>> a=[1 2 3;4 0 6;7 4 9];b=[5 2 1;4 2 6;3 4 8];c=a&bc =1 1 11 0 11 1 1>> a=[1 2 3;4 0 6;7 4 9];b=[5 2 1;4 2 6;3 4 8];c=a|bc =1 1 11 1 11 1 1>> a=[1 2 3;4 2 6;7 4 9];c=~ac =0 0 00 0 00 0 09.建立特殊数组:>> a=[1 2;4 2]a =1 24 2>> b=[a,eye(2,3);ones(3,2),rand(3)]b =1.00002.0000 1.0000 0 04.0000 2.0000 0 1.0000 01.0000 1.0000 0.7922 0.0357 0.67871.0000 1.0000 0.9595 0.8491 0.75771.0000 1.0000 0.6557 0.9340 0.7431>> a=[1 2 5;4 4 2]a =1 2 54 4 2>> a(:,3)=[]a =1 24 4实验四:一.矩阵,向量的建立与修改由文件生成和保存矩阵?myfileA =1 8 9 5 6 815 8 23 6 2 87 5 9 15 8 321 4 14 15 24 131 4 5 8 9 65 8 7 96 2?whoYour variables are:A?load txtfile.txt?whoYour variables are:A txtfile?save matfile?clear?who?load matfile?whoYour variables are:A txtfile?txtfiletxtfile =1 8 9 5 6 815 8 23 6 2 87 5 9 15 8 321 4 14 15 24 131 4 5 8 9 65 8 7 96 2由函数生成矩阵?eye(3)ans =1 0 00 1 00 0 1?eye(size(A))??? Undefined function or variable 'A'.?ones(3)ans =1 1 11 1 11 1 1?zeros(3)ans =0 0 00 0 00 0 0?rand(4)ans =0.9501 0.8913 0.8214 0.92180.2311 0.7621 0.4447 0.73820.6068 0.4565 0.6154 0.17630.4860 0.0185 0.7919 0.4057?randn(5)ans =-0.4326 1.1909 -0.1867 0.1139 0.2944 -1.6656 1.1892 0.7258 1.0668 -1.33620.1253 -0.0376 -0.5883 0.0593 0.71430.2877 0.3273 2.1832 -0.0956 1.6236-1.1465 0.1746 -0.1364 -0.8323 -0.6918ans =8 1 63 5 74 9 2符号矩阵的创建?exam=sym('[1,x/a,sin(x);y/x,1+1/y,tan(x/y);1=0,3+3,4*r]') exam =[ 1, x/a, sin(x)][ y/x, 1+1/y, tan(x/y)][ 1=0, 3+3, 4*r]矩阵的修改?A=rand(5)A =0.9355 0.3529 0.1987 0.7468 0.84620.9169 0.8132 0.6038 0.4451 0.52520.4103 0.0099 0.2722 0.9318 0.20260.8936 0.1389 0.1988 0.4660 0.67210.0579 0.2028 0.0153 0.4186 0.8381 ?A([1 3],:)=[]A =0.9169 0.8132 0.6038 0.4451 0.52520.8936 0.1389 0.1988 0.4660 0.67210.0579 0.2028 0.0153 0.4186 0.8381 ?A(:,[2 4])=[]A =0.9169 0.6038 0.52520.8936 0.1988 0.67210.0579 0.0153 0.8381?A(2:3,2:3)=eye(2)A =0.9169 0.6038 0.52520.8936 1.0000 00.0579 0 1.0000二.矩阵向量的基本运算矩阵的运算?myfileA =1 8 9 5 6 815 8 23 6 2 87 5 9 15 8 321 4 14 15 24 131 4 5 8 9 65 8 7 96 2txtfile =1 8 9 5 6 815 8 23 6 2 87 5 9 15 8 321 4 14 15 24 131 4 5 8 9 65 8 7 96 2 ?A+txtfileans =2 16 18 10 12 1630 16 46 12 4 1614 10 18 30 16 642 8 28 30 48 262 8 10 16 18 1210 16 14 18 12 4 ?A-txtfileans =0 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 0 ?A*txtfileans =335 225 430 383 316 216 464 395 676 646 500 359 483 257 530 516 574 372 583 530 828 873 900 607 303 181 345 350 363 225 379 215 462 379 384 282 ?2*A+3*txtfileans =5 40 45 25 30 4075 40 115 30 10 4035 25 45 75 40 15105 20 70 75 120 655 20 25 40 45 3025 40 35 45 30 10 ?A'ans =1 15 7 21 1 58 8 5 4 4 89 23 9 14 5 75 6 15 15 8 96 2 8 24 9 68 8 3 13 6 2?txtfile/Aans =1.0000 0 0 0 0 00 1.0000 0 0 0 0-0.0000 0.0000 1.0000 -0.0000 0.0000 -0.00000.0000 -0.0000 0.0000 1.0000 -0.0000 -0.0000-0.0000 0.0000 -0.0000 -0.0000 1.0000 0.0000-0.0000 0.0000 -0.0000 -0.0000 0.0000 1.0000?det(A)ans =172040?inv(A)ans =0.3072 -0.1053 0.1918 0.0939 -0.5346 -0.10130.2260 -0.0823 0.0307 0.0287 -0.3019 0.0988-0.5492 0.2256 -0.3248 -0.1105 0.7671 0.19830.2514 -0.0955 0.2864 0.0190 -0.3329 -0.1780-0.5344 0.1472 -0.3939 -0.0534 0.7305 0.29480.7221 -0.2091 0.4271 0.1121 -0.8341 -0.4190?rank(A)ans =6三.解线性方程组奇次线性方程组求解?A=[1 -2 3 -4;0 1 -1 1;-1 0 -1 2;1 -3 4 -5];?a=null(A)a =0.5890 -0.5655-0.0117 0.57720.5655 0.58900.5772 0.0117?a=null(A,'r')a =-1 21 -11 00 1四.矩阵特征值与特征向量求矩阵特征值与特征向量?AA =1 -23 -40 1 -1 1-1 0 -1 21 -3 4 -5?[v,d]=eig(A)v =0.5843 0.7746 -0.2420 0.7674-0.1711 -0.2582 0.4132 -0.5231-0.2420 0.2582 -0.5843 -0.27890.7555 0.5164 -0.6552 0.2443d =-4.8284 0 0 00 -0.0000 0 00 0 0.8284 00 0 0 0.0000?a=eig(A)a =-4.8284-0.00000.82840.0000?B=[1 2 3 4 ;5 2 4 3;4 2 5 3;5 1 4 2]B =1 2 3 45 2 4 34 25 35 1 4 2?[vl,dl]=eig(A,B)vl =0.1553 1.0000 0.0694 0.2301-1.0000 -0.5758 1.0000 -0.6150-0.2213 -0.1515 -0.2466 -1.00000.7085 0.4242 -0.3920 -0.3850dl =-4.1242 0 0 00 -0.0000 0 00 0 4.6242 00 0 0 0.0000四、实验总结在这次实习中,实习了MATLAB基本数学运算和线性代数实验,通过这次实习,掌握了MATLAB的一些基本数学运算,并且熟悉了矩阵的各种运算,线性方程组的求解及特征值特征向量的求解。

(2021年整理)Matlab实验报告

(2021年整理)Matlab实验报告

Matlab实验报告编辑整理:尊敬的读者朋友们:这里是精品文档编辑中心,本文档内容是由我和我的同事精心编辑整理后发布的,发布之前我们对文中内容进行仔细校对,但是难免会有疏漏的地方,但是任然希望(Matlab实验报告)的内容能够给您的工作和学习带来便利。

同时也真诚的希望收到您的建议和反馈,这将是我们进步的源泉,前进的动力。

本文可编辑可修改,如果觉得对您有帮助请收藏以便随时查阅,最后祝您生活愉快业绩进步,以下为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/23.两个函数之间的操作求和(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上机实习实习报告十八周我们进行了本学期的第二次实习——Matlab上机实习。

本次实习一星期,共实习了五个部分。

分别是:一,Matlab基本操作;二,数值数组及其运算;三,Matlab图形绘制基础;四:SIMULINK仿真基础;五:句柄图形和图形用户界面制作。

下面分别介绍各部分内容:实习内容一Matlab基本操作1.实习目的:●掌握Matlab的启动和退出;●熟悉的命令窗口;●熟悉其它常用窗口和工具栏;2.输入一段程序: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’)输出结果:3.了解文件中的一些功能和指令.例如:clear 清除当前工作区中所有变量clc 清除指令窗内容clf 清除图形窗口cd 设置当前工作目录exit,quit 退出MATLAB实习内容二数值数组及其运算1.实习目的●掌握一维数组的创建和寻访●掌握二维数组的创建和寻访●掌握MATLAB的矩阵和数组的运算●熟悉MATLAB关系操作和逻辑操作2.要求在闭区间[0,2pi]上产生50个等距采样的一维数组A,试用两种不同的指令实现。

要访问1到5元素如何实现;寻访7到最后一个元素如何实现;寻访第2.6.8个元素如何实现;寻访大于2的元素如何实现;给第3.5.9个元素赋值100如何实现:(1)>> A=linspace(0,2*pi,50);>> A=0:(50/49):50;>> AA =Columns 1 through 90 1.0204 2.0408 3.0612 4.0816 5.1020 6.1224 7.1429 8.1633Columns 10 through 189.1837 10.2041 11.2245 12.2449 13.2653 14.2857 15.3061 16.3265 17.3469Columns 19 through 2718.3673 19.3878 20.4082 21.4286 22.4490 23.4694 24.4898 25.5102 26.5306Columns 28 through 3627.5510 28.5714 29.5918 30.6122 31.6327 32.6531 33.6735 34.6939 35.7143Columns 37 through 4536.7347 37.7551 38.7755 39.7959 40.8163 41.8367 42.8571 43.8776 44.8980Columns 46 through 5045.9184 46.9388 47.9592 48.9796 50.0000(2)>> whos A>> A(1:5)ans =0 1.0204 2.0408 3.0612 4.0816(3)>> A(7:50)ans =Columns 1 through 86.12247.14298.16339.1837 10.2041 11.2245 12.2449 13.2653Columns 9 through 1614.2857 15.3061 16.3265 17.3469 18.3673 19.3878 20.4082 21.4286Columns 17 through 2422.4490 23.4694 24.4898 25.5102 26.5306 27.5510 28.5714 29.5918Columns 25 through 3230.6122 31.6327 32.6531 33.6735 34.6939 35.7143 36.7347 37.7551Columns 33 through 4038.7755 39.7959 40.8163 41.8367 42.8571 43.8776 44.8980 45.9184Columns 41 through 4446.9388 47.9592 48.9796 50.0000(4)>> A([2 6 8])ans =1.0204 5.1020 7.1429(5)>> poison=find(A>2)poison =Columns 1 through 143 4 5 6 7 8 9 10 11 12 13 14 15 16Columns 15 through 2817 18 19 20 21 22 23 24 25 26 27 28 29 30Columns 29 through 4231 32 33 34 35 36 37 38 39 40 41 42 43 44Columns 43 through 4845 46 47 48 49 50 (6)>> A([3 5 9])=100 (7)>> A=[1 2 3;4 5 6] A =1 2 3 4 5 63. 在指令窗口输入A=[1 2 3;4 5 6] 观察结果:分别用“全坐标”“单坐标”来寻访第2行第1列的元素,并说明二维数组元素的存储顺序。

matlab实习实验报告

matlab实习实验报告

实验报告: 分支语句一、 实验目的: 1.掌握使用分支语句;2.掌握自上而下程序设计技术方法.二、 实验内容及要求:1.实验内容:1).编写 MATLAB 语句计算 y(t)的值⎩⎨⎧<+≥+-=0530 53)(22t t t t t y 已知 t 从-5到 5 每隔0.5取一次值。

运用循环和选择语句进行计算。

2).用向量算法解决练习 1, 比较这两个方案的耗时。

三、 设计思路:1. 用循环和选择语句进行计算:1).定义自变量t :t=-5:0.5:5;2).用循环语句实现对自变量的遍历。

3).用选择语句实现对自变量的判断,选择。

4).将选择语句置入循环语句中,则实现在遍历中对数据的选择,从而实现程序的功能。

2. 用向量法实现:1).定义自变量t :t=-5:0.5:5;2).用 b=t>=0 语句,将t>=0得数据选择出,再通过向量运算y(b)=-3*t(b).^2 + 5; 得出结果。

3).用取反运算,选择出剩下的数据,在进行向量运算,得出结果。

四、 实验程序和结果1.实验程序clear%clctic;t=-5:0.5:5;for ii=1:size(t,2)columns.if(t(ii)<0)y(ii) = 3*t(ii)^2+5;elsey(ii)= -3*t(ii)^2+5;endendfigure(1);plot(t,y);title('Plot of y(t) and its derivative----(1)');xlabel('x');ylabel('y');grid on;toc; %Read the stopwatch timer,%prints the number of seconds required for the operation. clear%clctic;t=[-5:0.5:5];b=t>=0;y(b)=-3*t(b).^2 + 5;%b=t<0;y(~b)=3*t(~b).^2 + 5;figure(2);plot(t,y);title('Plot of y(t) and its derivative----(2)');xlabel('x');ylabel('y');grid on;toc;2.实验结果:>> clear>> y_tElapsed time is 0.998095 seconds.Elapsed time is 0.338708 seconds.>>五、实验总结:本次实验分支语句的使用和自上而下的程序设计技术。

matlab实习报告总结

matlab实习报告总结

MATLAB实习报告总结介绍本文是对我在MATLAB实习过程中的总结与反思。

通过这次实习,我对MATLAB的使用有了更深入的了解,并且对实际问题的解决能力也有了提升。

实习目标在实习开始之前,我设定了以下几个目标: 1. 熟悉MATLAB的基本语法和常用函数; 2. 学习如何使用MATLAB解决实际问题; 3. 提升自己的编程能力和解决问题的思维能力。

实习过程第一阶段:熟悉MATLAB基本语法在实习开始之初,我花了一段时间来学习MATLAB的基本语法。

通过阅读官方文档和参考书籍,我了解了MATLAB的变量定义、矩阵运算、函数定义等基本概念和用法。

我利用空闲时间编写了一些简单的代码来巩固这些知识。

第二阶段:学习MATLAB的常用函数在熟悉基本语法之后,我开始学习MATLAB的常用函数。

我学习了MATLAB 中的统计分析函数、数据可视化函数、图像处理函数等。

通过学习这些函数,我可以更高效地处理和分析数据,同时也能够更好地展示我的结果。

第三阶段:实际问题的解决在掌握了MATLAB的基本语法和常用函数后,我开始尝试使用MATLAB解决一些实际问题。

我选择了一个与我的专业相关的问题,通过MATLAB对数据进行分析和建模。

在这个过程中,我遇到了一些困难和挑战,但通过查阅资料和请教同学,我最终成功解决了这些问题。

第四阶段:编程能力和思维能力的提升通过这次实习,我的编程能力和解决问题的思维能力得到了提升。

在实际问题的解决过程中,我学会了如何分析问题、设计算法、编写代码并调试。

通过与同学的讨论和交流,我也学会了如何思考和解决问题的不同角度。

总结与反思这次MATLAB实习让我受益匪浅。

通过实际操作,我深入了解了MATLAB的使用方法和技巧。

我学会了如何使用MATLAB解决实际问题,并且在这个过程中提升了自己的编程能力和解决问题的思维能力。

然而,我也发现了自己的不足之处。

在实习过程中,我有时会陷入“死记硬背”的状态,只注重语法和函数的使用,而不去思考问题背后的本质,这导致我在解决一些复杂问题时遇到了困难。

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

实验二数据和函数的可视化
实验目的:
1.熟悉matlab软件中二维绘图和三维绘图的方法
2.学习使用matlab图形窗功能
实验内容:
基本命令
>> x=linspace(0,2*pi,100);
>> y=sin(x);
>> plot(x,y);
问题1:
1.当运行x=[1 5 3 7;3 6 8 4;9 6 1 5];y=[2 5 7 4;6 8 4 1;8 0 4 2];plot(x,y)时的结果:
2.当运行x=[1 5 3 7];y=[2 5 7 4;6 8 4 1;8 0 4 2];plot(x,y)时显示:
3. 运行x=[1 5 3 7;3 6 8 4;9 6 1 5];y=[2 5 7 4;6 8 4 1;8 0 4 2];m=cat(2,x,y);plot(m);显示:
例2:t=0:pi/20:pi;
plot(t,t.*cos(t),'-.r*');
hold on
plot(exp(t/100).*sin(t-pi/2),'--m0')
plot(sin(t*pi),':bs') %误把mo写成了m0 ^o^,导致如下结果:
??? Error using ==> plot
Error in color/linetype argument.
Error in ==> C:\MA TLAB6p5\work\Untitled2.m
On line 4 ==> plot(exp(t/100).*sin(t-pi/2),'--m0')
正确输入:
t=0:pi/20:pi;
plot(t,t.*cos(t),'-.r*');
hold on
plot(exp(t/100).*sin(t-pi/2),'--mo') plot(sin(t*pi),':bs')
如果没有hold on:
t=0:pi/20:pi;
plot(t,t.*cos(t),'-.r*');
plot(exp(t/100).*sin(t-pi/2),'--mo')
plot(sin(t*pi),':bs')
所以hold on 起到继续在当前图形上画图的作用
图形修饰与控制
axis square %将图形设置为正方形
axis equal %x,y轴单位刻度相等
title (‘字符串’)%图形标题、
axis([xmin,xmax,ymin,ymax])%x轴范围在xmin-xmax之间,y轴范围在ymin-ymax之间xlabel(‘字符串’)%x轴标注
ylabel(‘字符串’)%在(x,y)处标注说明文字
text(x,y,’字符串’)%在(x,y)处标注说明文字
grid on % 加网格线
grid off % 取消网格线
hold on %保持当前图形
hold off %解除hold on 命令
legend(‘first’,’second’,n) % 对一个坐标系上的两幅图形做出图例注解
subplot(m,n,p)将当前窗口分成m行n列区域,并指定在p区绘图
fill(x,y,’b’)%将x(1),y(1),(x(2),y(2)),…,(x(n),y(n)),(x(1),y(1))围成的封闭图形填充为蓝色
例三:subplot(2,2,1);x=0:pi/60:2*pi;plot(x,cos(x));
subplot(2,2,2);x=0:pi/60:2*pi;plot(x,sin(x));
subplot(2,1,2);plot(x,sin(x),':b',x,cos(x),'-r')
legend('sin(x)','cos(x)',1)
例四:
clf; %清除图片clear;
x=0:pi/60:2*pi;
y=sin(x);
x1=0:pi/60:pi/2;
y1=sin(x1);
plot(x,y,'-r')
hold on;
fill([x1,pi/2],[y1,0],'b')
问题2:将上面最后一句改为fill(x1,y1,'b'),fill([pi/2,x1,3*pi/2],[0,y1,0],'r') 效果
对数图形
格式:loglog(x,y) %对x轴pi,y轴的刻度用常用对数(以10为底)semilogx(x,y) %对x轴的刻度常用对数值,而y轴为线性刻度
semilogy(x,y)%对y轴的刻度常用对数值,而x轴为线性刻度
例五:
clear;
x=logspace(-1,2);
loglog(x,10*exp(x),'-s')
grid on
三维曲线图
格式:plot3(x,y,z,s)
说明:当x,y,z均为同维向量时,则plot3描点x(i),y(i),z(i)依次相连的空间曲线,若x,y 均为同维矩阵,x,y,z每一组相应列向量为坐标画出一条曲线。

S为’color-linestyle-marker’
例六绘制螺旋线
t=0:pi/60:10*pi;
x=sin(t);
y=cos(t);
plot3(x,y,t,'*-b')% *表示用*画出,-表示用线连接
%还有s 方格,:点;
grid on
三维网格图
格式:mesh(x,y,z,c) %画出颜色由C指定的三维网格图
说明:若x,y均为向量,n=length(x),z必须满足[m,n]=size(z),则空间中的点(x(j),y(i),z(i,j))
为所画曲面网线的交点,x对应于z点的列,y对应于z的行;若x,y,z均为同维矩阵,则空间中的点(x(i,j(,y(i,j),z(i,j))为所画曲面的网线的交点:矩阵C指定网线的颜色。

例七绘出带有底座的马鞍面
x=-8:8;
y=-8:8;
[x,y]=meshgrid(x,y);
z=(x.^2/4^2-y.^2/5^2);
mesh(x,y,z)
三维曲面图
格式:surf(x,y,z,c) %画出颜色由C指定的三维曲面。

说明:surf同mesh命令用法和使用格式相同,不同之处在于绘得的图形是一个彩色曲面而不是彩色网格。

C缺省时数据z同时为曲面高度,也是彩色数据。

R=8.31;
n=2;
p=(1:20);
v=(1:20);
[P,V]=meshgrid(p,v);
T=P.*V/n/R;
surf(P,V,T)
view(45,30)
练习与思考
1.输入x=[2,4,1,6,8];plot(x,'*-b'),显示:
2.在一幅图上画出两个周期的正弦曲线和余弦曲线,画出坐标轴,加上各种图注,并在正
弦曲线(0到2*pi)和横轴之间填充红色。

x=0:0.1:4*pi
plot(sin(x))
hold on
plot(cos(x))。

相关文档
最新文档