实验七在MATLAB中输入输出数据文件 (2)

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全部实验及答案

MATLAB全部实验及答案实验一、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*B就是线代里面的矩阵相乘 A.*B是对应位置的元素相乘（3）设a=10，b=20；求i=a/b=0.5与j=a\b=2？（4）设a=[1 -2 3;4 5 -4;5 -6 7]请设计出程序，分别找出小于0的矩阵元素及其位置（单下标、全下标的形式），并将其单下标转换成全下标。

clear,clca=[1 -2 3;4 5 -4;5 -6 7];[x,y]=find(a<0);c=[];for i=1:length(x)c(i,1)=a(x(i),y(i));c(i,2)=x(i);c(i,3)=y(i);c(i,4)=(y(i)-1)*size(a,2)+x(i);endc（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=？a^2= 22 16 1625 26 2326 24 28a.^2=1 4 99 16 425 4 9（7）有一段指令如下，请思考并说明运行结果及其原因clearX=[1 2;8 9;3 6];X( : ) 转化为列向量（8）使用三元组方法，创建下列稀疏矩阵2 0 8 00 0 0 10 4 0 06 0 0 0方法一：clear,clcdata=[2 8 1 4 6];ir=[1 1 2 3 4 ];jc=[1 3 4 2 1];s=sparse(ir,jc,data,4,4);full(s)方法二：不用三元组法clear,clca=zeros(4,4);a(1,[1,3])=[2,8];a(2,4)=1;a(3,2)=4;a(4,1)=6;a（9）写出下列指令的运行结果>> A = [ 1 2 3 ]; B = [ 4 5 6 ];>> C = 3.^A>> D = A.^B3、已知⎪⎭⎫ ⎝⎛+⋅=-334sin 234πt e y t 若需要计算t ∈[-1,1]，取间隔为0.01，试计算出相对应的y 值。

Matlab的使用方法及步骤详解

Matlab的使用方法及步骤详解一、Matlab简介Matlab是一种非常流行的科学计算软件，其全称为Matrix Laboratory（矩阵实验室）。

Matlab具有强大的数学计算和数据分析能力，广泛应用于工程、科学、经济等领域。

本文将详细介绍Matlab的使用方法及步骤。

二、安装与启动Matlab1. 下载与安装首先，访问MathWorks官方网站，找到适用于您操作系统的Matlab版本，并下载安装程序。

安装程序将引导您进行安装，按照提示完成即可。

2. 启动Matlab安装完成后，您可以在开始菜单或桌面上找到Matlab的启动图标。

点击启动图标，Matlab将打开并显示初始界面。

三、Matlab基本操作1. 工作区与编辑器Matlab的界面主要由工作区和编辑器组成。

工作区显示变量及其值，可用于查看和操作数据。

编辑器则用于编写和编辑Matlab脚本、函数等。

2. 脚本与命令窗口Matlab提供了两种主要的运行方式：脚本和命令窗口。

脚本是一系列命令的集合，可以一次性执行，适用于较复杂的计算任务。

命令窗口则可逐行输入命令并立即执行，用于快速测试和调试。

3. 基本算术和数学运算Matlab支持各种基本算术和数学运算，如加减乘除、幂运算、三角函数等。

可以直接在命令窗口输入表达式并执行。

四、数据操作与处理1. 数组的创建与操作在Matlab中，数组是最基本的数据结构之一。

可以使用多种方法创建数组，例如手动输入、加载外部文件、使用特定函数等。

一旦创建，可以对数组进行各种操作，如索引、切片、拼接等。

2. 矩阵运算Matlab对矩阵运算提供了强大的支持。

可以进行矩阵加减乘除、转置、求逆等运算。

矩阵运算在解决线性方程组、最小二乘拟合等问题时非常有用。

3. 数据可视化Matlab提供了丰富而强大的数据可视化功能。

使用plot、scatter、histogram等函数可以绘制各种类型的图表。

还可以对图表进行格式设置、添加标签、调整坐标轴等。

matlab数学实验

《管理数学实验》实验报告班级姓名实验1：MATLAB的数值运算【实验目的】（1）掌握MATLAB变量的使用（2）掌握MATLAB数组的创建，（3）掌握MA TLAB数组和矩阵的运算。

（4）熟悉MATLAB多项式的运用【实验原理】矩阵运算和数组运算在MA TLAB中属于两种不同类型的运算，数组的运算是从数组元素出发，针对每个元素进行运算，矩阵的运算是从矩阵的整体出发，依照线性代数的运算规则进行。

【实验步骤】（1）使用冒号生成法和定数线性采样法生成一维数组。

（2）使用MA TLAB提供的库函数reshape，将一维数组转换为二维和三维数组。

（3）使用逐个元素输入法生成给定变量，并对变量进行指定的算术运算、关系运算、逻辑运算。

（4）使用MA TLAB绘制指定函数的曲线图，将所有输入的指令保存为M文件。

【实验内容】（1）在[0,2*pi]上产生50个等距采样数据的一维数组，用两种不同的指令实现。

0:(2*pi-0)/(50-1):2*pi 或linspace(0,2*pi,50)（2）将一维数组A=1:18，转换为2×9数组和2×3×3数组。

reshape(A,2,9)ans =Columns 1 through 71 3 5 7 9 11 132 4 6 8 10 12 14Columns 8 through 915 1716 18reshape(A,2,3,3)ans(:,:,1) =1 3 52 4 6ans(:,:,2) =7 9 118 10 12 ans(:,:,3) =13 15 17 14 16 18（3）A=[0 2 3 4 ;1 3 5 0],B=[1 0 5 3;1 5 0 5]，计算数组A 、B 乘积，计算A&B,A|B,~A,A= =B,A>B 。

A.*Bans=0 0 15 121 15 0 0 A&Bans =0 0 1 11 1 0 0 A|Bans =1 1 1 11 1 1 1~Aans =1 0 0 00 0 0 1A==Bans =0 0 0 01 0 0 0A>=Bans =0 1 0 11 0 1 0（4）绘制y= 0.53t e -t*t*sin(t),t=[0,pi]并标注峰值和峰值时间,添加标题y= 0.53t e -t*t*sint ，将所有输入的指令保存为M 文件。

