《matlab基础及其应用教程》笔记
MATLAB基础教程笔记

MATLAB7.0基础教程清华大学笔记第一章MATLAB 7.0简介第二章MATLAB 7.0的安装和用户界面1.3.4.命令窗口查询函数查找具体的函数帮助help查找含有该字段的函数帮助lookfor第三章基本使用方法1.续行符…2.常用的操作命令Cd , clc , clear , clf , diary , dir , disp , echo , hold , load , pack , path , quit , save , type3.MATLAB 7.0 的常用常量Ans , beep , pi , eps , inf , nan , nargin , nargout , varagin , varagout , realmin , realmax , bitmax4.常用的函数Exp , log , log10 , fix , floor , ceil , round , rem , mod , sign , pow2 , sqrt , abs5.rem/mod(X,Y)当X,Y符号相同的时候,这两者一样,当符号不相同的时候才有差别,具体表现在rem结果的符号与被除数相同;mod结果的符号与除数相同6.表达式不能与复数I , j直接相连,只有数字能与复数I , j直接相连第四章数值计算功能1.等差向量的生成:(1)冒泡生成法,如:ans=[m:x:n] x代表步长(2) linspace 函数,如:ans=linspace(m,n,y), y代表向量个数,默认为100 2.向量与数的四则运算加法A+b 减法A-b3. 向量求点积:dot (向量维数要一致)向量求叉积:cross (向量维数要为3)向量求混合积:ans=dot(x1,cross(x2,x3)) (dot和cross的顺序不要搞错了)4. 矩阵的除法x=B\A 就是A*x=B的解,但是要求B与A行数一致x=B/A 就是x*A=B的解,但是要求B与A列数一致。
MATLAB笔记

