MATLAB第六次上机实验报告
北科大Matlab 第六次数学实验报告 (精华版)
《数学实验》报告实验名称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实验报告实验一典型环节的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上机实习实习报告上周我们进行了本学期第二个实习部分-Matlab上机实习。
本次实习一星期,共做了四个部分。
分别是:一,Matlab基本操作;二,数值数组及其运算;三,Matlab图形绘制基础;四,SIMULINK仿真基础。
下面分别介绍个部分实习内容。
实习内容一Matlab基本操作1,实习目的:·掌握Matlab的启动和退出;·熟悉Matlab的命令窗口;·熟悉其他常用窗口和工具栏;2,对本软件有了初步认识后,接下来就开始指导书要求做一些训练了在指令窗中输入下面一段程序,功能是::画出衰减振荡曲线y=(e-t/3)sin3t及其他的包络线y0=(e-t/3)。
T的取值范围是[0,4 ]。
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,通过下拉菜单File:Preferences改变设置指令窗中的字体、颜色等。
如:数据显示格式的设置:Format short 小数点后四位(缺省情况),如显示p为3.1416。
Format long 小数点后十四位,如显示p为3.14159265358979。
Formant bank 小数点后两位,如显示p为3.14。
Format short e 小数点后四位科学记数法,如显示p为3.1416e+000。
Format long e 小数点后十四五位科学记数法,如显示p为3.141592653589793e+000。
在指令窗中键入y1=2*sin(0.3*pi)/(1+sqrt(5))然后利用回调指令计算y2=2*y1*cos(0.3*pi)/(1+sqrt(5))y1 =0.5000y2 =0.18164,了解软件中的一些功能和命令。
包括熟悉历史指令窗(Command History)、熟悉当前目录浏览器(Current Diretory)、熟悉工作空间浏览器(Workspace Browser),以及一些指令如: Clear 清除当前工作区中的所有变量Clc 清除指令窗中内容(未清除当前工作区中的变量)Clf 清除图形窗口Cd 设置当前工作目录Exit,quit 退出Matlab实习内容二数值数组及其运算1,实习目的·掌握一维数组的创建和寻访;·掌握二维数组的创建和寻访;·掌握的Matlab矩阵和数组的运算;·熟悉Matlab关系操作和逻辑操作;2,实习内容和步骤(1)常量与变量变量名区分字母大小写;变量名必须以字母打头,其后可以是任意字母,数字,或下划线的组合。
matlab实验六报告
实验六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实验报告总结精选
matlab实验报告总结电气工程学院自动化102班 2012年12月21日实验一 MATLAB环境的熟悉与基本运算一、实验目的1.熟悉MATLAB开发环境2.掌握矩阵、变量、表达式的各种基本运算二、实验基本知识1.熟悉MATLAB环境MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器、文件和搜索路径浏览器。
2.掌握MATLAB常用命令变量与运算符变量命名规则如下:变量名可以由英语字母、数字和下划线组成变量名应以英文字母开头长度不大于31个区分大小写MATLAB中设置了一些特殊的变量与常量,列于下表。
MATLAB运算符,通过下面几个表来说明MATLAB的各种常用运算符表2 MATLAB算术运算符表3 MATLAB关系运算符表4 MATLAB逻辑运算符表5 MATLAB特殊运算的一维、二维数组的寻访表6 子数组访问与赋值常用的相关指令格式的基本运算表7 两种运算指令形式和实质内涵的异同表的常用函数表8 标准数组生成函数表9 数组操作函数三、实验内容1、新建一个文件夹2、启动,将该文件夹添加到MATLAB路径管理器中。
3、保存,关闭对话框4、学习使用help命令,例如在命令窗口输入help eye,然后根据帮助说明,学习使用指令eye5、学习使用clc、clear,观察command window、command history和workspace等窗口的变化结果。
6、初步程序的编写练习,新建M-file,保存,学习使用MATLAB的基本运算符、数组寻访指令、标准数组生成函数和数组操作函数。
注意:每一次M-file的修改后,都要存盘。
练习A:help rand,然后随机生成一个2×6的数组,观察command window、command history和workspace等窗口的变化结果。
学习使用clc、clear,了解其功能和作用。
答:clc是清除命令窗体内容 clear是清除工作区间输入C=1:2:20,则C表示什么?其中i=1,2,3,?,10。
MATLAB第6次实验报告
实验报告一、实验目的通过这些题目的训练,需要掌握静态数据成员以及静态成员函数的使用方法,友元函数的使用,以及掌握构造函数和析构函数,拷贝构造函数、内联成员函数、聚集的使用,熟练的使用类解决问题。
二、实验环境编译器: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实验报告(实验6)
学生实验报告一、实验目的1、掌握SIMULINK组件模块操作的基本方法;2、掌握常用连续系统模块的使用方法;3、掌握连续系统建模的基本方法;4、理解SIMULINK仿真的基本配置方法与内在运行机制。
二、实验原理Simulink是MATLAB的重要组成部分,提供建立系统模型、选择仿真参数和数值算法、启动仿真程序对该系统进行仿真、设置不同的输出方式来观察仿真结果等功能。
1、Simulink的基本模块Simulink的模块库提供了大量模块。
单击模块库浏览器中Simulink前面的“+”号,将看到Simulink模块库中包含的子模块库,单击所需要的子模块库,在右边的窗口中将看到相应的基本模块,选择所需基本模块,可用鼠标将其拖到模型编辑窗口。
同样,在模块库浏览器左侧的Simulink栏上单击鼠标右键,在弹出的快捷菜单中单击Open the ‘Simulink’ Libray 命令,将打开Simulink基本模块库窗口。
单击其中的子模块库图标,打开子模块库,找到仿真所需要的基本模块。
Simulink中几乎所有模块的参数都允许用户进行设置,只要双击要设置的模块或在模块上按鼠标右键并在弹出的快捷菜单中选择相应模块的参数设置命令就会弹出模块参数对话框。
该对话框分为两部分,上面一部分是模块功能说明,下面一部分用来进行模块参数设置。
同样,先选择要设置的模块,再在模型编辑窗口Edit菜单下选择相应模块的参数设置命令也可以打开模块参数对话框。
一个Simulink仿真模型的基本模块包括信源、信宿以及系统三个部分。
其中,信源可以是常数、正弦波、阶梯波等信号源,信宿可以是示波器、图形记录仪等,系统则是被研究系统的SIMULINK方框图。
系统、信源、信宿,可以从SIMULINK模块库中直接获得,也可以根据用户意愿用库中的模块构建而成。
2、连续系统的建模与仿真连续系统指的是可以用微分方程来描述的系统。
用于建模连续系统的模块:Simulink模块组中的Continous、Math以及Nonlinear模块库中。
MATLAB课后上机实验报告
《MATLAB程序设计与应用》实验指导书实验一 matlab 集成环境使用与运算基础1,先求下列表达式的值,然后显示matlab 工作空间的使用情况并保存全部变量。
(1)0122sin851z e =+程序:.>> z1=2*sin(85*pi/180)/(1+exp(2)) 结果: z1 =0.2375(2)22121(0.4552i z In x x +⎡⎤=+=⎢⎥-⎣⎦其中 程序:>> x=[2,1+2*i;-0.45,5];>> z2=0.5*log(x+sqrt(1+x*x)) 结果: z2 =0.7114 - 0.0253i 0.8968 + 0.3658i 0.2139 + 0.9343i 1.1541 - 0.0044i(3)0.3,9.2,8.2,...,8.2,9.2,0.3,23.0)3.0sin(23.03.03---=+++-=-a aIn a e e z a a 提示:利用冒号表达式生成a 向量,求各点函数值时用点乘运算。
程序:>> a=-3.0:0.1:30;>> z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin((a+0.3)*pi/180)+log((0.3+a)/2) 结果: z3 =1.0e+003 *Columns 1 through 40.0003 + 0.0031i 0.0003 + 0.0031i 0.0003 + 0.0031i 0.0002 + 0.0031iColumns 5 through 80.0002 + 0.0031i 0.0001 + 0.0031i 0.0001 + 0.0031i 0.0000 + 0.0031i Columns 9 through 12-0.0000 + 0.0031i -0.0001 + 0.0031i -0.0001 + 0.0031i -0.0002 + 0.0031i Columns 13 through 16-0.0003 + 0.0031i -0.0003 + 0.0031i -0.0004 + 0.0031i -0.0005 + 0.0031i Columns 17 through 20-0.0006 + 0.0031i -0.0007 + 0.0031i -0.0008 + 0.0031i -0.0009 + 0.0031i Columns 21 through 24-0.0010 + 0.0031i -0.0012 + 0.0031i -0.0014 + 0.0031i -0.0016 + 0.0031i Columns 25 through 28-0.0019 + 0.0031i -0.0023 + 0.0031i -0.0030 + 0.0031i -0.0370 Columns 29 through 32-0.0030 -0.0023 -0.0019 -0.0016 Columns 33 through 36-0.0014 -0.0012 -0.0010 -0.0009 Columns 37 through 40-0.0008 -0.0007 -0.0006 -0.0005 Columns 41 through 44-0.0004 -0.0003 -0.0003 -0.0002 Columns 45 through 48-0.0001 -0.0001 -0.0000 0.0000 Columns 49 through 520.0001 0.0001 0.0002 0.0002 Columns 53 through 560.0003 0.0003 0.0003 0.0004 Columns 57 through 600.0004 0.0005 0.0005 0.0005 Columns 61 through 640.0006 0.0006 0.0006 0.0007 Columns 65 through 680.0007 0.0007 0.0008 0.0008 Columns 69 through 720.0008 0.0008 0.0009 0.0009 Columns 73 through 760.0009 0.0010 0.0010 0.0010 Columns 77 through 800.0011 0.0011 0.0011 0.0011 Columns 81 through 840.0012 0.0012 0.0012 0.0013 Columns 85 through 880.0013 0.0013 0.0013 0.0014 Columns 89 through 920.0014 0.0014 0.0015 0.0015 Columns 93 through 960.0015 0.0016 0.0016 0.0016 Columns 97 through 1000.0017 0.0017 0.0017 0.0018 Columns 101 through 1040.0018 0.0018 0.0019 0.0019 Columns 105 through 1080.0020 0.0020 0.0020 0.0021 Columns 109 through 1120.0021 0.0022 0.0022 0.0023 Columns 113 through 1160.0023 0.0024 0.0024 0.0025 Columns 117 through 1200.0025 0.0026 0.0026 0.0027 Columns 121 through 1240.0027 0.0028 0.0029 0.0029 Columns 125 through 1280.0030 0.0031 0.0031 0.0032 Columns 129 through 1320.0033 0.0034 0.0034 0.0035 Columns 133 through 1360.0036 0.0037 0.0038 0.0039 Columns 137 through 1400.0040 0.0041 0.0042 0.0043 Columns 141 through 1440.0044 0.0045 0.0046 0.0047 Columns 145 through 1480.0049 0.0050 0.0051 0.0053 Columns 149 through 1520.0054 0.0056 0.0057 0.0059 Columns 153 through 1560.0060 0.0062 0.0064 0.0066 Columns 157 through 1600.0068 0.0069 0.0071 0.0074 Columns 161 through 1640.0076 0.0078 0.0080 0.0083 Columns 165 through 1680.0085 0.0088 0.0090 0.0093 Columns 169 through 1720.0096 0.0099 0.0102 0.0105 Columns 173 through 1760.0108 0.0112 0.0115 0.0119 Columns 177 through 1800.0123 0.0127 0.0131 0.0135 Columns 181 through 1840.0139 0.0144 0.0148 0.0153 Columns 185 through 1880.0158 0.0163 0.0168 0.0174 Columns 189 through 1920.0180 0.0185 0.0191 0.0198 Columns 193 through 1960.0204 0.0211 0.0218 0.0225 Columns 197 through 2000.0233 0.0241 0.0249 0.0257 Columns 201 through 2040.0265 0.0274 0.0284 0.0293 Columns 205 through 2080.0303 0.0313 0.0324 0.0335 Columns 209 through 2120.0346 0.0358 0.0370 0.0382 Columns 213 through 2160.0395 0.0409 0.0423 0.0437 Columns 217 through 2200.0452 0.0467 0.0483 0.0500 Columns 221 through 2240.0517 0.0534 0.0552 0.0571 Columns 225 through 2280.0591 0.0611 0.0632 0.0654 Columns 229 through 2320.0676 0.0699 0.0723 0.0748 Columns 233 through 2360.0773 0.0800 0.0827 0.0856 Columns 237 through 2400.0885 0.0915 0.0947 0.0979 Columns 241 through 2440.1013 0.1047 0.1083 0.1121 Columns 245 through 2480.1159 0.1199 0.1240 0.1282 Columns 249 through 2520.1326 0.1372 0.1419 0.1467 Columns 253 through 2560.1518 0.1570 0.1624 0.1679 Columns 257 through 2600.1737 0.1796 0.1858 0.1921 Columns 261 through 2640.1987 0.2055 0.2125 0.2198 Columns 265 through 2680.2273 0.2351 0.2431 0.2514 Columns 269 through 2720.2600 0.2689 0.2781 0.2876 Columns 273 through 2760.2974 0.3076 0.3180 0.3289 Columns 277 through 2800.3401 0.3517 0.3637 0.3761 Columns 281 through 2840.3889 0.4021 0.4158 0.4299 Columns 285 through 2880.4446 0.4597 0.4753 0.4915 Columns 289 through 2920.5082 0.5254 0.5433 0.5617 Columns 293 through 2960.5807 0.6004 0.6208 0.6418 Columns 297 through 3000.6636 0.6861 0.7093 0.7333 Columns 301 through 3040.7581 0.7838 0.8103 0.8377 Columns 305 through 3080.8660 0.8952 0.9254 0.9567 Columns 309 through 3120.9890 1.0223 1.0568 1.0924 Columns 313 through 3161.1292 1.1673 1.2066 1.2472Columns 317 through 3201.2892 1.3326 1.3774 1.4237Columns 321 through 3241.4715 1.5210 1.5721 1.6249Columns 325 through 3281.6794 1.7357 1.7940 1.8541Columns 329 through 3311.9163 1.98052.0468(4)⎪⎩⎪⎨⎧=<≤<≤<≤+--=5.2:5.0:0,322110,121,2224t t t t t t t t z 其中提示:用逻辑表达式求分段函数值。
MATLAB上机实验实验报告
MATLAB上机实验一一、实验目的初步熟悉MATLAB 工作环境,熟悉命令窗口,学会使用帮助窗口查找帮助信息。
命令窗口二、实验内容(1) 熟悉MATLAB 平台的工作环境。
(2) 熟悉MATLAB 的5 个工作窗口。
(3) MATLAB 的优先搜索顺序。
三、实验步骤1. 熟悉MATLAB 的5 个基本窗口①Command Window (命令窗口)②Workspace (工作空间窗口)—③Command History (命令历史记录窗口)④Current Directory (当前目录窗口)⑤Help Window (帮助窗口)(1) 命令窗口(Command Window)。
在命令窗口中依次输入以下命令:>>x=1>> y=[1 2 34 5 67 8 9];>> z1=[1:10],z2=[1:2:5];>> w=linspace(1,10,10);>> t1=ones(3),t2=ones(1,3),t3=ones(3,1)>> t4=ones(3),t4=eye(4)x =1z1 =1 2 3 4 5 6 7 8 9 10 t1 =1 1 11 1 11 1t2 =1 1 1t3 =111t4 =1 1 11 1 11 1 1t4 =1 0 0 00 1 0 00 0 1 00 0 0 1思考题:①变量如何声明,变量名须遵守什么规则、是否区分大小写。
答:(1)变量声明1.局部变量每个函数都有自己的局部变量,这些变量只能在定义它的函数内部使用。
当函数运行时,局部变量保存在函数的工作空间中,一旦函数退出,这些局部变量将不复存在。
脚本(没有输入输出参数,由一系列MATLAB命令组成的M文件)没有单独的工作空间,只能共享调用者的工作空间。
当从命令行调用,脚本变量存在基本工作空间中;当从函数调用,脚本变量存在函数空间中。
2.全局变量在函数或基本工作空间内,用global声明的变量为全局变量。
Matlab上机报告
《MATLAB 课程》课程实习报告(一) 实习日期:2016-9-13学院:数学与统计学院学院 专业:数学与应用数学 姓名:王子璇 学号:201413230061, 先求下列表达式的值,然后显示MATLAB 工作空间的使用情况,并保存全部变量。
(1)21185sin 2ez +=代码:>> z1=2*sin(85*pi/180)/(1+exp(2))截图:(3)()0.3,9.2,8.2,...,8.29.2,0.3,233.0ln 3.0sin 23.03.03---=+++-=-a a a e e z a a 提示:利用冒号表达式生成a 向量,求各点的函数值时用点乘运算。
代码:>> a=-3.0:0.1:3.0;>> z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)截图:(4)5.2:5.0:0其中,32,1221,110,2224=⎪⎩⎪⎨⎧<≤+-<≤-<≤=t t t t t t t t z提示:用逻辑表达式求分段函数值。
代码:>> t=0:0.5:2.5;>> z4=t.^2.*(t>=0&t<1)+(t.^2-1).*(t>=1&t<2)+(t.^2-2*t+1).*(t>=2&t<3)截图:2, 设有矩阵A 和B⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=25242322212019181716151413121110987654321A ,⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡--=11134079423096171603B (1) 求它们的乘积C 。
(2) 将矩阵C 的右下角23⨯子矩阵赋给D 。
(3) 查看MATLAB 工作空间的使用情况。
代码:(1)>> A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25]>> B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]>> C=A*B(2)>> F=size(C);>> D=C(F(1)-2:F(1),F(2)-1:F(2));>> whos截图:(1)(2)3,下面是一个线性方程组:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡52.067.095.061141514111121321x x x(1) 求方程的解。
MATLAB技术上机实践报告
目录上机实践一MATLAB在“高等数学”中的应用练习(4学时)(1)上机练习——MATLAB解常微分方程(SJ101.m)(2)上机练习——MATLAB二元函数作图(SJ102.m)(3)上机编程——微分运算和积分运算(SJ103.m)(4)上机编程——线性方程组的数值解法(SJ104.m)上机实践二MATLAB在“力学”中的应用练习(4学时)(1)上机练习——质点直线运动的位置、速度和加速度(SJ201.m)(2)上机练习——轻质杆的斜抛运动(SJ202.m)(3)上机编程——质点的抛体运动(SJ203.m)(4)上机编程——物体在恒力作用下的直线运动(SJ204.m)上机实践三MATLAB在“电磁学”中的应用练习(4学时)(1)上机练习——电磁波的发射方法一(SJ301.m)(2)上机练习——电磁波的发射方法二(SJ302.m)(3)上机编程——静电场场强分布计算(SJ303.m)(4)上机编程——带电粒子在电场中的运动(SJ304.m)上机实践四MATLAB在“振动和波动”中的应用练习(4学时)(1)上机练习——旋转矢量法(SJ401.m)(2)上机练习——用惠更斯作图法确定波阵面(SJ402.m)(3)上机编程——同方向简谐振动的合成(SJ403.m)(4)上机编程——简谐波横波的演示(SJ404.m)上机实践五MATLAB在“量子物理”中的应用练习(4学时)(1)上机练习——黑体辐射随波长的变化规律(SJ501.m)(2)上机练习——势垒和隧道效应(SJ502.m)(3)上机编程——双缝干涉图样的动画模拟(SJ503.m)上机实践六MATLAB在“电子信息技术”中的应用练习(4学时)(1)上机练习——含受控源的电阻电路计算(SJ601.m)(2)上机练习——方波分解为多次正弦波之和(SJ602.m)(3)上机编程——电阻电路的计算(SJ603.m)(4)上机编程——连续信号的MATLAB描述(SJ604.m)上机实践一 MATLAB 在“高等数学”中的应用练习(4学时)【上机实践目的】使用MATLAB 软件编写相应的M 文件进行计算模拟高等数学中的一些基本问题。
MATLAB上机实验报告
MATLAB上机实验报告实验一、实验目的--------------------------------------------------------------------------------------------------------3二、实验内容--------------------------------------------------------------------------------------------------------3三、实验步骤--------------------------------------------------------------------------------------------------------31. M 文件的编辑----------------------------------------------------------------------------------------------32. 程序流程控制结构---------------------------------------------------------------------------------------31) for 循环结构-----------------------------------------------------------------------------------------42) while 循环结构----------------------------------------------------------------------------------------5练习-----------------------------------------------------------------------------------------------------------53)if-else-end分支结构---------------------------------------------------------------------------------54)switc-case结构----------------------------------------------------------------------------------------63. 子函数和参数传递------------------------------------------------------------------------------------------6练习-----------------------------------------------------------------------------------------------------------74. 局部变量和全局变量--------------------------------------------------------------------------------------8 MATLAB上机练习6.1 M 文件--------------------------------------------------------------------------------------------8(1)脚本文件-----------------------------------------------------------------------------------------------------8(2)函数文件-----------------------------------------------------------------------------------------------------9(3)函数调用-----------------------------------------------------------------------------------------------------9 6.2 MATLAB 的程序控制结构( 1) for 循环----------------------------------------------------------------------------------------------------9(2) while 循环----------------------------------------------------------------------------------------------------12(3)if语句---------------------------------------------------------------------------------------------------------12(4)switch 语句-------------------------------------------------------------------------------------------------13(5)try 语句-----------------------------------------------------------------------------------------------------14(6)程序流的控制--------------------------------------------------------------------------------------------14 6.3 数据的输入与输出-----------------------------------------------------------------------------------14(1)键盘输入语句------------------------------------------------------------------------------------------14(2)屏幕输出语句---------------------------------------------------------------------------------------------14 6.4 MATLAB文件操作(1)fopen 语句-----------------------------------------------------------------------------------------------------15(2) fclose 语句-----------------------------------------------------------------------------------------------------15 6.5 面向对象编程--------------------------------------------------------------------------------------------15(1)创建类目录----------------------------------------------------------------------------------------------------15(2) 建立类的数据结构-----------------------------------------------------------------------------------------15(3)创建类的基本方法---------------------------------------------------------------------------------------------16 6.6 MATLAB 程序优化-----------------------------------------------------------------------------------------17 6.7程序调试------------------------------------------------------------------------------------------------------17 6.8 习题----------------------------------------------------------------------------------------------------------17实验五MATLAB 程序设计一、实验目的掌握MATLAB 程序设计的主要方法,熟练编写MATLAB 函数。
matlab上机实习报告
Matlab上机实习报告学院班级学号姓名指导老师:2011年12月25日第十六周我们进行了本学期第二个实习部分-Matlab上机实习。
时间过得很快,一周的时间转眼间就过去了。
通过实习我们对Matlab有了基本的了解。
本次实习一星期,共有五个部分。
分别是:一,Matlab基本操作;二,数值数组及其运算;三,Matlab图形绘制基础;四,SIMULINK仿真基础;五,句柄图形和图形用户界面制作。
实习内容一Matlab基本操作1,实习目的:·掌握Matlab的启动和退出;·熟悉Matlab的命令窗口;·熟悉其他常用窗口和工具栏;2,对本软件有了初步认识后,接下来就开始按照指导书要求做一些训练了在指令窗中输入下面一段程序,功能是::画出衰减振荡曲线y=(e-t/3)sin3t及其他的包络线y0=(e-t/3)。
T的取值范围是。
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,通过下拉菜单File:Preferences改变设置指令窗中的字体、颜色等。
如:数据显示格式的设置:Format short 小数点后四位(缺省情况),如显示p为3.1416。
Format long 小数点后十四位,如显示p为3.14159265358979。
Formant bank 小数点后两位,如显示p为3.14。
Format short e 小数点后四位科学记数法,如显示p为3.1416e+000。
Format long e 小数点后十四五位科学记数法,如显示p为3.141592653589793e+000。
在指令窗中键入y1=2*sin(0.3*pi)/(1+sqrt(5))然后利用回调指令计算y2=2*y1*cos(0.3*pi)/(1+sqrt(5))y1 =0.5000y2 =0.18164,了解软件中的一些功能和命令。
matlab实验内容答案解析
实验报告说明:matlab 课程实验需撰写8个实验报告,每个实验报告内容写每次实验内容中标号呈黑体大号字显示的题目。
第一次实验内容:实验一 MATLAB 运算基础一、实验目的1.熟悉启动和退出MATLAB 的方法。
2.熟悉MATLAB 命令窗口的组成。
3.掌握建立矩阵的方法。
4.掌握MATLAB 各种表达式的书写规则以及常用函数的使用。
二、实验内容1.先求下列表达式的值,然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1)22sin 8511z 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:0.5:2.5 提示:用逻辑表达式求分段函数值。
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 工作空间使用情况。
4.完成下列操作:(1)求[100,999]之间能被21整除的数的个数。
(完整word)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/2 3.两个函数之间的操作求和(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上机实验报告实验1MATLAB初步1. 举例说明format命令的功能;Format命令⽤来调整数据的显⽰格式。
>> format long %15位定点表⽰,若改成short则为5位定点表⽰>> 2^ans =2. 说明clc,clf和clear命令的区别;clc即clear command(清屏)clear清理内存所有变量clf即clear figure(清理图形窗⼝3. 举例说明help命令的功能;在help命令后加⼊所要查询的函数或类名,将会返回所查询的对象的使⽤说明。
例如:>> help cosCOS Cosine.COS(X) is the cosine of the elements of X.See also acos, cosd.Overloaded functions or methods (ones with the same name in other directories)help sym/Reference page in Help browserdoc cos4. 如何观察和调⽤历史命令;在MTALAB界⾯的左下⾓有个Command History窗⼝,⾥⾯可以看到历史命令,如果双击该条历史命令即可实现调⽤。
5. 如何设置当前路径。
在MATLAB界⾯上⽅的⼯具栏有个Current Folder(当前⽬录),更改该⽬录就可以设置当前路径。
实验2矩阵⽣成与元素标识1. 矩阵的创建=32943253036120031A ,=757065605550454035302520151051B (尽可能⽤最简单⽅法);创建矩阵A : >> A=[1 3 0 0;2 61 3 0;3 5 2 3;4 9 2 3]创建矩阵B : >>B=[1 5 10 15;20 25 30 35;40 45 50 55;60 65 70 75]2. ⽣成⼀个单位矩阵C ,其结构与B 相同;>> eye(size(B))3.找出矩阵A 中所有⼤于8的元素位置(⾏号和列号);>> C=find(B>35 & B<65)>> B(C)=B(C)+95. 提取矩阵B 的主对⾓组元素。
matlab上机实验报告
上海电力学院数值计算方法上机实习报告院系:能源与机械工程学院专业年级:动力机械及工程2012级学生姓名:***学号:ys**********指导教师:***2012年12月26日数值计算方法上机实习题1. 设⎰+=105dx xx I nn , (1) 由递推公式nI I n n 151+-=-,从0I 的几个近似值出发,计算20I ; 解:I 0=1015x d x +⎰=0.1823计算I20编辑matlab 命令如下:I=0.1823 for n=1:1:20, I=-5*I+1/n;fprintf('%.1d %.4f\n',n ,I); end结果:(2) 粗糙估计20I ,用nI I n n 51511+-=-,计算0I ; 解:I 20= 20105xx d x +⎰ 使用复合中点公式进行积分,相应的matlab 程序如下:I=0;for h=0:0.001:1,m=h+0.0005;I=I+0.001*m^20/(5+m);fprintf('%.1d %.4f\n',m ,I);enddisp(I);for k=1:20,n=21-k;I=0.2*(1/n-I);fprintf('%.1d %.4f\n',n ,I);enddisp(I)结果:程序结束时输出两个I值,第一个表示I20,第二个表示I0;分别为I20=0.0082I0=0.1823(3)分析结果的可靠性及产生此现象的原因(重点分析原因)。
从上述计算中分析得到如果先得到I0,再从I0由递推公式得到I20,I20结果跟精确值相比误差很大;如果先估算I20,在从I20有递推公式得到I0,I0的结果跟精确值相比近似相等。
原因分析:如果从I 0推I 20的近似值,需要用到递推公式I n =-5I n-1+1/n ,I0本身结果是有误差的;经过递推公式计算20次,就等于误差被认为的放大5的20次方倍,所以得到的I 20与其精确值相差甚远。
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课程实验报告学院:资源与安全工程学院专业班级:城地1503班学号: 020******* 姓名:沈超锋指导教师:向淑晃2016年12月目录实验一 MATLAB工作环境熟悉及简单命令的执行 (2)一、实验目的 (2)二、实验内容 (2)三、实验步骤: (2)实验二 MATLAB的数值计算功能(一) (6)一、实验目的 (6)二、实验内容 (6)三、实验步骤 (6)实验三 MATLAB的数值计算功能(二) (11)一、实验目的: (11)二、实验内容 (11)三、实验步骤 (11)实验四 MATLAB程序设计入门 (15)一、实验目的 (15)二、实验内容 (15)三、实验步骤 (15)实验五 MATLAB的符号运算功能 (19)一、实验目的 (19)二、实验内容 (19)三、实验步骤 (19)实验六 MATLAB的可视化功能 (22)一、实验目的 (22)二、实验内容 (22)三、实验步骤 (22)数学建模MATLAB上机实验报告实验一MATLAB工作环境熟悉及简单命令的执行第九周10月29日一、实验目的:熟悉MATLAB的工作环境,学会使用MATLAB进行一些简单的运算。
二、实验内容:MATLAB的启动和退出,熟悉MATLAB的桌面(desktop),包括菜单(Menu),工具条(toolbar),命令窗口(Command Window),历史命令窗口、工作空间(Workspace)等;完成一些基本的矩阵操作;学习使用在线帮助系统。
三、实验步骤:(一)MATLAB的工作环境1、启动MATLAB,熟悉MATLAB的桌面。
MATLAB的工作环境主要由命令窗(Command Window)、若干个图形窗(Figure window)、文本编辑器(File Editor)和文件管理器组成。
各视窗之间的切换可用快捷键Alt+Tab。
2、命令窗命令窗是人们与MATLAB作人机对话的主要环境,通过它可以键入MATLAB的各种命令并读出相应的结果。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
MATLAB第六次上机实验报告0210901班学号2009210811 姓名:李贤凤Quiz 3.1>> a=20;>> b=-2;>> c=0;>> d=1;>> a>bans =1>> b>dans =>> a>b&c>dans =>> a==bans =>> a&b>cans =>> ~~bans =1>> a=2;>> b=[1 -2;-0 10];>> c=[0 1;2 0];>> d=[-2 1 2;0 1 0];>> ~(a>b)ans =0 00 1>> a>c&b>cans =1 00 1>> c<=d??? Error using ==> <=Matrix dimensions must agree.>> a=2;>> b=3;>> c=10;>> d=0;>> a*b^2>a*cans =>> d|b>aans =1>> (d|b)>aans =>> a=20;>> b=-2;>> c=0;>> d='Test';>> isinf(a/b)ans =>> isinf(a/c)Warning: Divide by zero.(Type "warning off MATLAB:divideByZero" to suppress this warning.)ans =1>> a>b&ischar(d)ans =1>> isempty(c)ans =Quiz3.21.% Script file : Sqrt_x.m% Purpose:% This program is used to calculate the square root of a randem number% Record of revisions:% Date Programmer Description of change % === ======== ================ % 10/22/2010 lixianfeng Original code% Define variables :% x --x is a randem number% y --y is a function of the sqrt(x)% Prompt the user for the value xx=input('Enter the number of x:');% Calculate the function sqrt(x)if x>=0y=sqrt(x);fprintf('The square root of x is:%f',y);elsedisp('It is an error message about the square root function ');disp('The value of sqrt is 0');endexample:test aEnter the number of x:5The square root of x is: 2.236068Enter the number of x:-2It is an error message about the square root functionThe value of sqrt is 02. % Script file : Scores.m% Purpose:% This program is used to calculate the numerator / denominator% Record of revisions:% Date Programmer Description of change% ==== ========== =================== % 10/22/2010 Li xianfeng Original code% Define variables :% fun --fun is calculated as numerator / denominator% Prompt the user for the value numerator and denominatornumerator=input('Enter the numerator:');denominator=input('Enter the denominator:');% Calculate the function numerator / denominatorif abs(denominator)>=1.0E-300fun=numerator / denominator ;fprintf('The numerator / denominator is:%f',fun);elsedisp('Divide by 0 error ');endexample: >> test cEnter the numerator:4Enter the denominator:6The numerator / denominator is0.666667>> test cEnter the numerator:3Enter the denominator:-5E-400Divide by 0 error3. % Script file :cost.m% Purpose:% This program is used to calculate the cost per mile for a rented vehicle% Record of revisions:% Date Programmer Description of change% ==== ========== =================== % 10/22/2010 Li xianfeng Original code% Define variables :% d --d is the distance% aver --the average cost per mile% Prompt the user for the value distanced=input('Enter the distance:');% Calculate the functionif d>300disp('The cost per mile for a tented vehicleis $0.20 ');elseif d>100aver=(100*0.50+(d-100)*0.30)/d;disp(['The average cost per mile is ' num2str(aver)]);elseif d>=0disp('The average cost per mile is $0.50');elsedisp('The distance should be 0 or positive number ');endexample: Enter the distance:325The cost per mile for a tented vehicleis $0.20>> test3Enter the distance:158The average cost per mile is 0.42658>> test3Enter the distance:38The average cost per mile is $0.50>> test3Enter the distance:-23The distance should be 0 or positive number4. if volts>125disp('WARNING:High voltage on line.');if volts<105disp('WARNING:Low voltage on line.');elsedisp('Line voltage is within tolerances.');endAns:There is no …Input‟ sentence and the second “if” should be replaced by “elseif” or adda “end” at the end of program..5. color = 'yellow';switch (color)case 'red',disp('stop now!');case 'yellow',disp('prepare to stop.');case 'green',disp('proceed through intersection.');otherwise,disp('illegal color encountered.');endAns: prepare to stop.6. if temperature>100disp('human body temperature exceeded.');elseif temperature>37disp('boiling point of water exceeded.');endAns: Undefined function or variable 'temperature'.3.8 Exercises:3.1% Script file :tan(theta).m% Purpose:% This program is used to evaluate the tangent function% tan_theta=sin(theta)/cos(theta)% Record of revisions:% Date Programmer Description of change% ==== ========== ================% 10/22/2010 Li xianfeng Original code% Define variables :% theta --the angle% tan --the tangent function% sin --the tangent function% cos --the tangent function% Prompt the user for the value thetatheta=input('Enter the theta:');% Calculate the functionif cos(theta)<1e-20disp('error message');elsetan_theta=sin(theta)/cos(theta);disp(['The tan_theta is ' num2str(tan_theta)]);endexample: Enter the theta:5error message>> test3_1Enter the theta:0.05The tan_theta is 0.0500423.2 Ans:The logical order is wrong,the last”elseif” shoude be corrected by “else”and there lack aintput sentence.if temp>103.0disp('Temperature dangerously high');elseif temp>99.5disp('Temperature slightly high');elseif temp>97.5disp('Temperature normal');else temp<97.5disp('Temperature below normal');endexample: enter the Temperature:88temp =88Temperature below normalenter the Temperature:98temp =98Temperature normalenter the Temperature:101temp =101Temperature slightly highenter the Temperature:121temp =121Temperature dangerously high3.3 % Script file :cost.m% Purpose:% This program is used to calculate the cost of sending a package% Record of revisions:% Date Programmer Description of change% ==== ========== =================== % 10/22/2010 Li xianfeng Original code% Define variables :% weight --the weight of package% cost --the cost of sending a package% Prompt the user for the value weightweight=input('Enter the weight:');% Calculate the functionif weight>100disp('the package is not accepted!!!');elseif weight>70cost=10+3.75*68+(weight-70);disp(['The cost of sending a package is $' num2str(cost)]);elseif weight>2cost=10+3.75*(weight-2);disp(['The cost of sending a package is $' num2str(cost)]);elseif weight>0cost=10;disp('The cost of sending a package is $10.00' );elsedisp('the weight can‟t be a negative number.')end3.4 % Script file :funxy.m% Purpose:% This program is solving the function f(x,y) for a user-specified x and y% ,where f(x,y) is defined as:% -% | x+y x>=0 and y>=0% f(x,y)= | x+y^2 x>=0 and y<0% | x^2+y x<0 and y>=0% | x^2+y^2 x<0 and y<0% -% Record of revisions:% Date Programmer Description of change % ==== ========== =============== % 10/22/2010 LI xianfeng Original code% Define variables :% x --the first independent variable% y --the second independent variable% fun --resulting function% Prompt the user for the value xx=input('Enter the x coefficient:');y=input('Enter the y coefficient:');% Calculate the function fun(x,y) based on the signs of x and y.if x>=0 & y>=0fun=x+y;elseif x>=0 & y<0fun=x+y^2 ;elseif x<0 & y>=0fun=x^2+y;elsefun=x^2+y^2;endendend% write the value of the functiondisp(['the value of the function is ' num2str(fun)]);3.5 % Script file :lnx.m% Purpose:% This program is solving the function ln(1/(1-x)) for a user-specified x and y% Record of revisions:% Date Programmer Description of change% ==== ========== =================== % 10/22/2010 Li xianfeng Original code% Define variables :% x --the independent variable% y --the function of ln(1/(1-x))% ln --the natural logarithm,the logarithm to the base e.% Prompt the user for the value xx=input('Enter the x :');% Calculate the function ln(1/(1-x)) based upon the signs of x .if x<1y=log(1/(1-x));disp(['the value of ln(1/(1-x)) is;' num2str(y)]);elsedisp('the variable must be positive number.');endexample: Enter the x :-9the value of ln(1/(1-x)) is;-2.3026Enter the x :8the variable must be positive number。