实验7串行接口输入输出实验

北京林业大学11学年—12学年第 2 学期计算机组成原理实验任务书专业名称：计算机科学与技术实验学时： 2课程名称：计算机组成原理任课教师：张海燕实验题目：实验七串行接口输入输出实验实验环境：TEC-XP+教学实验系统、PC机实验内容1．串行接口输入输出；2．串行接口扩展。

实验目的学习串行口的正确设置与使用。

实验要求1．实验之前认真预习，明确实验的目的和具体实验内容，做好实验之前的必要准备。

2．想好实验的操作步骤，明确通过实验到底可以学习哪些知识，想一想怎么样有意识地提高教学实验的真正效果；3．在教学实验过程中，要爱护教学实验设备，记录实验步骤中的数据和运算结果，仔细分析遇到的现象与问题，找出解决问题的办法，有意识地提高自己创新思维能力。

4．实验之后认真写出实验报告，重点在于预习时准备的内容，实验数据，运算结果的分析讨论，实验过程、遇到的现象和解决问题的办法，自己的收获体会，对改进教学实验安排的建议等。

善于总结和发现问题，写好实验报告是培养实际工作能力非常重要的一个环节，应给以足够的重视。

必要知识串行接口是计算机主机和某些设备之间实现通信，硬件造价比较低廉、标准化程度比较高的一种输入输出接口线路，缺点是通信的速度比较低。

从在程序中使用串行接口芯片的角度看，接口芯片内有用户可以访问的4个寄存器，分别是接收CPU送来数据的输出数据缓冲寄存器，向CPU提供数据的输入数据缓冲寄存器，接收CPU发来的控制命令的控制寄存器，向CPU提供接口运行状态的状态寄存器，必须有办法区分这4个寄存器。

接口芯片中还有执行数据串行和并行转换的电路，接口识别电路等。

串行接口用于执行数据的输入输出操作。