MATLAB“关键符(词)”的说明A.1 A.1MATLAB的标点及符号A.1.1 A.1.1算术运算符Arithmetic operators+ 加;正号- 减;负号* 矩阵乘.* 数组乘\ 矩阵左除/ 矩阵右除.\ 数组左除./ 数组右除^ 矩阵幂.^ 数组幂A.1.2 A.1.2关系运算符 Relational operators== 等于~= 不等于< 小于> 大于<= 小于等于>= 大于等于A.1.3 A.1.3逻辑运算符 Logical operators& 逻辑与| 逻辑或~ 逻辑非A.1.4 A.1.4特殊符号Special characters, 逗号; 分号空格. 小数点号、构架域号: 冒号... 续行号' ' 单引号共轭转置号.' 转置号= 赋值号_ 下连符! 调用DOS操作指令号( ) 圆括号[ ] 方括号[ ] 空阵{ } 花括号% 注释号%# 编译注记A.2 A.2MATLAB的函数及指令Functions and CommandsA.2.1 A.2.1 A aabs 绝对值、模、字符的ASCII码值acos 反余弦acosh 反双曲余弦acot 反余切acoth 反双曲余切acsc 反余割acsch 反双曲余割align 启动图形对象几何位置排列工具all 所有元素非零为真angle 相角ans 表达式计算结果的缺省变量名any 所有元素非全零为真area 面域图argnames 函数M文件宗量名asec 反正割asech 反双曲正割asin 反正弦asinh 反双曲正弦assignin 向变量赋值atan 反正切atan2 四象限反正切atanh 反双曲正切autumn 红黄调秋色图阵 axes 创建轴对象的低层指令 axis 控制轴刻度和风格的高层指令 A.2.2 A.2.2 B b bar二维直方图 bar3 三维直方图bar3h 三维水平直方图 barh 二维水平直方图 base2decX 进制转换为十进制bin2dec 二进制转换为十进制blanks 创建空格串 bone 蓝色调黑白色图阵 box 框状坐标轴breakwhile 或for 环中断指令brighten 亮度控制A.2.3 A.2.3 C ccapture (5.3版以前)捕获当前图形 cart2pol 直角坐标变为极或柱坐标 cart2sph 直角坐标变为球坐标 cat 串接成高维数组 caxis 色标尺刻度 cd 指定当前目录 cdedit 启动用户菜单、控件回调函数设计工具 cdf2rdf 复数特征值对角阵转为实数块对角阵 ceil 向正无穷取整 cell 创建元胞数组 cell2struct 元胞数组转换为构架数组 celldisp 显示元胞数组内容 cellplot 元胞数组内部结构图示 char 把数值、符号、内联类对象转换为字符对象 chi2cdf2χ分布累计概率函数chi2inv 2χ分布逆累计概率函数chi2pdf 2χ分布概率密度函数 chi2rnd 2χ分布随机数发生器 chol Cholesky 分解clabel 等位线标识 cla 清除当前轴class 获知对象类别或创建对象 clc 清除指令窗 clear清除内存变量和函数clf 清除图对象 clock 时钟 collect 符号计算中同类项合并 colmmd 列最小度排序 colorbar 色标尺colorcube 三浓淡多彩交叉色图矩阵colordef 设置色彩缺省值 colormap 色图colspace 列空间的基close关闭指定窗口colperm 列排序置换向量comet彗星状轨迹图comet3 三维彗星轨迹图 compass 射线图 compose 求复合函数cond(逆)条件数condeig 计算特征值、特征向量同时给出条件数 condest 范 –1条件数估计 conj 复数共轭 contour等位线contourf 填色等位线contour3 三维等位线contourslice 四维切片等位线图conv 多项式乘、卷积cool 青紫调冷色图copper 古铜调色图cos 余弦cosh 双曲余弦cot 余切coth 双曲余切cplxpair 复数共轭成对排列csc 余割csch 双曲余割cumsum 元素累计和cumtrapz 累计梯形积分cylinder 创建圆柱A.2.4 A.2.4 D ddblquad 二重数值积分deal 分配宗量deblank 删去串尾部的空格符dec2base 十进制转换为X进制dec2bin 十进制转换为二进制dec2hex 十进制转换为十六进制deconv 多项式除、解卷delaunay Delaunay 三角剖分del2 离散Laplac ian差分demo MATLAB演示det 行列式diag 矩阵对角元素提取、创建对角阵diary MATLAB指令窗文本内容记录diff 数值差分、符号微分digits 符号计算中设置符号数值的精度dir 目录列表disp 显示数组display 显示对象内容的重载函数dlinmod 离散系统的线性化模型dmperm 矩阵Dulmage-Mendelsohn 分解dos 执行DOS 指令并返回结果double 把其他类型对象转换为双精度数值drawnow 更新事件队列强迫MATLAB刷新屏幕dsolve 符号计算解微分方程A.2.5 A.2.5 E eecho M文件被执行指令的显示edit 启动M文件编辑器eig 求特征值和特征向量eigs 求指定的几个特征值end 控制流FOR等结构体的结尾、数组每维最后元素下标eps 浮点相对精度error 显示出错信息并中断执行errortrap 错误发生后程序是否继续执行的控制erf 误差函数erfc 误差补函数erfcx 刻度误差补函数erfinv 逆误差函数errorbar 带误差限的曲线图etime 计算时间间隔etree 消去树etreeplot 画消去树eval 串演算指令evalin 跨空间串演算指令exist 检查变量或函数是否已定义exit 退出MATLAB环境exp 指数函数expand 符号计算中的展开操作expint 指数积分函数expm 常用矩阵指数函数expm1 Pade法求矩阵指数expm2 Taylor法求矩阵指数expm3 特征值分解法求矩阵指数eye 单位阵ezcontour 画等位线的简捷指令ezcontourf 画填色等位线的简捷指令ezgraph3 画表面图的通用简捷指令ezmesh 画网线图的简捷指令ezmeshc 画带等位线的网线图的简捷指令ezplot 画二维曲线的简捷指令ezplot3 画三维曲线的简捷指令ezpolar 画极坐标图的简捷指令ezsurf 画表面图的简捷指令ezsurfc 画带等位线的表面图的简捷指令A.2.6 A.2.6 F ffactor 符号计算的因式分解feather 羽毛图feedback 反馈连接feval 执行由串指定的函数fft 离散Fourier变换fft2 二维离散Fourier变换fftn 高维离散Fourier变换fftshift 直流分量对中的谱fieldnames 构架域名figure 创建图形窗fill 二维多边形填色图fill3 三维多边形填色图find 寻找非零元素下标findobj 寻找具有指定属性的对象图柄findstr 寻找短串的起始字符下标findsym 机器确定内存中的符号变量finverse 符号计算中求反函数fix 向零取整flag 红白蓝黑交错色图阵fliplr 矩阵的左右翻转flipud 矩阵的上下翻转flipdim 矩阵沿指定维翻转floor 向负无穷取整flops 浮点运算次数flow MATLAB提供的演示数据fmin 求单变量非线性函数极小值点(旧版)fminbnd 求单变量非线性函数极小值点fmins 单纯形法求多变量函数极小值点(旧版)fminunc 拟牛顿法求多变量函数极小值点fminsearch 单纯形法求多变量函数极小值点fnder 对样条函数求导fnint 利用样条函数求积分fnval 计算样条函数区间内任意一点的值fnplt 绘制样条函数图形fopen 打开外部文件for 构成for环用format 设置输出格式fourier Fourier 变换fplot 返函绘图指令fprintf 设置显示格式fread 从文件读二进制数据fsolve 求多元函数的零点full 把稀疏矩阵转换为非稀疏阵funm 计算一般矩阵函数funtool 函数计算器图形用户界面fzero 求单变量非线性函数的零点A.2.7 A.2.7 G ggallery 特殊测试矩阵gamma Γ函数gammainc 不完全Γ函数gammaln Γ函数的对数gca 获得当前轴句柄gcbo 获得正执行“回调”的对象句柄gcf 获得当前图对象句柄gco 获得当前对象句柄geomean 几何平均值get 获知对象属性getfield 获知构架数组的域getframe 获取影片的帧画面ginput 从图形窗获取数据global 定义全局变量gplot 依图论法则画图gradient 近似梯度gray 黑白灰度grid 画分格线griddata 规则化数据和曲面拟合gtext 由鼠标放置注释文字guide 启动图形用户界面交互设计工具A.2.8 A.2.8 H hharmmean 调和平均值help 在线帮助helpwin 交互式在线帮助helpdesk 打开超文本形式用户指南hex2dec 十六进制转换为十进制hex2num 十六进制转换为浮点数hidden 透视和消隐开关hilb Hilbert矩阵hist 频数计算或频数直方图histc 端点定位频数直方图histfit 带正态拟合的频数直方图hold 当前图上重画的切换开关horner 分解成嵌套形式hot 黑红黄白色图hsv 饱和色图A.2.9 A.2.9 I iif-else-elseif 条件分支结构ifft 离散Fourier反变换ifft2 二维离散Fourier反变换ifftn 高维离散Fourier反变换ifftshift 直流分量对中的谱的反操作ifourier Fourier反变换i, j 缺省的“虚单元”变量ilaplace Laplace反变换imag 复数虚部image 显示图象imagesc 显示亮度图象imfinfo 获取图形文件信息imread 从文件读取图象imwrite 把图象写成文件ind2sub 单下标转变为多下标inf 无穷大info MathWorks公司网点地址inline 构造内联函数对象inmem 列出内存中的函数名input 提示用户输入inputname 输入宗量名int 符号积分int2str 把整数数组转换为串数组interp1 一维插值interp2 二维插值interp3 三维插值interpn N维插值interpft 利用FFT插值intro MATLAB自带的入门引导inv 求矩阵逆invhilb Hilbert矩阵的准确逆ipermute 广义反转置isa 检测是否给定类的对象iscell 若是元胞数组则为真iscellstr 若是字符串组成的元胞数组则为真ischar 若是字符串则为真isequal 若两数组相同则为真isempty 若是空阵则为真isfinite 若全部元素都有限则为真isfield 若是构架域则为真isglobal 若是全局变量则为真ishandle 若是图形句柄则为真ishold 若当前图形处于保留状态则为真isieee 若计算机执行IEEE规则则为真isinf 若是无穷数据则为真isletter 若是英文字母则为真islogical 若是逻辑数组则为真ismember 检查是否属于指定集isnan 若是非数则为真isnumeric 若是数值数组则为真isobject 若是对象则为真isprime 若是质数则为真isreal 若是实数则为真isspace 若是空格则为真issparse 若是稀疏矩阵则为真isstruct 若是构架则为真isstudent 若是MATLAB学生版则为真iztrans 符号计算Z反变换A.2.10 A.2.10 I j , K kjacobian 符号计算中求Jacobian 矩阵jet 蓝头红尾饱和色jordan 符号计算中获得Jordan标准型keyboard 键盘获得控制权kron Kronecker乘法规则产生的数组A.2.11 A.2.11 L llaplace Laplace变换lasterr 显示最新出错信息lastwarn 显示最新警告信息leastsq 解非线性最小二乘问题(旧版)legend 图形图例length 数组长度light 创建光对象lighting 照明模式line 创建线对象lines 采用plot 画线色linmod 获连续系统的线性化模型linmod2 获连续系统的线性化精良模型linspace 线性等分向量ln 矩阵自然对数load 从MAT文件读取变量log 自然对数log10 常用对数log2 底为2的对数loglog 双对数刻度图形logm 矩阵对数logspace 对数分度向量lookfor 按关键字搜索M文件lower 转换为小写字母lsqnonlin 解非线性最小二乘问题lu LU分解A.2.12 A.2.12 M mmad 平均绝对值偏差magic 魔方阵maple 运作 Maple格式指令mat2str 把数值数组转换成输入形态串数组material 材料反射模式max 找向量中最大元素mbuild 产生EXE文件编译环境的预设置指令mcc 创建MEX或EXE文件的编译指令mean 求向量元素的平均值median 求中位数menuedit 启动设计用户菜单的交互式编辑工具mesh 网线图meshz 垂帘网线图meshgrid 产生“格点”矩阵methods 获知对指定类定义的所有方法函数mex 产生MEX文件编译环境的预设置指令mfun 对MAPLE经典函数实施数值计算mfunlis 能被mfun计算的MAPLE经典函数列表mhelp 引出 Maple的在线帮助min 找向量中最小元素mkdir 创建目录mkpp 逐段多项式数据的明晰化mod 模运算more 指令窗中内容的分页显示movie 放映影片动画moviein 影片帧画面的内存预置mtaylor 符号计算多变量Taylor级数展开A.2.13 A.2.13 N nndims 求数组维数NaN 非数(预定义)变量nargchk 输入宗量数验证nargin 函数输入宗量数nargout 函数输出宗量数ndgrid 产生高维格点矩阵newplot 准备新的缺省图、轴nextpow2 取最接近的较大2次幂nnz 矩阵的非零元素总数nonzeros 矩阵的非零元素.norm 矩阵或向量范数normcdf 正态分布累计概率密度函数normest 估计矩阵2范数norminv 正态分布逆累计概率密度函数normpdf 正态分布概率密度函数normrnd 正态随机数发生器notebook 启动MATLAB和Word的集成环境null 零空间num2str 把非整数数组转换为串numden 获取最小公分母和相应的分子表达式nzmax 指定存放非零元素所需内存A.2.14 A.2.14 O oode113 非Stiff 微分方程变步长解算器ode15s Stiff 微分方程变步长解算器ode23 非Stiff 微分方程变步长解算器ode23s Stiff 微分方程解算器ode23t 适度Stiff 微分方程解算器ode23tb Stiff 微分方程解算器ode45 非Stiff 微分方程变步长解算器odefile ODE 文件模板odeget 获知ODE 选项设置参数odephas2 ODE 输出函数的二维相平面图odephas3 ODE 输出函数的三维相空间图odeplot ODE 输出函数的时间轨迹图odeprint 在MATLAB指令窗显示结果odeset 创建或改写 ODE选项构架参数值ones 全1数组optimset 创建或改写优化泛函指令的选项构架参数值orient 设定图形的排放方式orth 值空间正交化A.2.15 A.2.15 P ppack 收集MATLAB内存碎块扩大内存pagedlg 调出图形排版对话框patch 创建块对象path 设置MATLAB搜索路径的指令pathtool 搜索路径管理器pause 暂停pcode 创建预解译P码文件pcolor 伪彩图peaks MATLAB提供的典型三维曲面permute 广义转置pi (预定义变量)圆周率pie 二维饼图pie3 三维饼图pink 粉红色图矩阵pinv 伪逆plot 平面线图plot3 三维线图plotmatrix 矩阵的散点图plotyy 双纵坐标图poisscdf 泊松分布概率分布函数poisspdf 泊松分布累计概率分布函数poissinv 泊松分布逆累计概率分布函数poissrnd 泊松分布随机数发生器pol2cart 极或柱坐标变为直角坐标polar 极坐标图poly 矩阵的特征多项式、根集对应的多项式poly2str 以习惯方式显示多项式poly2sym 双精度多项式系数转变为向量符号多项式polyder 多项式导数polyfit 数据的多项式拟合polyval 计算多项式的值polyvalm 计算矩阵多项式pow2 2的幂ppval 计算分段多项式pretty 以习惯方式显示符号表达式print 打印图形或SIMULINK模型printsys 以习惯方式显示有理分式prism 光谱色图矩阵procread 向MAPLE输送计算程序profile 函数文件性能评估器propedit 图形对象属性编辑器pwd 显示当前工作目录A.2.16 A.2.16 Q qquad 低阶法计算数值积分quad8 高阶法计算数值积分quit 推出MATLAB 环境quiver 二维方向箭头图quiver3 三维方向箭头图A.2.17 A.2.17 R rrand 产生均匀分布随机数randn 产生正态分布随机数randperm 随机置换向量range 样本极差rank 矩阵的秩rats 有理输出rcond 矩阵倒条件数估计real 复数的实部reallog 在实数域内计算自然对数realpow 在实数域内计算乘方realsqrt 在实数域内计算平方根realmax 最大正浮点数realmin 最小正浮点数rectangle 画“长方框”rem 求余数repmat 铺放模块数组reshape 改变数组维数、大小residue 部分分式展开return 返回ribbon 把二维曲线画成三维彩带图rmfield 删去构架的域roots 求多项式的根rose 频数扇形图rot90 矩阵旋转90度rotate 绕指定的原点和方向旋转rotate3d 启动三维图形视角的交互设置功能round 向最近整数圆整rref 简化矩阵为梯形形式rsf2csf 实数块对角阵转为复数特征值对角阵rsums Riemann和A.2.18 A.2.18 S ssave 把内存变量保存为文件scatter 散点图scatter3 三维散点图sec 正割sech 双曲正割semilogx X轴对数刻度坐标图semilogy Y轴对数刻度坐标图series 串联连接set 设置图形对象属性setfield 设置构架数组的域setstr 将ASCII码转换为字符的旧版指令shading 色彩浓谈模式shg 使当前图形窗位于前台shiftdim 数组维序号左移重组sign 根据符号取值函数signum 符号计算中的符号取值函数sim 运行SIMULINK模型simget 获取SIMULINK模型设置的仿真参数simple 寻找最短形式的符号解simplify 符号计算中进行简化操作simset 对SIMULINK模型的仿真参数进行设置simulink 启动SIMULINK模块库浏览器sin 正弦sinh 双曲正弦size 矩阵的大小slice 立体切片图solve 求代数方程的符号解spalloc为非零元素配置内存sparse 创建稀疏矩阵spconvert 把外部数据转换为稀疏矩阵spdiags 稀疏对角阵spfun 求非零元素的函数值sph2cart 球坐标变为直角坐标sphere 产生球面spinmap 色图彩色的周期变化spline 样条插值spones 用1置换非零元素sprandsym 稀疏随机对称阵sprank 结构秩spring 紫黄调春色图sprintf 把格式数据写成串spy 画稀疏结构图sqrt 平方根sqrtm 平方根矩阵squeeze 删去大小为1的“孤维”sscanf 按指定格式读串stairs 阶梯图std 标准差stem 二维杆图stem3 三维杆图step 阶跃响应指令str2double 串转换为双精度值str2mat 创建多行串数组str2num 串转换为数strcat 接成长串strcmp 串比较strjust 串对齐strmatch 搜索指定串strncmp 串中前若干字符比较strrep 串替换strtok 寻找第一间隔符前的内容struct 创建构架数组struct2cell 把构架转换为元胞数组strvcat 创建多行串数组sub2ind 多下标转换为单下标subexpr 通过子表达式重写符号对象subplot 创建子图subs 符号计算中的符号变量置换subspace 两子空间夹角sum 元素和summer 绿黄调夏色图superiorto 设定优先级surf 三维着色表面图surface 创建面对象surfc 带等位线的表面图surfl 带光照的三维表面图surfnorm 空间表面的法线svd 奇异值分解svds 求指定的若干奇异值switch-case-otherwise 多分支结构sym2poly 符号多项式转变为双精度多项式系数向量sym 创建一个符号变量symmmd 对称最小度排序symrcm 反向Cuthill-McKee排序syms 创建多个符号对象symsum 符号计算求级数和A.2.19 A.2.19 T ttan 正切tanh 双曲正切taylortool 进行Taylor逼近分析的交互界面text 文字注释tf 创建传递函数对象tic 启动计时器title 图名toc 关闭计时器trapz 梯形法数值积分treelayout 展开树、林treeplot 画树图tril 下三角阵trim 求系统平衡点trimesh 不规则格点网线图trisurf 不规则格点表面图triu 上三角阵try-catch 控制流中的Try-catch结构type 显示M文件A.2.20 A.2.20 U uuicontextmenu 创建现场菜单uicontrol 创建用户控件uimenu 创建用户菜单unmkpp 逐段多项式数据的反明晰化unwrap 自然态相角upper 转换为大写字母A.2.21 A.2.21 V vvar 方差varargin 变长度输入宗量varargout 变长度输出宗量vectorize 使串表达式或内联函数适于数组运算ver 版本信息的获取view 三维图形的视角控制voronoi Voronoi多边形vpa 任意精度(符号类)数值A.2.22 A.2.22 W wwarning 显示警告信息what 列出当前目录上的文件whatsnew 显示MATLAB中Readme 文件的内容which 确定函数、文件的位置while 控制流中的While环结构white 全白色图矩阵whitebg 指定轴的背景色who 列出内存中的变量名whos 列出内存中变量的详细信息winter 蓝绿调冬色图workspace 启动内存浏览器A.2.23 A.2.23 X x , Y y , Z zxlabel X轴名xor 或非逻辑yesinput 智能输入指令ylabel Y轴名zeros 全零数组zlabel Z轴名zoom 图形的变焦放大和缩小ztrans 符号计算Z变换A.3 A.3SIMULINK的库模块A.3.1 A.3.1库模块Demos library 演示子库Simulink SIMULINK基本库A.3.2 A.3.2连续模块子库Continuous Continuous 连续模块子库Derivative 求导数模块Integrator 连续函数积分Memory 记忆模块State-Space 状态方程模块Transfer Fcn 传递函数模块A.3.3 A.3.3离散模块子库DiscreteDiscrete 离散模块子库Discrete Filter 离散滤波器模块Discrete-Time Integrator 离散时间积分模块Discrete Transfer Fcn 离散传递函数模块Discrete Zero-Pole 离散零极点增益模块Unit Delay 单位延迟模块Zero-Order Hold 零阶保持模块A.3.4 A.3.4解析函数和查表函数模块子库Functions & TablesFcnC语言格式的任意函数模块Functions & Tables 解析函数和查表函数模块子库MATLAB Fcn MATLAB 语言格式的任意函数Look-Up Table 一维查表函数模块Look-Up Table(2-D) 二维查表函数模块A.3.5 A.3.5一般数学函数子库MathAbs 取绝对值模块Combinatorial Logic 组合逻辑模块Gain 增益模块Logical 逻辑运算模块MinMax 取极大值或极小值的模块Math 一般数学函数子库Mux 复用模块Product 乘法器Relational 关系运算模块Sign 符号取值模块Slider 滑键增益模块Sum 求和模块A.3.6 A.3.6非线性模块子库NonlinearDead Zone 死区非线性模块Nonlinear 非线性模块子库Relay 继电器非线性模块Saturation 饱和非线性模块A.3.7 A.3.7信号和系统模块子库Signal & SystemsDemux 分用模块Enable 使能模块Ground 接地模块In1 输入端口模块Merge 汇合模块Out1 输出端口模块Signal & Systems 信号和系统模块子库SubSystem 子系统模块Trigger 触发模块Terminator 终端模块A.3.8 A.3.8信宿模块子库SinksDisplay 数值显示模块Scope 示波模块Sinks 信宿模块子库Stop 终止仿真To File 把数据保存为文件To Workspace 把数据写成为矩阵变量XY Graph 显示X-Y图形A.3.9 A.3.9信源模块子库SourcesClock 连续仿真时钟模块Constant 恒值输出模块From File 从文件读数据From Workspace 从内存读数据Pulse 脉冲发生器Signal Generator 信号发生器Sine Wave 正弦波输出Sources 信源模块子库Step 阶跃输出A.4 A.4图形对象的属性A.4.1 A.4.1低层对象创建指令和句柄获取指令axesfigurefindobjgcagcbogcflightlinepatchrectanglesurfacetextuicontrol uicontextmenuuimenuA.4.2图形对象属性AcceleratorAmbient AmbientStrength BackgroundColor BackingStoreBoxCallback CameraPosition CameraPositionMode CameraTarget CameraTargetMode CameraUpVector CameraUpVectorMode CameraViewAngle CameraViewAngleMode CDataCDataMappingCheckedChildrenColorColormapColorOrder CurrentAxes CurrentFigure CurrentPointCurvature DataAspectRatio DefaultUicontrolUnits DefaultUicontrolFontname DefaultUicontrolFontsize DefaultUicontrolHorizontal DiffuseStrength EdgeColor EdgeLightingEnableEraseModeFaceColorFaceLightingFaces FaceVertexCData FileNameFontNameFontAngleFontSize ForgroundColor HandleVisibility Horizontal HorizontalAlignment LabelLineStyleLineWidthListboxTop MarkerMarkerEdgeColorMarkerFaceColorMarkerSizeMaxMenuBarMeshStyleMinNameNextPlotNumberTitlePaperPositionPaperUnitsParentPlotBoxAspectRatioModePointerPointerShapeCDataPointerShapeHotSpotPositionProjectionRotationSeparatorSpecularColorReflectanceSpecularExponentSpecularsSpecularStrengthStringStyleTagTitleToolBarTypeUnitsUserDataValueVerticalAlignmentVerticesVisibleWindow ButtonDownFcnWindow ButtonMotionFcnWindow ButtonUpFcnXColor, YColor, ZColorXData, YData, ZDataXGrid, YGrid, ZGridXLabel, YLabel, ZLabelXLim, YLim, ZLimXTick, YTick, ZTtickXTickMode, YTickMode, ZTickMode XTickLabel, YTickLabel, ZTickLabel XTickLabelMode, YTickLabelMode, ZTickLabelModeXScale, YScale, ZScaleXDir, YDir, ZDirXAxisLocation, YAxisLocation, ZAxisLocation。
Matlab基础及其应用教程复习2

