北科大matlab第6次实验报告
MATLAb实验报告第六次
![MATLAb实验报告第六次](https://img.taocdn.com/s3/m/624ca7e14afe04a1b071de0b.png)
60 k< k<k80,即 k80.
也就是说,木工,电工,油漆工的日工资分别为 k元, k元,k元,其中 k80.
为了简便起见,可取k= 72,于是木工,电工,油漆工的日工资分别为62元, 64元, 72元.
试建立数学模型确定该交通网络未知部分的具体流量。
3.详细设计:
由网络流量假设,所给问题满足如下线性方程组
由Matlab软件的编辑器构造M函数文件:jt. m如下:
A5 =[0 1 -1 1 0 0 0 0 0 0
0 0 0 1 1 0 0 0 0 0
0 0 0 0 0 -1 1 0 0 0
1 1 0 0 0 0 0 0 0 0
其对应的非齐次方程组为
赋值给自由未知量(x5,x8)为(0,0)得非齐次方程组的特解
4:实验结果:
方程组的通解为X = k1 η1 + k2 η2 + x,其中k1,k2为任意常数,X的每一个分量即为交通网络未知部分的具体流量. X有无穷多解.
增广矩阵的秩rank(B5)= 8
增广矩阵阶梯形最简形式为:
B6 =[1 0 0 0 1 0 0 0 0 0 800
0 1 0 0 -1 0 0 0 0 0 0
0 0 1 0 0 0 0 0 0 0 200
0 0 0 1 1 0 0 0 0 0 500
0 0 0 0 0 1 0 1 0 0 800
0 0 0 0 0 0 1 1 0 0 1000
B5 =[A5 ';b]';
Rank(A5);
Rank(B5);
MATLAB实验报告
![MATLAB实验报告](https://img.taocdn.com/s3/m/2459908db04e852458fb770bf78a6529647d35b4.png)
MATLAB实验报告一、实验目的本次 MATLAB 实验旨在深入了解和掌握 MATLAB 软件的基本操作和应用,通过实际编程和数据处理,提高解决问题的能力,培养编程思维和逻辑分析能力。
二、实验环境本次实验使用的是 MATLAB R2020a 版本,运行在 Windows 10 操作系统上。
计算机配置为英特尔酷睿 i5 处理器,8GB 内存。
三、实验内容(一)矩阵运算1、矩阵的创建使用直接输入、函数生成和从外部文件导入等方式创建矩阵。
例如,通过`1 2 3; 4 5 6; 7 8 9` 直接输入创建一个 3 行 3 列的矩阵;使用`ones(3,3)`函数创建一个 3 行 3 列元素全为 1 的矩阵。
2、矩阵的基本运算包括矩阵的加减乘除、求逆、转置等。
例如,对于两个相同维度的矩阵`A` 和`B` ,可以进行加法运算`C = A + B` 。
3、矩阵的特征值和特征向量计算通过`eig` 函数计算矩阵的特征值和特征向量,加深对线性代数知识的理解和应用。
(二)函数编写1、自定义函数使用`function` 关键字定义自己的函数,例如编写一个计算两个数之和的函数`function s = add(a,b) s = a + b; end` 。
2、函数的调用在主程序中调用自定义函数,并传递参数进行计算。
3、函数的参数传递了解值传递和引用传递的区别,以及如何根据实际需求选择合适的参数传递方式。
(三)绘图功能1、二维图形绘制使用`plot` 函数绘制简单的折线图、曲线等,如`x = 0:01:2pi; y = sin(x); plot(x,y)`绘制正弦曲线。
2、图形的修饰通过设置坐标轴范围、标题、标签、线条颜色和样式等属性,使图形更加清晰和美观。
3、三维图形绘制尝试使用`mesh` 、`surf` 等函数绘制三维图形,如绘制一个球面`x,y,z = sphere(50); surf(x,y,z)`。
(四)数据处理与分析1、数据的读取和写入使用`load` 和`save` 函数从外部文件读取数据和将数据保存到文件中。
MATLAB实验报告
![MATLAB实验报告](https://img.taocdn.com/s3/m/42c5135d26284b73f242336c1eb91a37f011325d.png)
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 第六次数学实验报告 (精华版)](https://img.taocdn.com/s3/m/2f3e19d9240c844769eaee63.png)
《数学实验》报告实验名称Matlab高等数学相关运算学院专业班级姓名学号2013年6月一、 【实验目的】1. 学会用Matlab 进行常微分方程的求解、随机试验和统计作图;2. 掌握相关运算处理的基本指令和参数设置 二、 【实验任务】 P168 习题24 求解微分方程yxx y cos sin ,=。
P168 习题27用数值方法求解析下列微分方程,用不同颜色和线形将y 和y ’画在同一个图形窗口里:t y y y 2-1-t ,,,=+初始时间:0t =0;终止时间:π=f t ;初始条件:2.0|1.0|00,====t t y y 。
P190 习题15描绘以下数组的频数直方图:6.8, 29.6,33.6,35.7, 36.9, 45.2, 54.8, 65.8, 43.4, 53.8, 63.7, 69.9, 70.7, 79.5, 97.9, 139.4, 157.0 P190 习题16若样本为85,86,78,90,96,82,80,74 求样本均值、标准差、中位数、极差和方差。
三、 【实验程序】 习题24:>>dsolve('Dy=x*sin(x)/cos(y)','x') 习题27:function xdot=exf(t,x) u=1-2*t;xdot=[0,1;1,-t]*x+[0,1]'*u;clf;t0=0;tf=pi;x0t=[0.1;0.2];[t,x]=ode23('exf',[t0,tf],x0t);y=x(:,1)Dy=x(:,2)plot(t,y,'r-', t,Dy,'b*')legend('y','Dy')xlabel('t轴')习题15:clf;load A.txt;figure(1)hist(A,5)figure(2)hist(A,10)figure(3)hist(A,20)习题16:B=[85 86 78 90 96 82 80 74];disp(' 样本均值标准差中位数极差方差') C=[mean(B),std(B),median(B),range(B),var(B)]四、【实验结果】习题24:习题27:习题15:习题16:五、【实验总结】通过这最后一次实验,我学习了怎么用Matlab作常微分方程的求解、概率论与数理统计的相关计算,感受到了Matlab软件的强大与方便。
(最新版)MATLAB实验报告
![(最新版)MATLAB实验报告](https://img.taocdn.com/s3/m/315fed3afe00bed5b9f3f90f76c66137ee064ff5.png)
(最新版)MATLAB实验报告实验一典型环节的MATLAB仿真一、实验目的1.熟悉MATLAB桌面和命令窗口,初步了解SIMULINK功能模块的使用方法。
2.通过观察典型环节在单位阶跃信号作用下的动态特性,加深对各典型环节响应曲线的理解。
3.定性了解各参数变化对典型环节动态特性的影响。
二、SIMULINK的使用MATLAB中SIMULINK是一个用来对动态系统进行建模、仿真和分析的软件包。
利用SIMULINK功能模块可以快速的建立控制系统的模型,进行仿真和调试。
1.运行MATLAB软件,在命令窗口栏“>>”提示符下键入simulink命令,按Enter 键或在工具栏单击按钮,即可进入如图1-1所示的SIMULINK仿真环境下。
2.选择File菜单下New下的Model命令,新建一个simulink 仿真环境常规模板。
3.在simulink仿真环境下,创建所需要的系统。
以图1-2所示的系统为例,说明基本设计步骤如下:1)进入线性系统模块库,构建传递函数。
点击simulink下的“Continuous”,再将右边窗口中“Transfer Fen”的图标用左键拖至新建的“untitled”窗口。
2)改变模块参数。
在simulink仿真环境“untitled”窗口中双击该图标,即可改变传递函数。
其中方括号内的数字分别为传递函数的分子、分母各次幂由高到低的系数,数字之间用空格隔开;设置完成后,选择OK,即完成该模块的设置。
3)建立其它传递函数模块。
按照上述方法,在不同的simulink 的模块库中,建立系统所需的传递函数模块。
例:比例环节用“Math”右边窗口“Gain”的图标。
4)选取阶跃信号输入函数。
用鼠标点击simulink下的“Source”,将右边窗口中“Step”图标用左键拖至新建的“untitled”窗口,形成一个阶跃函数输入模块。
5)选择输出方式。
用鼠标点击simulink下的“Sinks”,就进入输出方式模块库,通常选用“Scope”的示波器图标,将其用左键拖至新建的“untitled”窗口。
matlab实验报告
![matlab实验报告](https://img.taocdn.com/s3/m/b5d04b94ec3a87c24028c4eb.png)
MATLAB实验报告姓名:学号: 0905130129专业班级:通信工程1301班学院:信息科学与工程学院指导老师:陈科文、支国明、张金焕、周扬目录实验一熟悉MATLAB环境 (3)[实验目的] (3)[实验原理] (3)[实验内容] (4)[实验结果] (5)[实验遇到的问题及解决] (11)实验二数值数组创建、应用及可视化 (12)[实验目的] (12)[实验原理] (12)[实验内容] (13)[实验结果] (15)[实验遇到的问题及解决] (20)实验三字符串数组的使用、简单脚本文件和函数的编写 (21)[实验目的] (21)[实验原理] (21)[实验内容] (24)[实验结果] (26)[实验遇到的问题及解决] (35)实验四数据可视化方法 (35)[实验目的] (35)[实验原理] (35)[实验内容] (37)[实验结果] (39)[实验遇到的问题及解决] (46)总结 (46)实验一熟悉MATLAB环境[实验目的]1 .熟悉MATLAB 主界面,并学会简单的菜单操作。
2 .学会简单的矩阵输入与信号输入。
3 .掌握部分绘图函数。
[实验原理]MATLAB 是以复杂矩阵作为基本编程单元的一种程序设计语言。
它提供了各种矩阵的运算与操作,并有较强的绘图功能。
用户第一次使用MATLAB 时,建议首先在屏幕上键入demo 命令,它将启动MATLAB 的演试程序,用户可在此演示程序中领略MATLAB 所提供的强大的运算与绘图功能。
也可以键入help 进行进一步了解。
MATLAB 启动界面如图所示:操作界面主要的介绍如下:指令窗(Command Window ),在该窗可键入各种送给MATLAB 运作的指令、函数、表达式,并显示除图形外的所以运算结果。
历史指令窗(Command History ),该窗记录已经运行过的指令、函数、表达式;允许用户对它们进行选择复制、重运行,以及产生M 文件。
工作空间浏览器(Workspace Browser ),该窗口罗列出MATLAB 工作空间中所有的变量名、大小、字节数;并且在该窗中,可对变量进行观察、编辑、提取和保存。
matlab实验六报告
![matlab实验六报告](https://img.taocdn.com/s3/m/aa893d1452d380eb62946dec.png)
实验六MATLAB数据可视化实验者:祝松年级:机设092 学号:09405701002一、实验目的掌握 MATLAB 二维、三维图形绘制,掌握图形属性的设置和图形修饰;掌握图像文件的读取和显示。
二、实验内容(1) 二维图形绘制。
(2) 三维曲线和三维曲面绘制。
(3) 图像文件的读取和显示。
三、实验步骤1 .二维图形绘制(1) 二维图形绘制主要使用函数 plot 。
>>clear all;>> x=linspace(0,2*pi,100);>> y1=sin(x);>> plot(x,y)>> hold on % 保持原有的图形>> y2=cos(x)>> plot(x,y)注: hold on 用于保持图形窗口中原有的图形, hold off 解除保持。
(2) 函数 plot 的参数也可以是矩阵。
>>close all % 关闭所有图形窗口>>x=linspace(0,2*pi,100);>> y1=sin(x);>> y2=cos(x);>> A=[y1 ; y2]'; % 把矩阵转置>> B=[x ; x]'>> plot(B,A)(3) 选用绘图线形和颜色。
>> close all % 关闭所有图形窗口>> plot(x,y1,'g+',x,y2, 'r:')>> grid on % 添加网格线(4) 添加文字标注。
>> title(' 正弦曲线和余弦曲线 ')>> ylabel(' 幅度 ')>> xlabel(' 时间 ')>> legend('sin(x)', 'cos(x)')>> gtext('\leftarrowsinx') % 可用鼠标选择标注的位置,%\leftarrow 产生左箭头,‘\’为转义符(5) 修改坐标轴范围。
MATLAB第6次实验报告
![MATLAB第6次实验报告](https://img.taocdn.com/s3/m/8f14de076edb6f1aff001f39.png)
实验报告一、实验目的通过这些题目的训练,需要掌握静态数据成员以及静态成员函数的使用方法,友元函数的使用,以及掌握构造函数和析构函数,拷贝构造函数、内联成员函数、聚集的使用,熟练的使用类解决问题。
二、实验环境编译器:Visual c++6.0操作系统:Windows XP三、实验题目基础题1:题目:定义一个Cat类,拥有静态数据成员HowManyCats,记录Cat的个体数目;静态成员函数GetHowMany(),存取HowManyCats。
设计程序测试这个类,体会静态数据成员和静态成员函数的用法。
分析:定义一个名为Cat类,静态数据成员HowManyCats,记录Cat的个体数目,静态成员函数GetHowMany(),输入一只Cat的名字就让HowManyCats加1,调用静态成员函数GetHowMany()输出数目即可。
核心代码:#include<iostream>using namespace std;class Cat//定义一个Cat类{public:Cat(char Name)//定义一个构造函数,形参为Cat的名字{CatName=Name;HowManyCats++;}static void GetHowMany(){cout<<"the number of cats: "<<HowManyCats<<endl;}private:char CatName;static int HowManyCats;};int Cat::HowManyCats=0;//静态数据成员声明void main(){char _name;int W=0;{ cout<<"please input the name of one cat(一个字母):"<<endl;cin>>_name;W=1;Cat c(_name);cout<<"if you want to contiue to input cat,please input 0 or else input 1."<<endl;cin>>W;}while(W==0);Cat::GetHowMany();//静态成员函数调用}程序运行及截图:总结:掌握静态成员函数的定义与实现,以及静态数据成员的声明、定义和初始化方法,学会使用类解决实际问题。
matlab实践报告
![matlab实践报告](https://img.taocdn.com/s3/m/b1eaf39549649b6648d747b1.png)
竭诚为您提供优质文档/双击可除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大作业](https://img.taocdn.com/s3/m/3dfba910c850ad02de8041b1.png)
《数学实验》报告实验名称MATLAB在研究物体振动方面的应用学院专业班级姓名学号2015年 1月一、【实验目的】物体振动这样一个看似简单但又包含着很多复杂计算的运动中,在人为的计算时是很难精确的实现,而通过MATLAB可以处理诸多科学中的许多问题,利用它来研究物理学中的机械振动,不仅特别方便还非常有效。
二、【实验任务】本列举振动的一些实例,用matlab语言编制计算机程序进行仿真以达到研究简谐振动以及振动的合成,振动的计算以及受迫振动。
三、【实验程序】(一)简谐振动介绍最简单和最基本的振动是简谐振动.任何复杂的振动,都可以看成为许多简谐振动的合成.1.特点质点作简谐振动的条件是:在任何时候所受到的力与质点离开平衡位置的位移成正比,其指向与位移相反,始终指向平衡位置.所受的力与位移的关系表示为(1)式中为正的常数.对于弹簧振子,就是弹簧劲度系数2.运动的微分方程及其解根据牛顿第二定律,作简谐振动的质点的微分方程写成即(2)式中。
如下面的(3)和(4)所示,是简谐振动的圆频率。
微分方程(2)的解是(3)或(4)式(7.3)也可以表为复数形式(5)但要约定取其实数部分.利用三角公式,很容易导出A ,和B,C之间的关系即(6)3.速度和加速度作简谐振动的质点,它的速度和加速度很容易得到.只要将(7.3)对时间分别求导一次和求导两次即可,(7)(8)式(1)、(2)、(3)、(4)、(5)都是判别一个系统是否作简道振动的依椐.4.圆频率、周期和频率之间的关系,,(9),,三者不是独立的,只要知道其中一个,就可以由(7.9)求出其余两个。
它们是由振动系统的固有性质决定,常称为固有圆频率,固有周期和固有频率.5.振幅和初周相(3)中和是两个积分常数,可由初始条件决定.将初始条件:“,,”代入(3)和(7),得(10)解得(11)求解质点作简谐振动的具体运动情况,也就是要确定(7.3)中的,,三个值.其中和由初始条件决定,因此一般来说,首先必须确定初始值和,而根据(7.10)或(7.11)求出和值.至于(或或),它是由系统固有性质决定的,与初始情况无关.例如对于弹簧振子,,完全由弹簧劲度系数和物体质量所决定.弹簧的大(即所谓硬的弹簧),振动的圆频率也就大。
MATLAB实验报告
![MATLAB实验报告](https://img.taocdn.com/s3/m/a9a321dc0408763231126edb6f1aff00bed57077.png)
MATLAB实验报告实验⼀ MATLAB 环境的熟悉与基本运算⼀、实验⽬的及要求1.熟悉MATLAB 的开发环境;2.掌握MATLAB 的⼀些常⽤命令;3.掌握矩阵、变量、表达式的输⼊⽅法及各种基本运算。
⼆、实验内容1、熟悉MATLAB 的开发环境: ① MATLAB 的各种窗⼝:命令窗⼝、命令历史窗⼝、⼯作空间窗⼝、当前路径窗⼝。
②路径的设置:建⽴⾃⼰的⽂件夹,加⼊到MATLAB 路径中,并保存。
? 设置当前路径,以⽅便⽂件管理。
2、学习使⽤clc 、clear,了解其功能与作⽤。
3、矩阵运算:已知:A=[1 2;3 4]; B=[5 5;7 8]; 求:A*B 、A 、*B,并⽐较结果。
4、使⽤冒号选出指定元素:已知:A=[1 2 3;4 5 6;7 8 9]; 求:A 中第3列前2个元素;A 中所有列第2,3⾏的元素; 5、在MATLAB 的命令窗⼝计算: 1))2sin(π2) 5.4)4.05589(÷?+ 6、关系及逻辑运算1)已知:a=[5:1:15]; b=[1 2 8 8 7 10 12 11 13 14 15],求: y=a==b ,并分析结果 2)已知:X=[0 1;1 0]; Y=[0 0;1 0],求: x&y+x>y ,并分析结果 7、⽂件操作1)将0到1000的所有整数,写⼊到D 盘下的data 、txt ⽂件 2)读⼊D 盘下的data 、txt ⽂件,并赋给变量num 8、符号运算1)对表达式f=x 3-1 进⾏因式分解2)对表达式f=(2x 2*(x+3)-10)*t ,分别将⾃变量x 与t 的同类项合并 3)求3(1)xdz z +?三、实验报告要求完成实验内容的3、4、5、6、7、8,写出相应的程序、结果实验⼆ MATLAB 语⾔的程序设计⼀、实验⽬的1、熟悉 MATLAB 程序编辑与设计环境2、掌握各种编程语句语法规则及程序设计⽅法3、函数⽂件的编写与设计4、了解与熟悉变量传递与赋值⼆、实验内容1.编写程序,计算1+3+5+7+…+(2n+1)的值(⽤input 语句输⼊n 值)。
北科大_matlab第六次作业_常微分方程
![北科大_matlab第六次作业_常微分方程](https://img.taocdn.com/s3/m/ccb149d35fbfc77da269b171.png)
《数学实验》报告实验名称常微分方程学院计算机与通信工程学院专业班级姓名学号2015年6月一、 【实验目的】1、学习常微分方程符号解法和数值解法;2、学习频数直方图的作法。
二、 【实验任务】1.P168 T24 求解微分方程yxx y cos sin '=。
2.P168 T27用数值方法求解下列微分方程,用不同颜色和线型将y 和y ’画在同一个图形窗口里:t y ty y 21'''-=-+初始时间:00=t ;终止时间:π=f t ;初始条件:1.0|0==t y ,2.0|'0==t y 。
3.P190 T15描绘以下数组的频数直方图:6.8,29.6,33.6,35.7,36.9,45.2,54.8,65.8,43.4,53.8,63.7,69.9,70.7,79.5,97.9,139.4,157.0三、 【实验程序】1.P168 T24y=dsolve('Dy=x*sin(x)/cos(y)','x')2.P168 T27先将高阶微分方程转化为一阶微分方程,令''''',2121x y x y x y x =⇒===,即原微分方程化为:⎩⎨⎧-++-==tx tx x x x 21''12221 写成矩阵形式为:)21(10t -110)21(10t -110'''2121t t x x x x x -⎥⎦⎤⎢⎣⎡+⎥⎦⎤⎢⎣⎡=-⎥⎦⎤⎢⎣⎡+⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡=xu t xdot t u ⎥⎦⎤⎢⎣⎡+⎥⎦⎤⎢⎣⎡-=-=10110,21x 放入函数exf.m 中,命令如下:[t,x]=ode23('exf',[t0,tf],x0t)(1)将导数表达式的右端写成exf.m 函数文件:function [xdot ] = exf( t,x ) u=1-2*t;xdot=[0 1;1 -t]*x+[0 1]'*u; end(2)主程序如下: clft0=0,tf=pi,x0t=[0.1,0.2]; [t,x]=ode23('exf',[t0,tf],x0t); y=x(:,1);y1=x(:,2); plot(t,y,'r-',t,y1,'b--'); legend('y','y 的一阶导数');3.P190 T15A=[6.8,29.6,33.6,35.7,36.9,45.2,54.8,65.8,43.4,53.8,63.7,69.9,70.7,79.5,97.9,139.4,157.0] hist(A,12)四、 【实验结果】1.P168 T24 y =asin(C3 + sin(x) - x*cos(x))2.P168T273.P190T15五、【实验总结】1、学会了常微分方程符号解法和数值解法;2、学会了频数直方图的作法。
matlab实验内容答案
![matlab实验内容答案](https://img.taocdn.com/s3/m/5cf6fc5ac77da26924c5b012.png)
m a t l a b实验内容答案(总16页) -本页仅作为预览文档封面,使用时请删除本页-实验报告说明:matlab 课程实验需撰写8个实验报告,每个实验报告内容写每次实验内容中标号呈黑体大号字显示的题目。
第一次实验内容:实验一 MATLAB 运算基础一、实验目的1.熟悉启动和退出MATLAB 的方法。
2.熟悉MATLAB 命令窗口的组成。
3.掌握建立矩阵的方法。
4.掌握MATLAB 各种表达式的书写规则以及常用函数的使用。
二、实验内容1.先求下列表达式的值,然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1)22sin8511z e ︒=+ (2)12ln(2z x =,其中2120.455i +⎡⎤=⎢⎥-⎣⎦(3)0.30.33sin(0.3), 3.0, 2.9, 2.8,,2.8,2.9,3.02a ae e z a a --=+=--- 提示:利用冒号表达式生成a 向量,求各点的函数值时用点乘运算。
(4)2220141122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪--≤<⎩,其中t =0::提示:用逻辑表达式求分段函数值。
2.已知12344347873657A -⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,131203327B -⎡⎤⎢⎥=⎢⎥⎢⎥-⎣⎦求下列表达式的值:(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 和B12345678910111213141516171819202122232425A ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦, 30161769023497041311B ⎡⎤⎢⎥-⎢⎥⎢⎥=-⎢⎥⎢⎥⎢⎥⎣⎦(1) 求它们的乘积C 。
(2) 将矩阵C 的右下角3×2子矩阵赋给D(3) 查看MATLAB 工作空间使用情况。
MATLAB实验报告
![MATLAB实验报告](https://img.taocdn.com/s3/m/fcbf9edc360cba1aa811dadd.png)
MATLAB 实验报告一、实验名称实验6 程序设计(1)二、实验目的:熟悉和掌握MA TLAB 的m 文件程序和程序设计。
三、实验内容:1.for 循环结构%m 文件的使用:新建一个m 文件将下面的语句复制进去,保存成m 文件,假如文件名是“xunhuan ”,则只需在matlab 窗口中输入“xunhuan ”然后回车即可执行其中的内容,关于m 文件的使用方法下同for i=1:10; %i 的数值是从1到10,默认间隔是1 x(i)=i; %矩阵x 第i 个数的数值赋值i end; %结束循环x %输出矩阵x 的值2.while 循环结构Fibonaci 数组的元素满足Fibonacci 规则:);,2,1k (,a a a 1k k 2k =+=++且1a a 21==。
现要求该数组中第一个大于10000的元素。
a(1)=1;a(2)=1;i=2; %赋初始值while a(i)<=10000 %设定语句结束条件a(i+1)=a(i-1)+a(i); %利用给出的前后数的关系赋值 i=i+1; %循环控制变量每次自加1 end; %结束循环i ,a(i) %输出符合条件的i 值和a(i)值3.if-else-end 结构cost=10;number=12; %赋初始值 if number>8 %如果+条件 sums=number*0.95*cost; %计算公式end,sums %结束循环,输出sums 的值4.编写一个脚本文件和函数文件并调用,说明其调用方式以及各自的工作空间的不同。
脚本文件stat1.m%脚本文件%求阵列x 的平均值 %[m,n]=size(x); if m==1 m=n; ends=sum(x);mean=s/m;函数文件stat2.mfunction[mean1]=stat2(x) %函数文件%求阵列x的平均值%调用格式为%[mean]=stat2(x)%[m,n]=size(x);if m==1m=n;ends=sum(x);mean=s/m然后在,MATLAB下执行这两个文件,在MA TLAB中输入>>clear all>>x=rand(4,4)+2;>>stat1执行后检查基本工作空间的变量情况:>>whosName Size Bytes Classm 1x1 8 double arraymean 1x4 32 double arrayn 1x1 8 double arrays 1x4 32 double arrayx 4x4 128 double arrayGrand total is 26 elements using 208 bytes这说明,在脚本文件中产生的变量都保存在基本工作空间。
MATLAB实验报告实验六
![MATLAB实验报告实验六](https://img.taocdn.com/s3/m/36b57aa668dc5022aaea998fcc22bcd126ff4234.png)
MATLAB实验报告实验六实验六1.x=linspace(0,2*pi,101);y=(0.5+3*sin(x)./(1+x.^2));plot(x,y)2.x=-pi:pi/100:pi;y1=x.^2;y2=cos(2*x);y3=y1.*y2;plot(x,y1,'r:',x,y2,'g-',x,y3,'b--')x=-pi:pi/100:pi;y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(1,3,1);plot(x,y1,'r:');title('y1=x.^2'); subplot(1,3,2);plot(x,y2,'g-');title('y2=cos(2*x)'); subplot(1,3,3);plot(x,y1,'b--');title('y3=y1.*y2');x=-pi:pi/100:pi;y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1,'r:',x,y2,'g-',x,y3,'b--'); subplot(2,2,2);bar(x,y1,'r');title('y1=x^2');subplot(2,2,3);bar(x,y2,'g');title('y2=cos(2*x)');subplot(2,2,4);bar(x,y3,'b');title('y3=y1.*y2');x=-pi:pi/100:pi;y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1,'r:',x,y2,'g-',x,y3,'b--'); subplot(2,2,2); stairs(x,y1,'r:');title('y1=x^2');subplot(2,2,3);stairs(x,y2,'g-');title('y2=cos(2*x)');subplot(2,2,4);stairs(x,y3,'b--');title('y3=y1.*y2');x=-pi:pi/100:pi;y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1,'r:',x,y2,'g-',x,y3,'b--'); subplot(2,2,2); stem(x,y1,'r:');title('y1=x^2');subplot(2,2,3);stem(x,y2,'g-');title('y2=cos(2*x)');subplot(2,2,4);stem(x,y3,'b--');title('y3=y1.*y2');x=-pi:pi/100:pi;y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1,'r:',x,y2,'g-',x,y3,'b--'); subplot(2,2,2);fill(x,y1,'r:');title('y1=x^2');subplot(2,2,3);fill(x,y2,'g-');title('y2=cos(2*x)');subplot(2,2,4);fill(x,y3,'b--');title('y3=y1.*y2');3.x=-5:0.01:5;y=(x+sqrt(pi))./(exp(2)).*(x<=0)+0.5*log(x+sqrt(1+x.*x)).*(x>0); plot(x,y4.theta=0:pi/100:2*pi; a=input('请输⼊a:');b=input('请输⼊b:');n=input('请输⼊n:'); rho=a*sin(b+n*theta); polar(theta,rho,'g')a=b=n=2a=b=2,n=4a=n=2,b=4a=4,b=n=2采⽤控制变量法的办法,固定两个参数,变动第三个参数观察输出图象的变化。
(2021年整理)Matlab实验报告
![(2021年整理)Matlab实验报告](https://img.taocdn.com/s3/m/8cf3f79969eae009591becdb.png)
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_数学实验报告1~6次(全)
![北科大Matlab_数学实验报告1~6次(全)](https://img.taocdn.com/s3/m/3e83eb37580216fc700afd15.png)
《数学实验》报告实验名称 Matlab 基础知识学院专业班级姓名学号2014年 6月一、【实验目的】1.认识熟悉Matlab这一软件,并在此基础上学会基本操作。
2.掌握Matlab基本操作和常用命令。
3.了解Matlab常用函数,运算符和表达式。
4.掌握Matlab工作方式和M文件的相关知识。
5.学会Matlab中矩阵和数组的运算。
二、【实验任务】P16 第4题编写函数文件,计算1!nkk =∑,并求出当k=20时表达式的值。
P27第2题矩阵A=123456789⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,B=468556322⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,计算A*B,A.*B,并比较两者的区别。
P27第3题已知矩阵A=5291⎡⎤⎢⎥⎣⎦,B=1292⎡⎤⎢⎥⎣⎦,做简单的关系运算A>B,A==B,A<B,并做逻辑运算(A==B)&(A<B),(A==B)&(A>B)。
P34 第1题用11114357π=-+-+……公式求π的近似值,直到某一项的绝对值小于-610为止。
三、【实验程序】P16 第4题function sum=jiecheng(n) sum=0;y=1;for k=1:nfor i=1:ky=y*i;endsum=sum+y;endsumP27第2题>>A=[1 2 3;4 5 6;7 8 9]>>B=[4 6 8;5 5 6;3 2 2]>>A*BP27第3题>> A=[5 2;9 1];B=[1 2;9 2];>>A>B>>A==B>>A<B>> (A==B)&(A<B)>> (A==B)&(A>B)P34 第1题t=1;pi=0;n=1;s=1;while abs(t)>=1e-6pi=pi+t;n=n+2;s=-s;t=s/n;endpi=4*pi;四、【实验结果】P16 第4题P27第2题两者的区别:A*B是按正规算法进行矩阵的计算, A.*B是对应元素相乘。
[IT计算机]MATLAB实验报告
![[IT计算机]MATLAB实验报告](https://img.taocdn.com/s3/m/1286c93d4a35eefdc8d376eeaeaad1f34693113e.png)
MATLAB实验报告姓名:学号:班级:专业:指导老师:实验一MATLAB运算基础一、实验目的:1、熟悉启动和退出MATLAB的方法。
2、熟悉MATLAB命令窗口的组成。
3、掌握建立矩阵的方法、4、掌握MATLAB各种表达式的书写规则以及常用函数的使用。
二、实验内容:1.1、>> t=0:0.5:2.5t =0 0.5000 1.0000 1.5000 2.0000 2.5000>>y=t.^2.*((t>=0)&(t<1))+(t.^2-1).*((t>=1)&(t>=1)&(t<2))+(t.^2-2*t+1).*((t>=2)&(t <3))y =0 0.2500 0 1.2500 1.0000 2.2500分析:该实验用逻辑表达式求出了分段函数的值。
1.2、>> p=rem([100:999],21)==0;>> sum(p)ans =43>> p1=rem([100:999],21),A=find(p1==0),length(A)ans =43分析:该实验求出了[100,999]之间能被21整除的数的个数为43.1.3、>> ch='ABcdefGHd',k=find(ch>='A'&ch<='Z'),ch(k)=[]ch =ABcdefGHdk =1 2 7 8ch =cdefd分析:该实验建立了一个字符串向量,删除了其中的大写字母。
1.4、A=[1 2 3;4 5 6;7 8 9],[i,j]=find(A>=5),for n=1:length(i)m(n)=A(i(n),j(n))endm>> mm =7 5 8 6 9分析:该实验求出了矩阵A中大于或等于5的元素为7 5 8 6 9.1.5、>>a11=input('all='),a12=input('a12='),a21=input('a21='),a22=input('a22='),A=[a 11,a12;a21,a22],DA=det(A),IN=inv(A),EA=eig(A)a11=23a12=12a21=11a22=32A =23 1211 32DA =604IN =0.0530 -0.0199-0.0182 0.0381EA =15.161039.8390分析:该实验求出了矩阵A中的行列式的值、逆和特征值。
matlab实验报告
![matlab实验报告](https://img.taocdn.com/s3/m/fd1137c489eb172ded63b73a.png)
《MATLAB及其工程应用》实验报告专业:班级:学号:姓名:指导老师:实验一实验名称:Practice how to initialize/save/load the MATLAB实验时间:第六周周三7-8节 实验目的:1、熟练的掌握MATLAB 变量的使用规则,学会创建变量。
2、熟练的掌握数组的相关知识并能加以运用。
实验要求:1、熟练命名合法的变量名。
2、能够应对关于数组的相关问题,如:数组的大小,子数组等。
实验步骤:1.4 假设u=1和v=3,用MATLAB 语句编写下列语句 a. v u 34 b.22)(2v u v +- c.333u v v - d.234v π 解:打开matlab7.0程序,新建一个M-file ,M-file 中输入:u=1;v=3;a=4*u/(3*v);b=2*v^(-2)/(u+v)^2;c=v^3/(v^3-u^3);d=4*pi*v^2/3;abcd然后点击运行,将得到以下结果:a =0.4444b =0.0139c =1.0385d =37.6991实验二实验名称:Practice how to use the MATLAB operators.实验时间:第七周周三7-8节实验目的:1、熟练的掌握数组和矩阵的使用规则。
2、熟练的掌握数组和矩阵的相关知识并能加以运用。
实验要求:1、熟练命名合法的变量名。
2、能够应对关于数组和矩阵的相关问题。
实验步骤:2.1看下面的数组回答有关问题⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=0.01.10.01.54.13.14.03.01.01.24.38.26.61.10.00.65.31.20.01.11array a. array1的大小是多少?b. array1(4,1)的值是多少?c. array1(:,1:2)的大小和值为多少?d. array1([1 3],end)的大小和值为多少?解:打开matlab7.0程序,新建一个M-file ,M-file 中输入:array1=[1.1,0.0,2.1,-3.5,6.0;0.0,1.1,-6.6,2.8,3.4;2.1,0.1,0.3,-0.4,1.3;-1.4,5.1,0.0,1.1,0.0];size(array1) %2.11aarray1(4,1) %2.1barray1(:,1:2) %2.1csize(array1(:,1:2)) %2.1carray1([1 3],end) %2.1dsize(array1([1 3],end)) %2.1d然后点击运行,将得到以下结果:ans =4 5ans =-1.4000ans =1.1000 00 1.10002.1000 0.1000-1.4000 5.1000ans =4 2ans =6.00001.3000ans =2 12.3 写出下面的数组的大小和内容.注意后面的数组可能根据前面数组的定义.a. a=1:2:5;b. b=[a' a' a'];c. c=b(1:2:3,1:2:3);d. d=a+b(2,:);e. w=[zeros(1,3) ones(3,1)' 3:5'];f. b([1 3],2)=b([3 1],2);解:>> a=1:2:5a =1 3 5>> size(a)ans =1 3>> b=[a' a' a']b =1 1 13 3 35 5 5>> size(b)ans =3 3>> c=b(1:2:3,1:2:3)c =1 15 5>> size(c)ans =2 2>> d=a+b(2,:)d =4 6 8>> size(d)ans =1 3>> w=[zeros(1,3) ones(3,1)' 3:5'w =0 0 0 1 1 1 3 4 5>> size(w)ans =1 9>> b([1 3],2)=b([3,1],2)b =1 5 13 3 35 1 5>> size(b)ans =3 32.4数组定义如下,写下面的子数组的内容⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=0.01.10.01.54.13.14.03.01.01.24.38.26.61.10.00.65.31.20.01.11array a. array1(3,:)b. array1(:,3)c. array1(1:2:3,[3 3 4])d. array1([1 1],:)解:>> array1=[1.1 0.0 2.1 -3.5 6.0;0.0 1.1 -6.6 2.8 3.4;2.1 0.1 0.3 -0.4 1.3 ;-1.45.1 0.0 1.1 0.0]array1 =1.1 02.1 -3.5 60 1.1 -6.6 2.8 3.42.1 0.1 0.3 -0.4 1.3-1.4 5.1 0 1.1 0>> array1(3,:)ans =2.1 0.1 0.3 -0.4 1.3>> array1(:,3)ans =2.1-6.60.3>> array1(1:2:3,[3 3 4])ans =2.1 2.1 -3.50.3 0.3 -0.4>> array1([1 1],:)ans =1.1 02.1 -3.5 61.1 02.1 -3.5 6实验三实验名称:Practice MATLAB display and plotting functions.实验时间:第九周周三7-8节实验目的:1、掌握三种常见函数:disp ,fprintf ,plot2、掌握画图函数的使用方法实验要求:1、熟练掌握plot 函数的使用方法2、能运用画图函数画一些较简单的图形实验步骤:2.10 球的位置和速度.如果一静止小球在离地h 0的地方以初速度v 0做垂直运动,其等等式为o o h t v gt t h ++=221)( (2.10) o v gt t v +=)( (2.11)其中g 为重力加速度(9.81m/(s*s)),h (t )代表在t 时刻小球的高度.v (t )代表在时刻t 小球的速度.编写一个MATLAB 程序,计算出每一秒钟的速度和高度,并打印出h ,v 关于时间t 的函数.确保在你的图中有合适的标签.解:v0=input('v0=');h0=input('h0=');t=-10:0.01:10;h=0.5.*(-9.81).*t.^2+v0.*t+h0;v=(-9.81).*t+v0;plot(t,h,'r--',t,v,'g');axis([-10,10,-50,50]);xlabel('t');ylabel('h and v');legend('h','v')结果:v0=10h0=102.13双曲余弦.双曲余弦的定义如下2cosh xx e e x -+=编写一个程序,计算出用户指定的x 的值对应的双曲余弦值.用这个程序计算3.0的双曲余弦值.和MATLAB 中的内建函数cosh(x)得到的值是否完全相同.用MATLAB 打印出这个函数的图象.当x 为何值时,这个函数有最小值?最小值为多少?解:注:以下结果均为x=3.0时的值。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
《数学实验》报告
实验名称数学实验
学院自动化
专业班级
2015年4月
一、 【实验目的】
使用MATLAB 求解常微分方程,包括常微分方程的符号解法以及数值解法 二、 【实验任务】 1、求解微分方程y
cos xsinx
y =
'。
2、用数值方法求解下列微分方程,用不同颜色和线形将y 和y '画在同一个图形
窗口里:
2t -1y -y t y ='+''
初始时间:0t 0=;终止时间:π=f t ;初始条件:1.0|y 0t == 2.0|y 0t ='=。
三、 【实验程序】
1、
y = dsolve('Dy = x * sin(x) / cos(y)', 'x') 2、
函数exf.m :
function xdot = exf(t,x)
xdot = [0 1;1 -t]*x + [0;1]*(1 - 2*t);
主程序:
t0 = 0; tf = pi; x0t = [0.1;0.2]; [t, x] = ode23('exf', [t0, tf], x0t); y = x(:,1); y2 = x(:,2);
plot(t, y, 'r-', t, y2, 'b--') legend('数值积分解y', '数值积分解y ’')
四、 【实验结果】 1、
y =
asin(C3 + sin(x) - x*cos(x))
2、
五、【实验总结】
1、对于高阶常微分方程,需要先将它转化为一阶常微分方程组,即状态方程,所以求解常微分方程时,需要先用学过的知识对方程进行降阶处理;
2、运用数值解法求解微分方程时需要建立模型的函数文件,然后运用MATLAB中函数调用的知识在主程序中调用子程序,最终解出微分方程的解;
3、微分方程有数值积分解也有解析解,两者有各自优缺点,不同的求解方法适用于不同的场合,具体方程可以将两者互相对比,画出对比的图形;
4、
→解微分方程→用图形显示。