一次输入或输出操作通常需要两个操作步骤完成，第一步是为接口芯片提供入出端口地址，即把指令寄存器低位字节的内容（8位的IO端口地址）经过内部总线和运算器部件写进地址寄存器AR，第二步是执行输入或输出操作，若执行输入指令IN，则应从接口芯片读出一个8位的数据并经过数据总线DB和内部总线IB写进寄存器堆中的R0寄存器，若执行OUT指令，则需要把R0寄存器的内容经过内部总线IB和数据总线DB写入接口芯片。

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操作环境。

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=（x-1）/(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三、设计提示1．初次接触Matlab应该注意函数表达式的文本式描述。

实验七变量、伪指令的使用及字符串的输入输出

（1）输入程序并检查无误，经汇编、连接后产生正确的可执行文件 SY7-1.EXE。（2）执行可执行文件，可以看到执行正常，但是没有任何输出，可以使用 DEBUG来调试。（3）用 DEBUG 的 U CS:0000 查看 MOV AX,XXXX(DATA)语句，即得到数据段段地址 DS=XXXX。（4）用 G 命令执行，调试过程如下： C:\ ＞ DEBUG SY7-1.EXE 一 U 12AF:0000 B8AE12 12AF:0003 8ED8 12AF:0005 BE0000 12AF:0008 8A04 12AF:000A B90900 12AF:000D 46 12AF:000E 3A04 12AF:0010 7302 12AF:0012 8A04 INC CMP JNB MOV MOV MOV MOV DS,AX MOV AL,[SI] MOV SI AL,[SI] 0014 AL,[SI] CX,0009 SI,0000 AX,12AE
实验注意事项 1．指令语句以 80X86 指令系统的助记符为基础构成,经汇编后将产生相应的机器代码而构成目标程序，供机器执行。 2．伪指令语句是为汇编程序和连接程序提供一些必要控制信息的、由伪指令构成的管理性语句，其对应的伪操作是在汇编过程中完成，汇编后不产生机器代码。实验与思考 1．已知字符串‘ABCDEFGH’位于数据段开始的区域中，编程序将该字符串转换为小写后，存于其后。 2．编写程序完成如下功能：输入一个字符串，在下一行原样输出。
Hale Waihona Puke 12AF:0014 E2F7 LOOP 000D 12AF:0016 A20A00 MOV [000A],AL 12AF:0019 B44C MOV AH,4C 12AF:001B CD21 INT 21 12AF:001D C45EF4 LES BX,[BP-0C] 一 G=0 19 AX=1290 BX=0000 CX=0000 DX=0000 SP=0000 BP=0000 SI=0009 DI=0000 DS=12AE ES=129E SS=12AE CS=12AF IP=0019 OV UP EI PL NZ NA PO NC 12AF:0019 B44C MOV AH,4C 一 D 0000 000f;最大数为 90H 12AE:0000 60 50 20 40 80 90 00 10-30 70 90 00 00 00 00 00 2. 数据段开始区域中，连续存放着10_个无符号数，编程序找出这 _10_个数中奇数的个数，并将其存到该数据区的后面。参考程序如下：（程序名：SY7-2.ASM） data segment buf1 db 64h,52h,25h,47h,86h,98h,30h,17h,34h,79h buf2 db ? data ends code segment assume cs:code,ds:data start:

MATLAB实验指导

MATLAB基础教程实验指导书实验一：Desktop操作桌面基础一、实验目的及要求1、熟悉MATLAB系统的安装流程，掌握MATLAB的启动和退出。

2、掌握MATLAB系统的各命令窗口的功能，熟悉常用选项和工具栏的功能和用途。

3、熟悉简单程序的输入、运行、调试及结果的显示过程。

二、实验内容1、认识MATLAB集成环境：熟悉个操作窗口的功能和用途。

掌握File(文件)、Edit(编辑)、View(显示)、Web(网络)、Window(窗口)和Help(帮助) 等菜单命令的使用。

2、启动和退出MATLAB（1）启动MATLAB的M文件。

在启动MATLAB时，系统可自动执行主M文件matlabrc.m，在matlabrc.m的末尾还会检测是否存在startup.m，如存在则会自动执行它。

在网络系统中，matlabrc.m保留给系统管理员，而各个用户可利用startup.m进行初始设置。

(2)、终止或退出MATLAB。

quit命令可终止MATLAB，但不保存工作空间的内容。

为保存工作空间的内容，可使用save命令。

1、利用save、load命令，保存和恢复工作空间。

用clear命令可清空工作空间。

（1）、工作空间中的变量可以用save命令存储到磁盘文件中。

（2）、用load命令可将变量从磁盘文件读入MATLAB的工作空间。

（3）、用clear命令可清除工作空间中现存的变量。

4、MATLAB的所有图形工具窗体都可以嵌入MATLAB窗体(Dock)，也可以从MATLAB窗体中弹出(Undock)，例如在MATLAB默认的图形窗体环境下，单击命令行窗体左上角按钮，就可以将MATLAB命令行窗体弹出。

要求分别将命令行窗体(Command Window)、命令行历史窗体(Command History)、当前路径查看器(Current Directory)、工作空间浏览器(Workspace Browser)、帮助(Help)、MATLAB性能剖析工具(Profiler) 从MATLAB窗体中弹出和嵌入MATLAB窗体(Dock)。

matlab实验内容答案解析

实验报告说明:matlab课程实验需撰写8个实验报告,每个实验报告内容写每次实验内容中标号呈黑体大号字显示的题目。

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

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

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

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

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

(1) Zl =2 sin85。

1 + e22-0.45,0.3“ _ -0.3d提示：利用冒号表达式生成a向量，求各点的函数值时用点乘运算。

t20</<1(4)Z4=$—11</<2 ，其中t = 0 : 0.5 : 2.52<t<3提示:用逻辑表达式求分段函数值。

2 •已知_12 34 -4_i 3 -rA =34 7 87,B = 2 0 33 65 73-2 7求下列表达式的值：(1) A+6二B和A-B+I(其中I为单位矩阵)。

(2) A*B和A.*B o(3) A^3 和A\3。

(4) A/B和B\A。

(5 ) [A , B]和[A([l f3],；)；B A2]。

3•设有矩阵A和B_12345・~301667891017-69A =1112131415,B =023-41617181920970212223242541311(1) 求它们的乘积C。

2 将矩阵C的右下角3x2子矩阵赋给D(3) 查看MATLAB I作空间使用情况。

4.完成下列操作：(1) 求［100,999］之间能被21整除的数的个数。

提示:先利用冒号表达式,再利用find和length函数。

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

提示：利用find函数和空矩阵。

第二次实验内容实验三选择结构程序设计一、实验目的1. 掌握建立和执行M文件的方法。

2. 掌握利用if语句实现选择结构的方法。

Matlab实验报告（1）

《Matlab语言与应用》课程实验报告*名：**班级：电信114学号：************指导老师：***二〇一三年十一月二十一日Matlab实验报告实验一一、实验问题求[12 + sin(2)×( 22 −4)]÷3^2的算术运算结果。

二、问题的分析该题主要熟悉Matlab环境下的对数学运算的熟悉，如何输入数据、建立函数输出结果。

三、上机实验结果如图四、实验的总结与体会通过本次实验，我学会了用Matlab来计算数学运算中的复杂技术。

我们也可以自己编写一个可以调用的函数，首先我们要了解Matlab语言函数的基本结构，结构如下：Function [返回变量列表]=函数名（输入变量列表）注释说明语句段，由%引导输入、返回变量格式的检测函数体语句例如：输入变量为k，返回的变量为m和s,其中s为前m项的和Function [m,s]=findsum(k)s=0;m=0;while(s<=k),m=m+1;s=s+m;end编写了函数，就可以将其存为findsum.m文件，这样就可以在Matlab环境中对不同的k值调用该函数了。

这样就可以灵活的实现我们想要的数据。

实验二一、实验问题二、问题分析输入矩阵时，空格或逗号表示间隔，分号表示换行，比如上面的矩阵A应写为A=[1,2,3;4,5,6;7,8,9]三、上机实验结果如图四、实验的总结与体会通过对本次上机实验了解到在Matlab中对矩阵的代数运算矩阵转置、矩阵的加减法运算、矩阵乘法、矩阵的左、右除、矩阵翻转、矩阵乘方运算、点运算等。

实验三一、实验要求画图，理解plot函数用法二、代码如下clear; clf;t=0:pi/20:2*pi;R=5;x=R*sin(t); y=R*cos(t);plot(x,y,'b:'), gridhold onrrr=[x;y;x+y];plot(rrr(1,:),'.','MarkerSize',10,'Color','r')plot(rrr(2,:),rrr(3,:),'o','MarkerSize',15,'Color','b'); axis([-8,20,-8,8]), % axis squarehold off三、Matlab运行结果如图实验四一、实验要求二、代码如下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运行结果如图实验五一、实验要求傅里叶频谱分析二、代码及分析(1)首先生成数据，包含50Hz和120Hz频率的正弦波x >>t = 0:.001:.25;>>x = sin(2*pi*50*t) + sin(2*pi*120*t);(2)再生成噪音信号yy = x + 2*randn(size(t));plot(y(1:50))title('Noisy time domain signal')(3)对y进行快速傅里叶变换Y = fft(y,256);(4)计算功率谱Pyy = Y.*conj(Y)/256;f = 1000/256*(0:127);plot(f,Pyy(1:128))title('Power spectral density')xlabel('Frequency (Hz)')(5)只查看200Hz以下频率段plot(f(1:50),Pyy(1:50))title('Power spectral density')xlabel('Frequency (Hz)')三、上机结果如下图实验六一、实验要求FIR数字滤波器设计代码如下clear;close allf=[0,0.6,0.6,1]; m=[0,0,1,1]; % 设定预期幅频响应b=fir2(30,f,m); n=0:30; % 设计FIR 数字滤波器系数subplot(3,2,1),stem(n,b,'.')xlabel('n'); ylabel('h(n)');axis([0,30,-0.4,0.5]),line([0,30],[0,0])[h,w]=freqz(b,1,256);subplot(3,2,2),plot(w/pi,20*log10(abs(h)));gridaxis([0,1,-80,0]),xlabel('w/pi'); ylabel('幅度(dB)');二、上级结果如图实验七二、实验要求用guide实验一个简单的加减乘除计算器二、实验步骤在命令行输入guide命令，进入guide界面新建一个空白guide文件在空白文件中设置好功能模块如图模块建立好后，就要把编写好的加减乘除代码加入到各自的回调函数中，见下图两个被加数代码如下图加模块代码如下图减模块代码如下图乘模块代码如下图除模块代码如下图各模块的回调函数加完后就可以运行了，运行结果如下图总结：Matlab一个高级的距阵/阵列语言，它包含控制语句、函数、数据结构、输入和输出和面向对象编程特点。

matlab数学实验报告

MATLAB数学实验报告指导老师：班级：小组成员：时间：201_/_/_Matlab第二次实验报告小组成员：1题目：实验四;MATLAB选择结构与应用实验目的：掌握if选择结构与程序流程控制;重点掌握break;return;pause语句的应用..问题：问题1：验证“哥德巴赫猜想”;即：任何一个正偶数n>=6均可表示为两个质数的和..要求编制一个函数程序;输入一个正偶数;返回两个质数的和..问题分析：由用户输入一个大于6的偶数;由input语句实现..由if判断语句判断是否输入的数据符合条件..再引用质数判断函数来找出两个质数;再向屏幕输出两个质数即可..编程：function z1;z2=geden;n=input'please input n'if n<6disp'data error';returnendif modn;2==0for i=2:n/2k=0;for j=2:sqrtiif modi;j==0k=k+1;endendfor j=2:sqrtn-iif modn-i;j==0k=k+1;endendif k==0fprintf'two numbers are'fprintf'%.0f;%.0f';i;n-ibreakendendend结果分析：如上图;用户输入了大于6的偶数返回两个质数5和31;通过不断试验;即可验证哥德巴赫猜想..纪录：if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰;更快的解决问题..2题目：实验四;MATLAB选择结构与应用实验目的：用matlab联系生活实际;解决一些生活中常见的实际问题..问题：问题四：在一边长为1的四个顶点上各站有一个人;他们同时开始以等速顺时针沿跑道追逐下一人;在追击过程中;每个人时刻对准目标;试模拟追击路线;并讨论.. （1）四个人能否追到一起（2）若能追到一起;每个人跑过多少路程（3）追到一起所需要的时间设速率为1问题分析：由正方形的几何对称性和四个人运动的对称性可知;只需研究2个人的运动即可解决此问题..编程：hold onaxis0 1 0 1;a=0;0;b=0;1;k=0;dt=0.001;v=1;while k<10000d=norma-b;k=k+1;plota1;a2;'r.';'markersize';15;plotb1;b2;'b.';'markersize';15;fprintf'k=%.0f b%.3f;%.3f a%.3f;%.3f d=%.3f\n';k;b1;b2;a1;a2;da=a+b1-a1/d*dt;b2-a2/d*dt;b=b+b2-a2/d*dt;-b1-a1/d*dt;if d<=0.001breakendendfprintf'每个人所走的路程为：%.3f';k*v*dtfprintf'追到一起所需要的时间为%.3f';k*dt结果分析：上图为2人的模拟运动路线;有对称性可解决所提问题..-上图为运算过程和运算结果..四个人可以追到一起;走过的路程为1.003;时间也为1.003.纪录：此题利用正方形和运动的对称性可以简便运算..3题目：实验八;河流流量估计与数据插值目的：由一些测量数据经过计算处理;解决一些生活实际问题..问题：实验八上机练习题第三题：瑞士地图如图所示;为了算出他的国土面积;做以下测量;由西向东为x轴;由南向北为y轴;从西边界点到东边界点划分为若干区域;测出每个分点的南北边界点y1和y2;得到以下数据mm..已知比例尺1:2222;计算瑞士国土面积;精确值为41288平方公里..测量数据如下：x=7.0 10.5 13.0 17.5 34 40.5 44.5 48 56 61 68.5 76.5 80.5 91 96 101 104 106 111.5 118 123.5 136.5 142 146 150 157 158 ;y1=44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68;y2=44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68;问题分析：先由题目给定的数据作出瑞士地图的草图;再根据梯形法;使用trapz语句;来估算瑞士国土的面积..编程：x=7.0 10.5 13.0 17.5 34 40.5 44.5 48 56 61 68.5 76.5 80.5 91 96 101 104 106 111.5 118 123.5 136.5 142 146 150 157 158;y1=44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68;y2=44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68;plotx;y1;'r.';'markersize';15;plotx;y2;'r.';'markersize';15;axis0 160 0 135grid;hold ont=7:158;u1=splinex;y1;t;u2=splinex;y2;t;plott;u1plott;u2s1=trapzt;u1;s2=trapzt;u2;s=s2-s1*2222*22222/10000000;fprintf'S=%.0f';s结果分析：上图为由所给数据绘制出的瑞士地图..上图为运算结果;计算出瑞士的国土面积为42472平方公里;与准确值41288较为接近..纪录：使用梯形分割的方法;trapz语句可以方便计算不规则图形面积;但存在一定误差..4题目：实验七：圆周率的计算与数值积分目的：将数值积分最基本的原理应用于matlab之中;解决一些与积分有关的问题..问题：实验七上机练习题第一题：排洪量某河床的横断面如图7.3所示;为了计算最大排洪量;需要计算其断面积;试根据所给数据m用梯形法计算其断面积..问题分析：河床断面可近似分割成若干曲边梯形;近似处理把它们当做梯形来计算面积可使问题得到简化..编程：clc;clear;x=0 4 10 12 15 22 28 34 40;y=0 1 3 6 8 9 5 3 0;y1=10-y;plotx;y1;'k.';'markersize';15;axis0 40 0 10;grid;hold ont=0:40;u=splinex;y1;t;plott;u;s=40*10-trapzt;u;fprintf's=%.2f\n';s结果分析：上图为河床的断面图..上图为计算结果面积约为180.70平方米..纪录：使用梯形法计算不规则图形面积十分简便易行..5题目：实验七：圆周率的计算与数值积分目的：使用matlab计算解决一些有关积分的问题..问题：实验七上机练习题第三题：从地面发射一枚火箭;在最初100秒内记录其加速度如下;试求火箭在100秒时的速度..Ts=0 10 20 30 40 50 60 70 80 90 100;Am/s*s=30.00 31.63 33.44 35.47 37.75 40.33 43.29 46.69 50.67 54.01 57.23;问题分析：加速度为速度的微分;已知微分求积分;类似于面积问题;可使用梯形法来计算..编程：clc;clear;x=0 10 20 30 40 50 60 70 80 90 100;y=30.00 31.63 33.44 35.47 37.75 40.33 43.29 46.69 50.6754.01 57.23;plotx;y;'k.';'markersize';15;axis0 100 20 60;grid;hold ons=0:10:100;z=splinex;y;s;plots;y;v=trapzx;y;fprintf'v=%.2f\n';v结果分析：上图为加速度变化图..上图为计算结果;求得火箭在100秒时速度约为4168.95m/s..纪录：梯形法可以推广解决许多已知微分求积分的其他问题..6题目：实验七：圆周率的计算与数值积分目的:计算曲线弧长闭曲线周长可使用微元法;ds=sqrtdx^2+dy^2;在转化微积分问题;累加即可得到结果..问题：实验七上机练习题第三题：计算椭圆想x^2/4+y^2=1的周长;使结果具有五位有效数字..问题分析：编程：s=0;dx=0.001;for x=0:0.001:1.999dy=1.-x+0.001.^2/4-1.-x.^2/4;ds=sqrtdx.^2+dy.^2;s=s+ds;ends=4*s;fprintf'the length is'fprintf'%.4f';s结果分析：上图为计算结果;给定椭圆的周长约为9.1823五位有效数字纪录：计算不规则曲线弧长;可使用微元法;划分为若干小的看做直角三角形;利用勾股定理解决..7题目：实验九人口预测与数据拟合目的：掌握一些曲线拟合的方法;了解曲线拟合常用函数..问题：用电压U=10v的电池给电容器充电;t时刻的电压Vt=U-U-V0exp-t/τ;其中V0是电容器的初始电压;τ是充电常数;由所给数据确定V0和τ..t=0.5 1 2 3 4 5 7 9;V=3.64 3.52 2.74 1.78 1.34 1.01 0.57 0.37;问题分析：题中已给出函数关系式;为指数函数曲线拟合;将所给函数式整理可得标准的exp形函数曲线;从而便于解决..编程：t=0.5 1 2 3 4 5 7 9;V=3.64 3.52 2.74 1.78 1.34 1.01 0.57 0.37;plott;V;'k.';'markersize';20;axis0 10 0 4;grid;hold onpause0.5n=8;a=sumt1:n;b=sumt1:n.*t1:n;c=sumlogV1:n;d=sumt1:n.*logV1:n;A=n a;a b;B=c;d;p=invA*Bx=0:10;y=expp1+p2*x;plotx;y;'r-';'linewidth';2结果分析：上图为电压与时间关系图..上图为计算结果;即U-V0=1.4766;所以V0=8.5234;-1/τ=-0.2835;所以τ=3.5273纪录：曲线拟合的一个重难点是选择合适的曲线函数;才能提高拟合度..8题目：实验七圆周率的计算与数值积分目的：拓展圆周率的各种计算方法;掌握其他数值的近似计算方法..问题：实验七练习2：计算ln2的近似值精确到10的-5次方（1）利用级数展开的方法来计算（2）利用梯形法计算（3）利用抛物线法问题分析：级数展开;梯形法;抛物线法是常见的近似运算方法..编程：1级数展开的方法clc;clear;n=0;r=1;p=0;k=-1;while r>=0.1e-5n=n+1;k=k*-1;p1=p+k/n;r=absp1-p;fprintf'n=%.0f;p=%.10f\n';n;p1;p=p1;end2梯形法clc;clear;f=inline'1./x';x=1:0.1:2;y=fx;p=trapzx;y;fprintf'p=%.6f\n';p3抛物线法clc;clear;f=inline'1./x';a=1;b=2;n=1;z=quadf;a;b;fprintf'z=%.10f\n';z结果分析：（1）级数展开的方法（2）梯形法3抛物线法纪录：级数展开法;梯形法;抛物线法;计算近似值时应合理利用..梯形法和抛物线法不易提高精确度;级数展开法可以提高精确度..9题目：实验八河流流量估计与数据插值目的：掌握求插值多项式的方法;并利用此计算近似值..问题：已知y=fx的函数表如下x=0.40 0.55 0.65 0.80 0.90 1.05;y=0.41075 0.57815 0.69675 0.88811 1.02652 1.25382;求四次拉格朗日插值多项式;并由此求f0.596问题分析：利用所给函数表可计算拉格朗日插值多项式..编程：function p=lagrangex;yL=lengthx;a=onesL;for j=2:La:;j=a:;j-1.*x';endx=inva*y';for i=1:Lpi=xL-i+1;endx=0.40 0.55 0.65 0.80 0.90 1.05;y=0.41075 0.57815 0.69675 0.88811 1.02652 1.25382; plotx;y;'k.';'markersize';15axis0 2 0 2grid;hold on;p=lagrangex;y;t=0:0.1:1.5;u=polyvalp;t;plott;u;'r-'a=polyvalp;0.596结果分析：上图为所求结果;估算值和插值多项式..纪录：插值多项式是一项十分实用的方法..10题目：求正整数n的阶乘：p=1*2*3*…*n=n;并求出n=20时的结果目的：练习使用循环变量解决数学问题问题：对程序：Clear;clc;n=20;p=1;for i=1:np=p*i;fprintf’i=%.0f;p=%.0f\n’;i;pend进行修改使它：利用input命令对n惊醒赋值问题分析：题中给出程序中“n=20”修改;使用input命令；讲题中的输出命令放出循环之外..编程：clear;clc;n=input'n=';p=1;for i=1:np=p*i;endfprintf'i=%.0f;p=%.0f\n';i;p结果：n=20i=20;p=2432902008176640000>>结果分析：使用input命令可以实现人机对话;使用户自由赋值；输出语句在程序中的位置对输出的结果有很大的影响；在循环内部可以在计算过城中不断输出结果；在循环之外则可以控制只输出最后结果..11题目：对于数列｛√2｝;n=1;2;…;求当其前n项和不超过1000时的n的值及合的大小..目的：运用条件循环解决文帝个项数的循环程序求解;问题：对程序：clear;clc;n=0;s=0;while s<=1000n=n+1;s=s+sqrtn;fprintf’n=%.0f;s=%.4f\n’;n;send问题分析：题中所给程序中的限制变量为上次循环之后的s;导致s超过上限后仍有一次的循环;若把循环变量改为这次的s;则可以避免这种情况的发生..编程：clear;clc;n=0;s=0;while s+sqrtn<=1000n=n+1;s=s+sqrtnfprintf'n=%.0f;s=%.4f\n';n;send结果：……s =970.8891n=128;s=970.8891s =982.2469n=129;s=982.2469s =993.6487n=130;s=993.6487>>结果分析：从结果中可以看出;最后一步为我们需要的答案;从这道题我们可以得出循环变量对一道编程的重要性..。

matlab中将数据输出保存为txt格式文件的方法 (1)

将matlab中数据输出保存为txt或dat格式总结网上各大论坛，主要有三种方法。

第一种方法：save（最简单基本的）具体的命令是：用save *.txt -ascii xx为变量*.txt为文件名,该文件存储于当前工作目录下，再打开就可以打开后,数据有可能是以指数形式保存的.例子：a =[17 24 1 8 15;23 5 7 14 16 ;4 6 13 20 22 ;10 12 19 21 3 ;11 18 25 2 9 ]；save afile.txt -ascii aafile.txt打开之后，是这样的：1.7000000e+0012.4000000e+001 1.0000000e+000 8.0000000e+000 1.5000000e+0012.3000000e+001 5.0000000e+000 7.0000000e+000 1.4000000e+001 1.6000000e+0014.0000000e+000 6.0000000e+000 1.3000000e+001 2.0000000e+001 2.2000000e+0011.0000000e+001 1.2000000e+001 1.9000000e+0012.1000000e+0013.0000000e+0001.1000000e+001 1.8000000e+0012.5000000e+001 2.0000000e+000 9.0000000e+000第二种方法：dlmwritedlmwrite('a.txt',a,'precision','%10.0f')或者是dlmwrite('a.txt',a,'delimiter', '\t')对于只有一行或者一列的数据，很适用，但是多行的，就乱了网上有很多这一类似的问题，但是都不是很理想第三种方法：fopen+fprintf下面主要介绍这种方法，可以解决以上问题：用fprintf命令：以上面的例子为例：第一种情况：>> a=[17 24 1 8 15;23 5 7 14 16 ;4 6 13 20 22 ;10 12 19 21 3 ;11 18 25 2 9 ];>> fid = fopen('b.txt','wt');fprintf(fid,'%g\n',a); # \n 换行fclose(fid);然后用写字板打开b.txt，内容如下：为列向量17234101124561218171319258142021215162239第二种情况：对上面的命令做一下改动：# \n 换行改为\t，table键>> fid = fopen('b.txt','w');fprintf(fid,'%g\t',a);fclose(fid);然后用写字板打开b.txt，内容如下：为行向量：17 23 4 10 11 24 5 6 12 18 1 7 13 19 25 8 14 20 21 2 15 16 22 3 9 第三种情况：综合上面的两个结果，我们编写以下命令：<A style="key: ,'wt');%写入文件路径[m,n]=size(a);for i=1:1:mfor j=1:1:nif j==nfprintf(fid,'%g\n',a(i,j));elsefprintf(fid,'%g\t',a(i,j));endendendfclose(fid);然后用写字板打开b.txt，内容如下：矩阵17 24 1 8 1523 5 7 14 164 6 13 20 2210 12 19 21 311 18 25 2 9说明：以上操作都是在当前的工作目录下完成！下面给出最一般的模型，大家可以试着自己操作,如果需要dat格式，直接把txt换为dat就可以fid=fopen('C:\Documents and Settings\cleantotal.ped','wt');%写入文件路径matrix=input_mattrix %input_matrix为待输出矩阵[m,n]=size(matrix);for i=1:1:mfor j=1:1:nif j==nfprintf(fid,'%g\n',matrix(i,j));elsefprintf(fid,'%g\t',matrix(i,j));endendendfclose(fid);×××××××××××××××××××××××××××××××××××××××××××××××××××××××××下面附了具体的matlab的fopen和fprintf函数具体解释，当然help一下是可以知道的，只是为了方便大家matlab中fopen函数在指定文件打开的实例如下：*1)“fopen”打开文件，赋予文件代号。

Matlab处理数据导出Paraview可读的vtk文件（二）

Matlab处理数据导出Paraview可读的vtk⽂件（⼆）由于我在⽤SPH⽅法仿真时⽤的是FORTRAN语⾔，并且没有找到直接输出vtk⽂件的代码，因此偷懒通过MATLAB转换⼀下数据。

⽤到的Matlab⼦程序可通过⼀下链接找到。

假设我每个粒⼦输出的格式按照下⾯输出的x(i), y(i), z(i), vx(i), vy(i), vz(i), p(i), rho(i), phase(i)含义为每个粒⼦的坐标，速度，压强，密度，相。

vtkwrite(‘output.vtk’,'unstructured_grid',x(:),y(:),z(:),'vectors','velocity',vx(:),vy(:),vz(:),'scalars', 'pressure',p(:),'scalars', 'rho',rho(:),'scalars', 'phase',phase(:));解释：1，输出⽂件名2，定义⾮结构⽹格，还有其他选项3，后⾯跟每个点的坐标4，‘vectors’定义向量，后⾯‘velocity’是⾃定义的名字，⽐如速度，然后后⾯跟速度的三个分量5，‘scalars’定义数值，后⾯引号中⾃定义名字，⽐如压强和密度等6，向量和数值最后在paraview中都有对应的显⽰⽅式。

将⽣成的vtk⽂件拖进paraview就可以在Representation中选择Point Gaussian。

在Coloring中选择要显⽰的数据⼆、显⽰向量1、选择菜单中filter->Alphabetical->Glyph或者点下⾯/2、设置你要显⽰的向量。

在Matlab中读取、写入二进制数据文件（转）

在Matlab中读取、写⼊⼆进制数据⽂件（转）Matlab中关于⼆进制数据⽂件操作所使⽤的语⾔是标准的C，有关数据⽂件的编程⽅法与LabVIEW类似。

1、⽂件的打开与关闭1）打开⽂件在读写⽂件之前，必须先⽤fopen函数打开或创建⽂件，并指定对该⽂件进⾏的操作⽅式。

fopen函数的调⽤格式为：fid=fopen（⽂件名，‘打开⽅式’）说明：其中fid⽤于存储⽂件句柄值，如果返回的句柄值⼤于0，则说明⽂件打开成功。

⽂件名⽤字符串形式，表⽰待打开的数据⽂件。

常见的打开⽅式如下：‘r’：只读⽅式打开⽂件（默认的⽅式），该⽂件必须已存在。

‘r+’：读写⽅式打开⽂件，打开后先读后写。

该⽂件必须已存在。

‘w’：打开后写⼊数据。

该⽂件已存在则更新；不存在则创建。

‘w+’：读写⽅式打开⽂件。

先读后写。

该⽂件已存在则更新；不存在则创建。

‘a’：在打开的⽂件末端添加数据。

⽂件不存在则创建。

‘a+’：打开⽂件后，先读⼊数据再添加数据。

⽂件不存在则创建。

另外，在这些字符串后添加⼀个“t”，如‘rt’或‘wt+’，则将该⽂件以⽂本⽅式打开；如果添加的是“b”，则以⼆进制格式打开，这也是fopen函数默认的打开⽅式。

2）关闭⽂件⽂件在进⾏完读、写等操作后，应及时关闭，以免数据丢失。

关闭⽂件⽤fclose函数，调⽤格式为：sta＝fclose(fid)说明：该函数关闭fid所表⽰的⽂件。

sta表⽰关闭⽂件操作的返回代码，若关闭成功，返回0，否则返回-1。

如果要关闭所有已打开的⽂件⽤fclose(‘all’)。

2、⼆进制⽂件的读写操作1）写⼆进制⽂件fwrite函数按照指定的数据精度将矩阵中的元素写⼊到⽂件中。

其调⽤格式为：COUNT＝fwrite（fid，A，precision）说明：其中COUNT返回所写的数据元素个数（可缺省），fid为⽂件句柄，A⽤来存放写⼊⽂件的数据，precision代表数据精度，常⽤的数据精度有：char、uchar、int、long、float、double等。

matlab中读取txt数据文件（txt文本文档）

matlab中读取txt数据⽂件（txt⽂本⽂档）根据txt⽂档不同种类介绍不同的读取数据⽅法⼀、纯数据⽂件（没有字母和中⽂，纯数字）对于这种txt⽂档，从matalb中读取就简单多了例如test.txt⽂件，内容为17.901 -1.1111 33.04517.891 -1.1286 33.04517.884 -1.1345 33.045可以在command window中输⼊load test.txt ，然后就会产⽣⼀个test的数据⽂件，内容跟test.txt中的数据⼀样；另⼀种⽅法是在file/import data....../next/finish 也可产⽣⼀个叫test的数据⽂件。

⼆、中英⽂和数据如test1.txt你好欢迎来到论坛1 11 111 11112 22 222 22223 33 333 33334 44 444 44445 55 555 5555这样的⽂件怎么读⼊数据呢？⽅法有多种，现举两个⽐较简单实⽤的。

⽅法⼀：file/import data....../next/finish>> whosName Size Bytes Classdata 5x4 160 double arraytextdata 4x1 300 cell arrayGrand total is 54 elements using 460 bytes>> datadata =1 11 111 11112 22 222 22223 33 333 33334 44 444 44445 55 555 5555>> textdatatextdata ='你好''欢迎来到''论坛'''⽅法⼆：[a1,a2,a3,a4]=textread('test1.txt','%s%s%s%s','headerlines',4)说明：%s可以是其他形式，跟读⼊的数据类型有关，⽐如这⾥也可以⽤%n,%f等。

《自动控制原理》Matlab求解控制系统数学模型实验

《自动控制原理》Matlab求解控制系统数学模型实验一、实验目的(1)熟练运用matlab软件，求解控制系统数学模型(2)掌握传递函数在matlab中的表达方法(3)掌握matlab求解拉氏变换和反变换(4)掌握matlab求系统极值点和零点判断系统稳定性二、实验仪器装配Matlab7.0的计算机三、实验原理传递函数在matlab中的表达方法控制系统的传递函数模型为：在MATLAB中，分子/分母多项式通过其系数行向量表示，即：num = [b0 b1 … bm]den = [a0 a1 … an]此时，系统的传递函数模型用tf函数生成，句法为：sys=tf(num, den)其中，sys为系统传递函数。

如：num = [1 5 0 2]; den = [2 3 15 8];则：sys=tf(num, den)输出为：Transfer function:传递函数的转换[num,den]=zp2tf(z,p,k)[z,p,k]=tf2zp(num,den)实际系统往往由多个环节通过串联、并联及反馈方式互连构成。

MATLAB提供的三个用于计算串联、并联及反馈连接形成的新系统模型的函数。

四、实验内容及步骤2、用MATLAB展求拉氏变换和反变换在MATLAB中，多项式通过系数行向量表示，系数按降序排列如要输入多项式：x4-12x3+25x+126>> p=[1 -12 0 25 126]-p = 1 -12 0 25 1263、连续系统稳定性分析的MATLAB函数roots函数：求多项式的根句法： r=roots(p)其中，r为由多项式根组成的列向量。

➢pole函数：计算系统的极点句法： p=pole(sys)其中，p为由极点组成的列向量zero函数：计算系统的零点句法： r=zero(sys) 或 [z, k]=zero(sys)其中，r为由多项式根组成的列向量。

k为零极点增益模型之增益pzmap函数：绘制零极点分布图句法： pzmap(sys) 或 [p,z] = pzmap(sys)五、实验原始数据记录与数据处理在MATLAB中，多项式通过系数行向量表示，系数按降序排列如要输入多项式：x4-12x3+25x+126>> p=[1 -12 0 25 126]-p = 1 -12 0 25 126六、实验结果与分析讨论七、结论掌握 MATLAB命令窗口的基本操作;掌握MATLAB 建立控制系统数学模型的命令及模型相互转换的方法;掌握了使用各种函数命令建立控制系统数学模型.八、实验心得体会（可略）通过该试验我们熟悉 MATLAB 实验环境，掌握 MATLAB命令窗口的基本操作;掌握MATLAB 建立控制系统数学模型的命令及模型相互转换的方法;掌握了使用各种函数命令建立控制系统数学模型:完成实验的范例题和自我实践，并记录结果;编写M文件程序，完成简单连接的模型等效传递函数，并求出相应的零极点。

数值分析拟合实验报告(3篇)

第1篇一、实验目的本次实验旨在通过数值分析方法对一组已知数据点进行拟合，掌握线性插值、多项式插值、样条插值等方法的基本原理和实现过程，并学会使用MATLAB进行数值拟合。

二、实验内容1. 线性插值线性插值是一种简单的插值方法，适用于数据点分布较为均匀的情况。

其基本原理是通过两个相邻的数据点，利用线性关系拟合出一条直线，然后通过该直线来估算未知的值。

2. 多项式插值多项式插值是一种较为精确的插值方法，通过构造一个多项式函数来逼近已知数据点。

其基本原理是利用最小二乘法求解多项式的系数，使得多项式在已知数据点上的误差最小。

3. 样条插值样条插值是一种更灵活的插值方法，通过构造一系列样条曲线来逼近已知数据点。

其基本原理是利用最小二乘法求解样条曲线的系数，使得样条曲线在已知数据点上的误差最小。

三、实验步骤1. 线性插值（1）在MATLAB中输入已知数据点，如：x = [1, 2, 3, 4, 5];y = [2, 4, 6, 8, 10];（2）使用MATLAB内置函数`linspace`生成插值点：xi = linspace(1, 5, 100);（3）使用MATLAB内置函数`interp1`进行线性插值：yi = interp1(x, y, xi, 'linear');（4）绘制插值曲线：plot(xi, yi, 'b-', x, y, 'ro');2. 多项式插值（1）在MATLAB中输入已知数据点，如：x = [1, 2, 3, 4, 5];y = [2, 4, 6, 8, 10];（2）使用MATLAB内置函数`polyfit`求解多项式系数：p = polyfit(x, y, 3);（3）使用MATLAB内置函数`polyval`进行多项式插值：yi = polyval(p, xi);（4）绘制插值曲线：plot(xi, yi, 'b-', x, y, 'ro');3. 样条插值（1）在MATLAB中输入已知数据点，如：x = [1, 2, 3, 4, 5];y = [2, 4, 6, 8, 10];（2）使用MATLAB内置函数`spline`进行样条插值：yi = spline(x, y, xi);（3）绘制插值曲线：plot(xi, yi, 'b-', x, y, 'ro');四、实验结果与分析1. 线性插值线性插值方法简单易行，但精度较低，适用于数据点分布较为均匀的情况。

实验七在MATLAB中输入输出数据文件

实验七、在MA TLAB中输入/输出数据文件
一、实验目的
掌握在MA TLAB中输入/输出数据文件的基本方法。

二、实验内容及步骤
1、已知一文本文件q.txt的内容为：Are you a student?
Do you like tea?
Are you from china?
要求从键盘输入该文本文件的名字，然后依次显示每一道题目并提示用户回答，并将用户的答案保存在文本文件qq.txt中。

(fgetl_examp.m)
2、编写程序，该程序能读取一个文本文件a.txt（该文本文件为一段文字），并能将文本文件中的小写字母转换为相应的大写字母而生成一个新的文本文件aa.txt。

(zhuanhuan.m)
3、下列程序执行后，变量b、position和c的值是多少？(pos1_examp.m) a=1:5;
fid=fopen('fdat.bin','w');
fwrite(fid,a,'int16');
status=fclose(fid);
fid=fopen('fdat.bin','r');
status=fseek(fid,6,'bof');
b=fread(fid,1,'int16');
position=ftell(fid);
status=fseek(fid,-4,'cof');
c=fread(fid,1,'int16');
status=fclose(fid);。

MATLAB如何导入excel数据

下面我就通过实例介绍一下MATLAB如何导入excel数据。

1.读取sheet1中的所有数据
我们以03版excel为例，假设excel中数据为
默认这些文本以及数字都放在sheet1中，我们将此excel选入当前工作目录（必要步骤），然后再MATLAB主窗口中输入下面命令：[NUM,TXT,R AW]=xlsread('example')，其中example是你的excel名，这里我所有
的数据都在example.xls中。

输入以上命令，回车
NUM返回的是excel中的数据，TXT输出的是文本内容，RAW输出的是未处理数据，这里因为文本占用的矩阵一共是8行6列，所以就输出了如上内容。

一般情况下，我们读取的都是excel中的数据，所以可以直接用[NUM]=xlsread('example')，只输出数据矩阵便可。

2.读取指定sheet中的数据
假如我们想读取第二个sheet中的数据，这里请注意：matlab读取exce l中的数据是按照sheet在excel中的排放顺序来的，sheet3现在排在
第二位，我们输入命令[NUM]=xlsread('example'，2),回车结果如下：
3.读取指定单元格中的数据
以sheet2中的数据为例
假设我们读取A3-D7之间的数据，那么我们输入命令[NUM]=xlsread('ex ample',2,'A3:D7'),回车，结果如下：
4.如果以上经验有用，请点击下方的有用按钮支持我的工作，谢谢！。