习题讲解
教材 P102 3.8习题 2、3、5。
练习
1.多项式的四则运算 2.插值与拟合。 已知: x0=[0,1,2,3,4,5,6]; y0=[0,2.24,2.83,3.00,2.83,2.24,0] 求(1)在区间[0,6]上对已知点分别进行线性插值,三次方 程组插值。 (2)在区间[0,6]上对已知点分别进行二次拟合,三次拟合, 并求出相应的拟合多项式。
n( x) x 2 4x 7
a( x) x 3 2 x 2 3x 4
>> a=[1 2 3 4]; >> b=[1 4 9 16]; >> e = conv(a,b) e= 1 6 20 50 75 84 64
b( x) x 3 4x 2 9x 16
例3.7 求商多项式及余多项式。
例3.16:取余弦曲线上11个点的自变量和函数值点作为已知 数据,再选取41个自变量点,分别用分段线性插值、三次方 程式插值和样条插值三种方法计算确定插值函数的值。
>> x=0:10; y=cos(x); >> xi=0:.25:10; >> y0=cos(xi); %精确值 >> y1=interp1(x,y,xi); %线性插值结果 >> y2=interp1(x,y,xi,'cubic'); %三次方程式插值结果 >> y3=interp1(x,y,xi,'spline'); %样条插值结果 >> plot(xi,y0,'o',xi,y1,xi,y2,'-.',xi,y3)
第4章 MATLAB符号运算
MATLAB基础及应用教程

4.1.1 脚本文件(Scripts)
当我们需要在命令窗进行大量的命令集合运行时, 直接从命令窗口输入比较麻烦, 这时 就可以将这些命令集合存放在一个脚本文件(Scripts)中,运行时只需要输入其文件名就 可以自动执行这些命令集合。需要注意的是,脚本文件运行所产生的变量都驻留在 MATLAB 的工作空间中,同时脚本文件也可以调用工作空间中的数据。因此,脚本文件所涉及的变量 是全局变量。前几章所涉及到的 M 文件都是这类脚本文件。 编辑一个脚本文件可以直接在命令窗口的左上角打开编辑窗进行编辑。 4.1.2 函数文件(func构成 (1)函数定义行: Function [输出参量]=gauss(输入参量) (2): 完成函数的功能。 (3)函数说明。 (4)函数行注。 从上面构成的情况看, 函数式文件实际上是完成输入参量与输出参量的转换, 这样的转换是 由函数文件名为 gauss 的文件来完成的。 函数体的功能必须说明清楚输入参量与输出参量的 关系。函数说明是用来解释该函数的功能的,函数行注是对程序行进行说明的。上面(1) 和(2)是必须的。 【例 4-1】分析下面函数文件。 %一个数列,任意项等于前两项之和,输入项数可以给出这个数列 function [a]=sul(n) if n==1
n 的最大数为 100,要求: (1) 保存你的 fibo.m 文件,当在命令窗调用 fibo 函数时,不论输入任何整数有正确的 输出。 (2) 做出 fibo 的二维离散函数图,n 取 1 到 10,图的函数值处用小圆圈并涂为黑色,请 保存你的图形。 (3) 用三次样条插值的方法对(2)中的 10 个点进行插值,自变量的分辨率为 0.01, 请保存你的图形。 (4) 将完成(3)工作的插值函数保存为 fib.m 文件, ,当在命令窗调用 fib 函数时,不论输 入任何具有两位小数且小于 10 大于 0 的数(如 5.45)时有正确的输出。 7. 设电子粒子束流从恒定磁场中某点以相同速率发射, 发射的方向与磁场方向的夹角很小, 观察不同方向入射的粒子束流的运动轨道。 (设磁场沿 Z 方向) 数学模型: 粒子流的速度初值为
Matlab基础与应用教程答案(人民邮电出版社)

第1章一、思考题4. (1) B=A(2:5,1:2:5)(2)A(7)=[](3) A=A+30(4)size(A)(5) t(find(t==0))=eps(6)t=reshape(x,3,4)(7)abs('matlab')(8)char(93)5. A=[97,67,34,10;-78,75,65,5;32,5,-23,-59;0,-12,54,7](1) B=A(1:3,:)C=A(:,1:2)D=A(2:4,3:4)E=B*C(2)A(find(A>=50 & A<=100))二、实验题1.x=-74/180*pi;y=-27/180*pi;sin(x*x+y*y)/sqrt(tan(abs(x+y)))+pi2.a=-3.0:0.1:3.0;exp(-0.3*a).*sin(a+0.3)3.x=[2,4;-0.45,5];log(x+sqrt(1+x.*x))/24. A=[3,54,2;34,-45,7;87,90,15];B=[1,-2,67;2,8,74;9,3,0]; (1)A*Bans =129 432 41977 -407 -1052402 591 12489A.*Bans =3 -108 13468 -360 518783 270 0(2) A^3ans =-28917 240246 -4368137883 -259101 27669171333 252504 38673A.^3ans =27 157464 839304 -91125 343658503 729000 3375(3)A/Bans =-5.8845 5.3549 -0.20286.3554 -5.6596 4.3293-6.6325 6.2078 9.0241B/Aans =30.2855 19.2643 -8.561433.4394 21.1547 -9.3974-0.7443 -0.3938 0.2830(4)[A,B]ans =3 54 2 1 -2 6734 -45 7 2 8 7487 90 15 9 3 0[A([1,3],:);B^2]ans =3 54 287 90 15600 183 -81684 282 72615 6 8255.a=1+2i;b=3+4i;c=exp(pi*i/6)c =0.8660 + 0.5000ic+a*b/(a+b)ans =1.6353 + 1.8462i第2章一、思考题3.s=0;for n=0:63s=s+2^n;enddisp(s)n=0:63;s=sum(2.^n)二、实验题1.x=input('输入一个四位整数:');y=[fix(x/1000),mod(fix(x/100),10),mod(fix(x/10),10),mod(x,10)] z=mod((y+7),10)x=z(3)*1000+z(4)*100+z(1)*10+z(2)2.gh=input('输入工号');h=input('输入工时');dj=84;if h>120gz=dj*120+1.15*dj*(h-120);elseif h<60gz=dj*h-700;elsegz=dj*h;endformat bank;display([gh,gz])3.循环结构n=input('input n:');s=0;for k=1:ns=s+1/k^2;enddisplay(sqrt(s*6))向量运算n=input('input n:');k=1:n;display(sqrt(sum(1./k.^2)*6))4.y=0;k=0;while y<3k=k+1;y=y+1/(2*k-1);enddisplay([k-1,y-1/(2*k-1)])5.x0=0;x=1;k=0;a=input('a=');b=input('b=');while abs(x-x0)>=1e-5 && k<500x0=x;x=a/(b+x0);k=k+1;enddisplay([k,x]);display([(-b+sqrt(b^2+4*a))/2,(-b-sqrt(b^2+4*a))/2]);6.y=fun(40)/(fun(30)+fun(20))(1)函数文件fun.mfunction f=fun(n)f=n+log(n^2+5);(2)函数文件fun.mfunction f=fun(n)a=1:n;f=sum(a.*(a+1));第3章一、思考题4.t=0:0.01:1.5;y=sqrt(3)/2*exp(-4*t).*sin(4*sqrt(3)*t+pi/3);5.x=-10:0.01:10;y=linspace(-6,6,size(x,2))z=x.^3+3*x.*y.^2;plot3(x,y,z)6.x=100:100:400;y=100:100:400;z=[636,697,624,478;698,712,630,478;680,674,598,412;662,626,552,334]; [X,Y]=meshgrid(x,y);mesh(X,Y,z)二、实验题1.(1)x=-10:0.01:10;plot(x,x-x.^3/3/2)(2)plot(x,exp(-x.*x/2)/2/pi)(3)x=-8:0.01:8;plot(x,sqrt((64-x.*x)/2))(4)t=0:pi/100:6*pi;plot(t.*sin(t),t.*cos(t))2. (1)x1=linspace(0,1,100);y1=2*x1-0.5;t=linspace(0,pi,100);x=sin(3*t).*cos(t);y=sin(3*t).*sin(t);plot(x1,y1,'r-',x,y,'b:');text(0.4,1,'y=2x-0.5');text(-0.4,0.5,'x=sin(3t)cos(t)');text(-0.4,0.3,'y=sin(3t)sin(t)');(2)subplot(1,2,1);scatter(x1,y1,10)title('y=2x-0.5');subplot(1,2,2);scatter(x,y,10)3.subplot(1,2,1);x=1:1:100; y=sin(1./x);plot(x,y)subplot(1,2,2);fplot('sin(1/x)',[1,100])4.subplot(2,2,1);bar(t,y);subplot(2,2,2);stairs(t,y);subplot(2,2,3);stem(t,y)subplot(2,2,4);semilogy(t,y);5.theta=linspace(-pi,pi,100);ro=5.*cos(theta)+4;polar(theta,ro);fi=linspace(0,2*pi,100);a=1r=a.*(1+cos(fi));polar(fi,r);6.t=0:pi/20:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t);z=t; plot3(x,y,z);t=0:0.01:1;x=t;y=t.^2;z=t.^3;plot3(x,y,z);7.x=-30:0.1:0;y=0:0.1:30;[x,y]=meshgrid(x,y);z=10.*sin(sqrt(x.^2+y.^2))./sqrt(1+x.^2+y.^2); meshc(x,y,z);8.x=linspace(-3,3,100);y=linspace(-3,3,100);[x y]=meshgrid(x,y);fxy=-5./(1+x.^2+y.^2);i=find(abs(x)<=0.8 & abs(y)<=0.5);fxy(i)=NaN;surf(x,y,fxy)9.u=linspace(1,10,100);v=linspace(-pi,pi,100);[u,v]=meshgrid(u,v);x=3.*u.*sin(v);y=2.*u.*cos(v);z=4*u.^2;x=3*u.*sin(v);y=2*u.*cos(v);z=4*u.^2;surf(x,y,z);shading interp;light('position',[1,0,1]);10.t=0:pi/100:2*pi;y=sin(t);comet(t,y)第4章一、思考题5.(1)A=eye(3);(2)C=100+(200-100)*rand(5,6);(3)D=1+sqrt(0.2)*randn(1,500);(4)E=ones(size(A));(5)A=A+30*eye(size(A));(6)B=diag(diag(A))二、实验题1.P=pascal(5);H=hilb(5);Dp=det(P);Dh=det(H);Kp=cond(P);Kh=cond(H);P矩阵的性能更好,因为Kp较小2.A=[1,-1,2,3;0,9,3,3;7,-5,0,2;23,6,8,3]B=[3,pi/2,45;32,-76,sqrt(37);5,72,4.5e-4;exp(2),0,97]A1=diag(A);B1=diag(B);A2=triu(A);B2=triu(B);A3=tril(A);B3=tril(B);rA=rank(A);rB=rank(B);nA=norm(A);nb=norm(B);cA=cond(A);cB=cond(B);3.A=[31,1,0;-4,-1,0;4,-8,-2];[V,D]=eig(A);4.A=diag([-1,-1,-1,-1],-1)+diag([-1,-1,-1,-1],1)+diag([2,2,2,2,2]) b=[1,0,0,0,0]';x1=inv(A)*b;x2=A\b;[L,U]=lu(A);x3=U\(L\b);[Q,R]=qr(a);[Q,R]=qr(A);x4=R\(Q\b)R=chol(A);x5=R\(R'\b)5.B=sparse(A);x1=inv(B)*b;x2=B\b;[L,U]=lu(B);x3=U\(L\b);第5章一、思考题3.A=randn(10,5);mean(A)std(A)max(max(A))min(min(A))sum(A,2)sum(sum(A))sort(A,1)sort(A,2,'descend')二、实验题1.A=rand(1,30000);mean(A)std(A)max(A)min(A)size(find(A>0.5))/size(A)2.h=[466,715,950,1422,1635];w=[7.04,4.28,3.40,2.52,2.13];hh=[500,900,1500];ww=interp1(h,w,hh,'spline')3.x=linspace(1,10,50);y=log(x);f=polyfit(x,y,5);yy=polyval(f,x);plot(x,y,'r-',x,yy,'g.')4.N=64; % 采样点数T=5; % 采样时间终点t=linspace(0,T,N); % 给出N个采样时间ti(I=1:N) x=3*exp(-t); % 求各采样点样本值xdt=t(2)-t(1); % 采样周期f=1/dt; % 采样频率(Hz)X=fft(x); % 计算x的快速傅立叶变换XF=X(1:N/2+1); % F(k)=X(k)(k=1:N/2+1) f=f*(0:N/2)/N; % 使频率轴f从零开始plot(f,abs(F),'-*') % 绘制振幅-频率图xlabel('Frequency');ylabel('|F(k)|')5.(1)p1=[1 2 0 0 7];p2=[1 -2];p3=[1 0 0 5 1];p12=conv(p1,p2);p=p12+[zeros(1,size(p12,2)-size(p3,2)),p3];roots(p)(2)A=[-1,4,3;2,1,5;0,5,6];Px=polyval(p,A)Pxm=polyvalm(p,A)6.(1)z=fzero('3*x-sin(x)+1',0)(2)建立函数文件myfun.mfunction F=myfun(X)x=X(1);y=X(2);F(1)=x*x+y*y-9;F(2)=x+y-1;在命令窗口中输入以下命令:x=fsolve(@myfun,[3,0]',optimset('Display','off')) 第6章一、思考题2.fx=inline('1./(1+x.^2)');[I,n]=quad(fx,-100000,100000,1e-3);[I,n]=quadl(fx,-100000,100000,1e-3);x=-100000:0.01:100000;y=1./(1+x.*x);trapz(x,y);3.(1)fx=inline('-2*y+2*x*x+2*x');[t,y]=ode23(fx,[0,0.5],1)(2)fx=inline('y-exp(x)*cos(x)');[t,y]=ode23(fx,[0,3],1)二、实验题1.for x=1:3fx=[x,x^2,x^3;1,2*x,3*x;0,2,6*x];diff(fx)end2.(1)x=0:0.01:1;y=x.^10+10.^x+1./log10(x);dy=diff(y)/0.01;(2)x=0:0.01:1;y=log(1+x);dy=diff(y,2)/0.01;plot(x(1:99),dy)3.(1)fx=inline('x.^2.*sqrt(2*x.*x+3)');quad(fx,1,5)(2)fx=inline('x./sin(x).^2');quad(fx,pi/4,pi/3)(3)fx=inline('abs(cos(x+y))');dblquad(fx,0,pi,0,pi)(4)syms x y ;fx=x*y;int(int(fx,y^2,y+2),-1,2) x 的积分区间为【0,2】时 fx=inline('x.*y'); dblquad(fx,0,2,-1,2)4.x=0.3:0.2:1.5;y=[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325]; trapz(x,y)5. (1)yp=inline('-(1.2+sin(10*x))*y'); [t,y]=ode23(yp,[0,5],1);(2)令''3,,21y x y x y x ='==,则可写出原方程的状态方程形式:123233221sin 31)1(2cos 5cos x t x x t t t x x x x x +--++='='='⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+-+-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡'''t x x x t t tx x x cos 00)1(2cos 51sin 311000103212321建立函数文件ztfun.mfunction y=ztfun(t,x)b=[0;0;cos(t)];y=[0,1,0;0,0,1;-1/(3+sin(t)),-1,5*cos(2*t)/(t+1)^2]*x+b;解微分方程组[t,y]=ode23(@ztfun,[0,5],[1;0;2]);6.建立函数文件ztfun.m function yy=ztfun(t,y)yy=[y(2)*y(3);-y(1)*y(3);-0.51*y(1)*y(2)]; 解微分方程组[t,y]=ode23(@ztfun,[0,5],[0;1;1])第7章一、思考题 3.(1)数值积分fx=inline('exp(x).*(1+exp(x)).^2'); quad(fx,0,log(2)) 符号积分f=sym('exp(x)*(1+exp(x))^2'); v=int(f,0,log(2)); eval(v) (2)略 二、实验题 1.A=sym('[1,2,3;x,y,z;3,2,1]') rank(A) inv(A) det(A)2.(1) y=sym('sqrt(x+sqrt(x+sqrt(x)))');y1=diff(y)y2=diff(y,'x',2) (2) syms x y;fxy=sin(x^2*y)*exp(-x^2-y); diff(diff(fxy,x),y)3.(1) syms x int(1/(1+x^4))(2) syms x tint((-2*x*x+1)/(2*x*x-3*x+1)^2,x,cos(t),exp(2*t))4.syms n xsymsum(1/(2*n+1)/(2*x+1)^(2*n+1),n,0,inf) symsum(1/(2*n+1)/(2*x+1)^(2*n+1),n,0,5)5.(1) syms xtaylor((exp(x)+exp(-x))/2,5,0) (2) syms a xtaylor(exp(-5*x)*sin(3*x+pi/3),5,a)6.(1)x=solve(sym('x^3+a*x+1=0'))(2) [x y]=solve(sym('sqrt(x^2+y^2)-100=0,3*x+5*y-8=0'))7. 方程转化为: π21'11't y y y y -=+=⎩⎨⎧符号解[y1,y11]=dsolve('Dy=y1,Dy1+y=1-t^2/pi','y(-2)=5,y1(-2)=5','t') 数值解编写函数文件ztfun.mfunction yy=ztfun(t,y) yy=[y(2);1-t^2/pi-y(1)]; 在命令窗口输入以下命令[t,y]=ode45(@ztfun,[-2,7],[-5;5]); t=linspace(-2,7,49) y2=y8.[x,y]=dsolve('Dx=3*x+4*y,Dy=-4*x+3*y','x(0)=0,y(0)=1')第9章二、实验题 1.(1)新建一个Blank GUI 。
MATLAB基础讲课笔记内容

1 概述70年代中期,Cleve Moler和他的同事在美国国家自然科学基金的资助下开发了LINPACK和EISPACK 的Fortran子程序库,这两个程序库是解线性方程和特征值问题的,代表当时矩阵计算软件的最高水平。
70年代末期,Cleve Moler 在新墨西哥大学给学生开线性代数,想让学生使用LINPACK和EISPACK程序库,但他不希望学生在编程上花太多时间,毕竟课程的目的不是学习编程。
因此,他为学生编写了LINPACK和EISPACK的接口程序,这程序取名为MATLAB,即MATrix LABoratory。
又过了几年,Celve Moler 应邀到另一所大学讲学,他把MATLAB带到了这所学校。
在此之后的一两年内,MATLAB开始受到欢迎,并成为应用数学界的术语。
1983年春天,工程师John Little 觉察到MATLAB的潜在应用天地是工程领域,他与Moler、Steve Bangert一起开发了第二代专业版MATLAB,从这开始,MATLAB增加了数据图形功能。
1984年,MathWorks 公司成立,MATLAB正式推向市场。
在此之前,国际上已经有大量的数值计算软件包,它们大多使用Fortran、C编写,它们有共同的缺点:适应面窄、可扩充性差、不开放等等,一个软件包只能解决一个局部问题,很难以推广应用。
MATLAB 的出现打破了这一局面,它的高度适应性、可扩充性等优良特性吸引了很多科学界人士,很多人开始在MATLAB上开发自己的应用工具箱,MATLAB一步步壮大起来。
MATLAB具有良好的可扩展性,MATLAB的函数大多为ASCII文件,可以直接编辑、修改,MATLAB的工具箱可以任意增减,任何人可以自己生成MATLAB工具箱,因此很多研究成果被直接作成MATLAB工具箱发布。
比如在小波变换刚刚出现时,就出现了MATLAB的小波工具箱,到5.0版,小波工具箱成为MATLAB 的标准工具箱。
MATLAB基础及其应用教程

3.1多项式 3.2插值和拟合
3.3数值微积分
3.4线性方程组的数 值解
3.5稀疏矩阵
3.6常微分方程的数 值解
3.7小结 3.8习题
3.1.1多项式的表达和创建 3.1.2多项式的四则运算 3.1.3多项式求值和求根运算 3.1.4多项式的构造
3.2.1多项式插值和拟合 3.2.2最小二乘法拟合
6.1 M文件
6.2 MATLAB的程序 控制结构
6.3数据的输入与输 出
6.4 MATLAB文件操 作
6.5面向对象
1
编程
6.6 MATLAB
2
程序优化
3
6.7程序调试
4
6.8小结
5
6.9习题
6.1.1局部变量与全局变量 6.1.2 M文件的编辑与运行 6.1.3脚本文件 6.1.4函数文件 6.1.5函数调用
MATLAB基础及其应用教程
读书笔记模板
01 思维导图
03 目录分析 05 读书笔记
目录
02 内容摘要 04 作者介绍 06 精彩摘录
思维导图
关键字分析思维导图
内容
矩阵
基础知识
符号
运算
习题
数值
教程
教学
数值 小结
变量
基础
文件
符号
函数
程序
第章
存储
内容摘要
《MATLAB基础及其应用教程》基于MATLAB 7.1版,详细介绍了MATLAB的基础知识、数值计算、符号运算、 图形处理、程序设计、SIMULINK仿真等内容,为配合教学,各章编写了教学提示、教学要求和爿题,书后附有上 机实验指导。
0 6
1.10习题
0 5
MATLAB基础及其应用教程周开利邓春晖课后答案

第三章习题及参考答案解答:>> p=[1 -1 -1];>> roots(p)ans =-0.61801.6180解答:取n=5,m=61>> x=linspace(0,2*pi,5); y=sin(x);>> xi=linspace(0,2*pi,61);>> y0=sin(xi);>> y1=interp1(x,y,xi);>> y2=interp1(x,y,xi,'spline');>> plot(xi,y0,'o',xi,y1,xi,y2,'-.');>> subplot(2,1,1); plot(xi,y1-y0);grid on>> subplot(2,1,2); plot(xi,y2-y0);grid on分段线性和三次样条插值方法与精确值之差取n=11,m=61>> x=linspace(0,2*pi,11); y=sin(x);>> xi=linspace(0,2*pi,61);>> y0=sin(xi);>> y1=interp1(x,y,xi);>> y2=interp1(x,y,xi,'spline');>> plot(xi,y0,'o',xi,y1,xi,y2,'-.');>> subplot(2,1,1); plot(xi,y1-y0);grid on>> subplot(2,1,2); plot(xi,y2-y0);grid on分段线性和三次样条插值方法与精确值之差解答:>> x=[0,300,600,1000,1500,2000];>> y=[0.9689,0.9322,0.8969,0.8519,0.7989,0.7491]; >> xi=0:100:2000;>> y0=1.0332*exp(-(xi+500)/7756);>> y1=interp1(x,y,xi,'spline');>> p3=polyfit(x,y,3);>> y3=polyval(p3,xi);>> subplot(2,1,1);plot(xi,y0,'o',xi,y1,xi,y3,'-.');>> subplot(2,1,2);plot(xi,y1-y0,xi,y3-y0);grid on插值和拟合方法相比较,都合理,误差也相近。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
《matlab基础及其应用教程》笔记
前提:已经学过matlab很久了,但是将近半年没有再碰,打算接下去继续学习matlab,并在帮忙增加国内关于matlab资料较少的simmechanics物理建模仿真部分,因此继续开始学matlab,打算重新看基本基础的书回顾一下,并做简单笔记。
plus:括号内是自己加的,自己的经验以及想法等
第一章:基础
1、b='Hi,Miss Black';可以直接用于表示字符串
2、cd可以设置路径;和DOS系统下一样,但是不够直观,推荐直接使用GUI的路径设置
3、clf:清除图像内容;cla:清除坐标内容(GUI会用到);clc,clear
4、lookfor:寻找函数名称;doc:直接打开某个函数的帮助文档(doc quad:打开积分的帮助文档)
第二章:matlab语言基础
1、matlab主要数据类型.pdf(matlab的数据类型,十分重要)
2、&与&&区别:&左右两边都比较;&&只要左边不对,就不进行比较(表示右边有错误也不会管)
3、各种函数比较的最好加上(),以防止因为优先级的问题而被吞
4、线性等分:linspace;对数等分:logspace
5、矩阵标志:①,单下标;②,多下标;③,:表示整行或者整列;④,删除直接使用[];
6、cat:行/列拼接;reshape:矩阵变换;repmat:行列拼接;
7、一般的函数:matlab函数(自己总结).xls里面有
8、多维数组(三维,图片数组)
第三章:matlab运算
1、polyadd:相加;conv:相乘/卷积;deconv:相除
2、polyval:多项式带入求解;polyvalm:矩阵带入求解
3、roots:多项式求根;
4、poly2sym:多项式转为字符形式;poly:构造多项式
5、interpX:插值与拟合(和《数值计算方法》对应,X不同而采取不同的差值拟合方法)
6、polyfit:生成多项式拟合(偶尔两端会出现Roung现象,使整个设计出现偏差)
7、diff:微分;quad:积分
8、cumsum:积分矩形计算法;trapz:积分梯形计算法;quad:辛普森积分法;quad1:科茨积分法;
9、rref:解矩阵方程
10、jacobi:雅克比迭代法;gseidel:Gauss-Seidel 迭代法;
11、sparse:稀疏矩阵;full:全矩阵;spdiags:对角稀疏矩阵
12、spalloc:分配存储空间(对稀疏矩阵而言)
13、Euler:欧拉法解微分(1:前插;2:后插)
14、ode23:二三阶龙科库塔公式解偏微分;ode45:三四阶龙科库塔公式
第四章:结构数组与细胞数组(感觉暂时用到比较少,所以只是简单浏览了一下)
1、结构数组和C语言的一样,后面加.
2、或者使用struct构建结构数组
3、结构数组相关函数.pdf
4、rmfield:删除域;gerfield:获得结构数组中的值(GUI中十分有用)
5、setfield:重新设定值
6、cell:构建细胞数组
7、cell2struc/struc2cell:结构数组与细胞数组的互相转换
8、mat2cell/cell2mat:细胞数组与普通数组的转换
第五章:符号运算
1、sym():定义单个符号常量;syms:定义多个符号常量
2、sym:也可以定义矩阵符号
3、合并同类项(collect)、多项式展开(expand)、因式分解(factor)、一般化简(simplify)、不定化简(simple)、通分(numden)和书写格式美化(pretty)。
4、limit,diff,int,taylor均用于符号运算(上面的是各种数值计算方法)
5、fourier:傅里叶变换;ifourier:反傅里叶变换;拉普拉斯变换,Z变换等可以查询help 文档
6、solve:线性方程解析解;dsolve:常微分解析解
第六章:MATLAB 程序设计
1、global:定义全局变量
2、tic,toc:计算运行时间
3、各种和C语言一样的循环,比较语句等等
4、break;return;pause;continue区别与比较
5、input:键盘输入;disp:屏幕显示
6、fprintf:文本格式显示;fscanf:文本格式读取
7、fwrite:二进制读取;fread:二进制存储
8、fgetl:数据文字行读取;fgets:数据文字行存储
9、matlab文件操作函数.pdf
10,fopen:打开文件(GUI使用较为广泛);fclose:关闭
11、6.5:面向对象编程——自己理解(自己定义一个可以再次调用的函数)
第七章:可视化
1、figure:画图形(可以调用图形句柄,并且进行修改,axis等都可以)
2、其余的画法,参数设置可以看toolbox的~(貌似图形的很难得使用到)
3、plot绘图函数的常用参数.pdf
4、subplot:画子图
5、plot3:三维曲线;surf:三维曲面
6、图像分为:索引,二进制,灰度以及RGB三种
7、imshow:显示图片
8、imread:加载/读取图片
9、ezplot:函数绘图
第八章:simulink
因为这个不久前仿真过,所以基本不用看。