[matlab笔记]绘图时dataTip的设置
2009-04-29 16:43:45| 分类：Computer | 标签：|字号大中小订阅
用matlab(我用的版本是2009)绘图后，有时候需要显示某个数据据点的信息，这时候可以用工具栏上的Data Cursor工具。如下图所示
选中数据点之后，会弹出一个Data Tip，Data Tip显示的就是当前被选中的数据点的基本信息，比如图中显示的是坐标值。现在就遇到一个问题了，如何自己定制这个Data Tip呢？特别是，我的数据有效数字太多时，默认情况下只能显示4位，如何让它显示的精度更高呢？
matlab提供了解决方案，不是Option，也不是Preferences，而是脚本。我这里把matlab的m 文件称作脚本。习惯了图形界面的人可能会觉得麻烦，但是我却觉得这是最自由的解决方案。下面以一个实例来演示操作过程：
1、绘图
x=rand(1000,1); y=rand(1000,1); plot(x,y,'*');
这样，就在[0,1]X[0,1]这样一个矩形内绘制了1000*1000个点。绘图的时候要指定一个点的图例，不然默认情况下是不画点的，这样也就没法点选了。
2、选中Data Cursor工具，这样，鼠标放到图中间会变成一个十字，这时候点击一个数据点会弹出一个Data Tip框。
可以看到显示只有4位有效数字。
3、右键菜单，里面有两项"Edit Text Update Function"和"Select Text Update Function"。分别表示编辑配置文件，选择配置文件（配置文件指的还是m文件）。在未编辑之前，当前图用的是默认配置文件。编辑之后，需要保存到一个位置，新保存的文件只对当前图有效。下次如果还想使用这个配置文件，就要用到"Select Text Update Function"了。
4、点选"Edit Text Update Function"，出现一个m文件编辑器。内容如下：
function output_txt = myfunction(obj,event_obj)
% Display the position of the data cursor
% obj Currently not used (empty)
% event_obj Handle to event object
% output_txt Data cursor text string (string or cell array of strings).
pos = get(event_obj,'Position');
output_txt = {['X: ',num2str(pos(1),4)],...
['Y: ',num2str(pos(2),4)]};
% If there is a Z-coordinate in the position, display it as well
if length(pos) > 2
output_txt{end+1} = ['Z: ',num2str(pos(3),4)];
end
眼尖的人应该一眼就看到了几个"4"，没错，正是它们限制了显示精度，改成更大的数就行了。这几行代码的意思应该很明显了。第一行是函数原型，以%开头的是注释，pos是变量，output_txt是返回值。二维图分两行显示X,Y坐标，如果是三维，则还会显示Z坐标。
5、保存为m文件。注意保存的文件一定要带有.m后缀，不然matlab识别不了，从而Data Tip 会显示为Error in custom datatip string function。这个错误说函数返回错误，实际上就是不能识别。
6、保存完了，就算完事了。以后可以直接改那个已经保存的文件，然后绘图后，选择这个文件。再次提醒，每次绘图时使用的都是默认配置，想要特定的效果，要么重新写，要么选择指定配置文件。
7、关于这个配置文件本身，还有许多值得探讨的地方。我这个例子只能更改显示精度，实际上，它还有许多可以自由定制的地方。比如，可以显示点的序号，即第几个点。在配置文件末尾加上一句：
output_txt{end+1} = ['index:', num2str(event_obj.DataIndex)];
平面图形直观图的画法 先观察下面的图形，总结投影变化规律。 投影规律： 1.平行性不变;但形状、长度、夹 角会改变； 2.平行直线段或同一直线上的两条 线段的比不变 3.在太阳光下，平行于地面 的直线在地面上的投影长不变 表示空间图形的平面图形，叫做 空间图形的直观图 画空间图形的直观图，一般都要 遵守统一的规则， 1．斜二测画法 我们常用斜二测画法画空间图形及水平放置的平面多边形的直 观图．斜二测画法是一种特殊的平行投影画法． 2．平面图形直观图的画法 斜二测画法的步骤： (1)在已知图形中取互相垂直的x 轴和y 轴，两轴相交于点O .画直观 图时，把它们画成对应的x ′轴和y ′轴，两轴交于点O ′，且使 ∠x ′O ′y ′＝_45°(或135°)_，它们确定的平面表示_水平面． (2)已知图形中平行于x 轴或y 轴的线段，在直观图中分别画成_ 平行
于x′轴或y′轴的线段． (3)已知图形中平行于x轴的线段，在直观图中保持原长度不变_，_垂直于x轴的线段，长度为原来的_一半_． 注意点： 1．斜二测画法中的“斜”和“二测”分别指什么？ 提示：“斜”是指在已知图形的xOy平面内垂直于x轴的线段，在直观图中均与x′轴成45°或135°；“二测”是指两种度量形式，即在直观图中，平行于x′轴或z′轴的线段长度不变；平行于y′轴的线段长度变为原来的一半。 2．圆的斜二测画法，其图形还是圆吗？ 提示：不是圆，是一个压扁了的“圆”，即椭圆。 3．立体图形直观图的画法 由于立体图形与平面图形相比多了一个z轴，因此，用斜二测画法画立体图形的直观图时，图形中平行于x轴、y轴或z轴的线段在直观图中分别画成平行于x′轴、y′轴或z′轴的线段．平行于x轴和z轴的线段，在直观图中长度不变，平行于y轴的线段，长度为原来的一半． 例１．用斜二测画法画水平放置的六边形的直观图 解：
《数值计算方法》实验指导 （Matlab 版） 肇庆学院数学与统计学学院 计算方法课程组
1. 实验名称 实验1 算法设计原则验证(之相近数相减、大数吃小数和简化计算步骤) 2. 实验题目 有效数字的损失. 123 )与1000个较小的数(3 10 15)的和，验证 大数吃小数的现象. (3)分别用直接法和秦九韶算法计算多项式 P(x) a 0x n a 1x n 1 在x =1.00037 处的值?验证简化计算步骤能减少运算时间. n 1 对于第(3)题中的多项式P (x )，直接逐项计算需要n (n 1) 2 1 次乘法 和n 次加法，使用秦九韶算法 P(x) (((a °x ajx a 2)x a . 则只需要n 次乘法和n 次加法. 3. 实验目的 验证数值算法需遵循的若干规则. 4. 基础理论 设计数值算法时，应避免两个相近的数相减、防止大数吃小数、简化计算步骤减少运算 次数以减少运算时间并降低舍入误差的积累. 两相近的数相减会损失有效数字的个数， 用一 《数值计算方法》实验 1报告 班级： 20xx 级 XXXXx 班 学号： 20xx2409xxxx 姓名： XXX 成绩: ⑴取 z 1016，计算z 1 Z 和 1/（、z 1 Z)，验证两个相近的数相减会造成 (2)按不同顺序求一个较大的数( a n 1 X a n
个大数依次加小数，小数会被大数吃掉，乘法运算次数太多会增加运算时间. 5.实验环境 操作系统：Win dows xp ；程序设计语言：Matlab 6.实验过程 (1)直接计算并比较； (2)法1 :大数逐个加1000个小数，法2 :先把1000个小数相加再与大数加； (3)将由高次项到低次项的系数保存到数组A[n]中，其中n为多项式次数. 7.结果与分析 (1)计算的~1V Z = _______________________________ ,1/( ~1 < z) ____________________ . 分析： (2)123逐次加1000个3 10 6的和是_________________________ ，先将1000个3 10 6相 加，再用这个和与123相加得_______________________ . 分析： (3)计算__________ 次的多项式： 直接计算的结果是___________________ ，用时___________________ ; 用秦九韶算法计算的结果是____________________ ，用时 ___________________ 分析：
三维图形 一． 三维曲线 plot3(x1,y1,z1,选项1,x2,y2,z2,选项2,…,xn,yn,zn,选项n) 其中每一组x,y,z 组成一组曲线的坐标参数，选项的定义和plot 函数相同。当x,y ,z 是同维向量时，则x,y,z 对应元素构成一条三维曲线。当x,y ,z 是同维矩阵时，则以x,y,z 对应列元素绘制三维曲线，曲线条数等于矩阵列数。 Example1.绘制三维曲线。 程序如下： clf, t=0:pi/100:20*pi; x=sin(t); y=cos(t); z=t.*sin(t).*cos(t); %向量的乘除幂运算前面要加点 plot3(x,y,z); title('Line in 3-D Space'); xlabel('X');ylabel('Y');zlabel('Z'); grid on; 所的图形如下： -1 1 X Line in 3-D Space Y Z 二． 三维曲面 1. 产生三维数据 在MATLAB 中，利用meshgrid 函数产生平面区域内的网格坐标矩阵。
语句执行后，矩阵X 的每一行都是向量x ，行数等于向量y 的元素的个数，矩阵Y 的每一列都是向量y ，列数等于向量x 的元素的个数。 2. 绘制三维曲面的函数 surf 函数和mesh 函数 example2. 绘制三维曲面图z=sin(x+sin(y))-x/10。 程序如下： clf, [x,y]=meshgrid(0:0.25:4*pi); %产生平面坐标区域内的网格坐标矩阵 z=sin(x+sin(y))-x./10; surf(x,y,z); axis([0 4*pi 0 4*pi -2.5 1]); title('surf 函数所产生的曲面'); figure; mesh(x,y ,z); axis([0 4*pi 0 4*pi -2.5 1]); title('mesh 函数所产生的曲面'); -2.5 -2-1.5-1-0.500.51surf 函数所产生的曲面
Matlab图形绘制经典案例 1、 三维曲线 >> t=0:pi/50:10*pi; >> plot3(sin(2*t),cos(2*t),t) >> axis square >> grid on
2、一窗口多图形>> t=-2*pi:0.01:2*pi; >> subplot(3,2,1)
>> plot(t,sin(t)) >> subplot(3,2,2) >> plot(t,cos(t)) >> subplot(3,2,3) >> plot(t,tan(t)) >> axis([-pi pi -100 100]) >> subplot(3,2,4) >> plot(t,cot(t)) >> axis([-pi pi -100 100]) >> subplot(3,2,5) >> plot(t,atan(t)) >> subplot(3,2,6) >> plot(t,acot(t))
3、图形样式、标注、题字(也可以利用菜单直接Insert) >> x=0:pi/20:2*pi;
>> plot(x,sin(x),'b-.') >> hold on >> plot(x,cos(x),'r--') >> hold on >> plot(x,sin(x)-1,'g:') >> hold on >> plot(x,cos(x)-1) >> xlabel('x'); >> xlabel('x轴'); >> ylabel('y轴'); >> title('图形样式、标注等'); >> text(pi,sin(pi),'x=\pi'); >> legend('sin(x)','cos(x)','sin(x)-1','cos(x)-1'); >> [x1,y1]=ginput(1) %利用鼠标定位查找线上某点的值x1 = 2.0893 y1 = -0.5000 >> gtext('x=2.5') %鼠标定位放置所需的值在线上
怎样画示意图题 一、显示空间关系 1．火车长100m ，车头距离桥头200m ，桥长200m ，火车从静止开始以a =1m/s 2的加速度运动，求火车过桥经历的时间。 2．长5.0m 的铁链悬于O 点，O 点下方距离铁链下方15m 处有一个（偏离O 点正下方少许）钉子。求铁链无初速释放后经过钉子的时间是多长？（g 取10m/s ） 3．1999年高考题 在光滑水平面上有一质量m =1.0×10-3 Kg 、电量q =1.0×10-10 C 的带正电小球，静止在O 点。以O 点为原点，在该水平面内建立直角坐标系Oxy 。现突然加一沿x 轴正方向、场强大小E =2.0×106V/m 的匀强电场，使小球开始运动。经过1.0s ，所加电场突然变为沿y 轴正方向，场强大小仍为E =2.0×106 V/s 匀强电场。再经过1.0s ，所加电场又突然变为另一个匀强电场，使小球在此电场作用下经1.0s 速度变为零。求此电场的方向及速度变为零时小球的位置。 4．2006年理综Ⅰ卷第23题 天空有近似等高的浓云层。为了测量云层的高度，在水平地面上与观测者的距离为d =3.0km 处进行一次爆炸，观测者听到由空气直接传来的爆炸声和由云层反射来的爆炸声时间上相差Δt ＝6.0s 。试估算云层下表面的高度。已知空气中的声速v =1 3 km/s 。 5．2007年理综Ⅰ卷第23题 甲乙两运动员在训练交接棒的过程中发现：甲经短距离加速后能保持9m /s 的速度跑完全程：乙从起跑后到接棒前的运动是匀加速的。为了确定乙起跑的时机，需在接力区前适当的位置设置标记。在某次练习中，甲在接力区前013.5m s =处作了标记，并以9m /s v =的速度跑到此标记时向乙发出起跑口令。乙在接力区的前端听到口令时起跑，并恰好在速度达到与甲相同时被甲追上，完成交接棒，已知接力区的长度为L =20m. 求：(1)此次练习中乙在接棒前的加速度a 。 (2)在完成交接棒时乙离接力区末端的距离. 二、把立体关系转化为平面关系： 6．如图所示，abcd 是一竖直的矩形导线框，线框面积为S ，放在磁感应强度为B 的均匀水平磁场中．ab 边在水平面内且与磁场方向成60?角．则通过导线框的磁通量等于 （ ） (A)BS (B) 12BS (C) 2 2BS (D) 32 BS
数值分析编程作业
2012年12月 第二章 14.考虑梯形电阻电路的设计，电路如下： 电路中的各个电流{i1,i2,…,i8}须满足下列线性方程组： 12 123 234 345 456 567 678 78 22/ 2520 2520 2520 2520 2520 2520 250 i i V R i i i i i i i i i i i i i i i i i i i i -= -+-= -+-= -+-= -+-= -+-= -+-= -+= 这是一个三对角方程组。设V=220V，R=27Ω，运用追赶法，求各段电路的电流量。Matlab程序如下： function chase () %追赶法求梯形电路中各段的电流量 a=input('请输入下主对角线向量a='); b=input('请输入主对角线向量b='); c=input('请输入上主对角线向量c='); d=input('请输入右端向量d='); n=input('请输入系数矩阵维数n='); u(1)=b(1); for i=2:n l(i)=a(i)/u(i-1); u(i)=b(i)-c(i-1)*l(i); end y(1)=d(1); for i=2:n y(i)=d(i)-l(i)*y(i-1); end x(n)=y(n)/u(n); i=n-1; while i>0 x(i)=(y(i)-c(i)*x(i+1))/u(i); i=i-1; end x 输入如下:
请输入下主对角线向量a=[0,-2,-2,-2,-2,-2,-2,-2]; 请输入主对角线向量b=[2,5,5,5,5,5,5,5]; 请输入上主对角线向量c=[-2,-2,-2,-2,-2,-2,-2,0]; 请输入方程组右端向量d=[220/27,0,0,0,0,0,0,0]; 请输入系数矩阵阶数n=8 运行结果如下： x = 8.1478 4.0737 2.0365 1.0175 0.5073 0.2506 0.1194 0.0477 第三章 14.试分别用(1)Jacobi 迭代法；(2)Gauss-Seidel 迭代法解线性方程组 1234510123412191232721735143231211743511512x x x x x ?????? ??????---????????????=--?????? --?????? ??????---?????? 迭代初始向量 (0)(0,0,0,0,0)T x =。 （1）雅可比迭代法程序如下： function jacobi() %Jacobi 迭代法 a=input('请输入系数矩阵a='); b=input('请输入右端向量b='); x0=input('请输入初始向量x0='); n=input('请输入系数矩阵阶数n='); er=input('请输入允许误差er='); N=input('请输入最大迭代次数N='); for i=1:n for j=1:n if i==j d(i,j)=a(i,j); else d(i,j)=0; end end end m=eye(5)-d\a; %迭代矩阵 g=d\b; x=m*x0+g; k=1; while k<=N %进行迭代 for i=1:5 if max(abs(x(i)-x0(i))) >er x=m*x+g; k=k+1;
Matlab绘图 强大的绘图功能是Matlab的特点之一，Matlab提供了一系列的绘图函数，用户不需要过多的考虑绘图的细节，只需要给出一些基本参数就能得到所需图形，这类函数称为高层绘图函数。此外，Matlab 还提供了直接对图形句柄进行操作的低层绘图操作。这类操作将图形的每个图形元素（如坐标轴、曲线、文字等）看做一个独立的对象，系统给每个对象分配一个句柄，可以通过句柄对该图形元素进行操作，而不影响其他部分。 本章介绍绘制二维和三维图形的高层绘图函数以及其他图形控制函数的使用方法，在此基础上，再介绍可以操作和控制各种图形对象的低层绘图操作。 一、二维绘图 二维图形是将平面坐标上的数据点连接起来的平面图形。可以采用不同的坐标系，如直角坐标、对数坐标、极坐标等。二维图形的绘制是其他绘图操作的基础。 （一）绘制二维曲线的基本函数 在Matlab中，最基本而且应用最为广泛的绘图函数为plot，利用它可以在二维平面上绘制出不同的曲线。 1．plot函数的基本用法
plot函数用于绘制二维平面上的线性坐标曲线图，要提供一组x 坐标和对应的y坐标，可以绘制分别以x和y为横、纵坐标的二维曲线。plot函数的应用格式 plot(x,y) 其中x,y为长度相同的向量，存储x坐标和y坐标。 例51 在[0 , 2pi]区间，绘制曲线 程序如下：在命令窗口中输入以下命令 >> x=0:pi/100:2*pi; >> y=2*exp(-0.5*x).*sin(2*pi*x); >> plot(x,y) 程序执行后，打开一个图形窗口，在其中绘制出如下曲线 注意：指数函数和正弦函数之间要用点乘运算，因为二者是向量。 例52 绘制曲线 这是以参数形式给出的曲线方程，只要给定参数向量，再分别求出x,y向量即可输出曲线： >> t=-pi:pi/100:pi; >> x=t.*cos(3*t); >> y=t.*sin(t).*sin(t); >> plot(x,y) 程序执行后，打开一个图形窗口，在其中绘制出如下曲线 以上提到plot函数的自变量x,y为长度相同的向量，这是最常见、最基本的用法。实际应用中还有一些变化。
第一讲Matlab 基本数值计算 一、矩阵 在Matlab 中，一个矩阵可以使数学意义上的矩阵，也可以是标量或者向量。对于一个标量（一个数）可以将之作为1?1 的矩阵，而向量（一行或一列）则可以认为是1?n 或者n ?1的矩阵。另外，一个0 ?0 矩阵在Matlab 中被认为是空矩阵，用“[]”表示。 1、矩阵的创建 矩阵的创建可以有以下几种形式 ⑴直接输入 >> A=[1 2 3;4 3 7;2 4 1] 注意：每行间的元素用逗号或空格分开，行与行之间用分号或回车分开，矩阵标示是一对中括号[ ]。 也可以采用数组编辑器（Array Editor）像在Excel 电子表格中据那样输入数据。 ⑵通过语句和函数产生 常用的特殊矩阵：zeros:全零矩阵,ones：全1 矩阵，eye：单位矩阵，rand:随机矩阵，diag:对角阵等。 例：>> A=ones(3,4) >> E=eye(3) >> D=diag([3 5 2]) ⑶对矩阵进行裁剪或拼接 ⑷从外部文件装入数据
外部数据文件可以是以保存的Matlab 工作空间，也可以是文本（.txt）文件，或者是电子表格创建的文件（.xls）. 例：已知一个文本格式的数据文件E:\Mathmodel\data1.txt >> load e:\Mathmodel\data1.txt 得到一个变量名与文件名相同的矩阵（data1）。注意：文件的扩展名不能省略。 例：已知一个Excel 文件的路径为E:\Mathmodel\data2.xls a.缺省操作: >> NUMBER=xlsread('E:\Mathmodel\data2.xls') >>[NUMBER,TXT]=xlsread('E:\Mathmodel\data2.xls') 默认操作是从第一个工作表（sheet1）中提取数据。 b.从指定的工作表（而不是第一个）中提取数据： >> NUMBER=xlsread('E:\Mathmodel\data2.xls','S2') 或者 >> NUMBER=xlsread('E:\Mathmodel\data2.xls',2) c.从指定的工作表中读取指定区域的数据： >> NUMBER=xlsread('E:\Mathmodel\data2.xls',2,'g3:i8') 2、Matlab 的矩阵运算 ⑴基本运算 矩阵的加（+）、减（-）、乘（*）、乘方（^）运算法则与代数中的定义完全一致。例如： >> A=[1 2;3 4];B=[3 1;4 8];
第2章牛顿插值法实现 参考文献：[1]岑宝俊. 牛顿插值法在凸轮曲线修正设计中的应用[J]. 机械工程师,2009,10:54-55. 求牛顿插值多项式和差商的MA TLAB 主程序： function[A,C,L,wcgs,Cw]=newpoly(X,Y) n=length(X);A=zeros(n,n);A(:,1) =Y'; s=0.0;p=1.0;q=1.0;c1=1.0; for j=2:n for i=j:n A(i,j)=(A(i,j-1)-A(i-1,j-1))/(X(i)-X(i-j+1)); end b=poly(X(j-1));q1=conv(q,b);c1=c1*j;q=q1; end C=A(n,n);b=poly(X(n));q1=conv(q1,b); for k=(n-1):-1:1 C=conv(C,poly(X(k)));d=length(C);C(d)=C(d)+A(k,k); end L(k,:)=poly2sym(C);Q=poly2sym(q1); syms M wcgs=M*Q/c1;Cw=q1/c1; （1）保存名为newpoly.m 的M 文件 （2）输入MA TLAB 程序 >> X=[242,243,249,250]; >> Y=[13.681,13.526,13.098,13.095]; >> [A,C,L,wcgs,Cw]=newpoly(X,Y) 输出3阶牛顿插值多项式L 及其系数向量C 差商的矩阵A ，插值余项wcgs 及其 ) ()()1(ξ+n n f x R 的系数向量Cw 。 A = 13.6810 0 0 0 13.5260 -0.1550 0 0 13.0980 -0.0713 0.0120 0 13.0950 -0.0030 0.0098 -0.0003 C = 1.0e+003 *
第二讲 MatLab图形绘制功能 教学目的：熟练掌握二维曲线作图；了解隐函数绘图；掌握三维曲线和三维曲面作图的方法和技巧；了解MATLAB图形修饰处理方法；学会一些简单动画制作手段． 教学重点：二维曲线的作图, 极坐标和参数作图; 三维曲线和曲面作图. 教学难点：点运算在作图中的应用; 曲线和曲面作图的差别; 隐函数作图.通过解说，演示以及举例子克服难点难关． 一、二维平面图形 表1基本绘图函数
hold on 命令用于在已画好的图形上添加新的图形. plot是绘制一维曲线的基本函数，但在使用此函数之前，我们需先定义曲线上每一点的x及y 坐标.下例可画出一条正弦曲线： x=0:0.001:10; % 0到10的1000个点的x坐标 y=sin(x); % 对应的y坐标 plot(x,y); % 绘图 Y=sin(10*x);plot(x,y,'r:',x,Y,'b') % 同时画两个函数
若要改变颜色，在坐标对后面加上相关字串即可： x=0:0.01:10; plot(x,sin(x),'r') 若要同时改变颜色及图线型态（Line style），也是在坐标对后面加上相关字串即可：plot(x,sin(x),'r*') 用axis([xmin,xmax,ymin,ymax])函数来调整图轴的范围：
axis([0,6,-1.5,1]) MATLAB 也可对图形加上各种注解与处理： xlabel('x 轴'); % x 轴注解 ylabel('y 轴'); % y 轴注解 title('余弦函数'); % 图形标题 legend('y = cos(x)'); % 图形注解 gtext('y = cos(x)'); % 图形注解 ,用鼠标定位注解位置 grid on; % 显示格线 fplot 的指令可以用来自动的画一个已定义的函数分布图，而无须产生绘图所需要的一组数据作为变数.其语法为fplot('fun',[xmin xmax ymin ymax])，其中 fun 为一已定义的函数名称，例如 sin , cos 等等；而 xmin , xmax , ymin , ymax 则是设定绘图横轴及纵轴的下限及上限. 以下的例子是将一函数 f(x)=sin(x)/x 在-20
小学数学空间与图形复习资料（二） A、图形的认识 (一)线与角 一、线 1、直线：直线没有端点；长度无限，无法比较长短；过一点可以画无数条直线，过两点只能画一条直线。 2、射线：射线只有一个端点；长度无限，无法比较长短。 3、线段：线段有两个端点，它是直线的一部分；长度有限；两点的连线中线段最短。 4、平行线：在同一平面内，不相交的两条直线叫做平行线。两条平行线间的垂线段长度都相等。 5、垂线：两条直线相交成直角时，这两条直线互相垂直，其中一条直线叫做另一条直线的垂线，相交的点叫做垂足。 点到直线的距离：从直线外一点到这条直线所画的垂线段的长度叫做这点到直线的距离。 二、角 1、角的定义：从一点引出两条射线，所组成的图形叫做角。这个点叫做角的顶点，这两条射线叫做角的边。 2、角的特点：角的大小与角两边的长短无关，与角两边叉开的大小有关。 3、角的分类： 锐角：小于900的角叫做锐角；直角：等于900的角叫做直角；钝角：大于900而小于1800的角叫做钝角。平角：角的两边成一条直线，这时所组成的角叫做平角，平角1800。周角：角的一边旋转一周，与另一边重合，周角是3600。注意：平角不能理解为一条直线，周角不能理解为一条射线。 4、角的度量：量角器中心点与顶点重合，角的一边与量角器的零刻度线重合。即点与点重合，边与边重合的量角方法。看量角器的度数，就需要看刻度线在哪边了。 （二）平面图形 一、长方形特征：对边相等，4个角都是直角的四边形；有2条对称轴。 二、正方形特征：4条边都相等，4个角都是直角的四边形；有4条对称轴。 三、三角形 1、特征：由三条线段围成的图形；三角形两边之和大于第三条边；三角形内角和是180度；三角形具有稳定性；三角形有三条高。 2、分类： （1）按角分锐角三角形：三个角都是锐角。直角三角形：有一个角是直角；等腰直角三角形的两个锐角都为45度，它有1条对称轴。钝角三角形：有一个角是钝角。（2）按边分任意三角形：三条边长度不相等。等腰三角形：有两条边长度相等；两个底角相等；有1条对称轴。等边三角形：三条边长度都相等；三个内角都是60度；有3条对称轴。 四、平行四边形特征：两组对边分别平行，相对的边平行且相等； 五、梯形特征：只有一组对边平行的四边形；等腰梯形有1条对称轴。
第二讲MATLAB的数值分析 2-1矩阵运算与数组运算 矩阵运算和数组运算是MATLAB数值运算的两大类型，矩阵运算是按矩阵的运算规则进行的，而数组运算则是按数组元素逐一进行的。因此，在进行某些运算（如乘、除）时，矩阵运算和数组运算有着较大的差别。在MATLAB中，可以对矩阵进行数组运算，这时是把矩阵视为数组，运算按数组的运算规则。也可以对数组进行矩阵运算，这时是把数组视为矩阵，运算按矩阵的运算规则进行。 1、矩阵加减与数组加减 矩阵加减与数组加减运算效果一致，运算符也相同，可分为两种情况： （1）若参与运算的两矩阵（数组）的维数相同，则加减运算的结果是将两矩阵的对应元素进行加减，如 A=[1 1 1;2 2 2;3 3 3]; B=A; A+B ans= 2 2 2 4 4 4 6 6 6 （2）若参与运算的两矩阵之一为标量（1*1的矩阵），则加减运算的结果是将矩阵（数组）的每一元素与该标量逐一相加减，如 A=[1 1 1;2 2 2;3 3 3]; A+2 ans= 3 3 3 4 4 4 5 5 5 2、矩阵乘与数组乘 （1）矩阵乘 矩阵乘与数组乘有着较大差别，运算结果也完全不同。矩阵乘的运算符为“*”，运算是按矩阵的乘法规则进行，即参与乘运算的两矩阵的内维必须相同。设A、B为参与乘运算的 =A m×k B k×n。因此，参与运两矩阵，C为A和B的矩阵乘的结果，则它们必须满足关系C m ×n 算的两矩阵的顺序不能任意调换，因为A*B和B*A计算结果很可能是完全不一样的。如：A=[1 1 1;2 2 2;3 3 3]; B=A;
A*B ans= 6 6 6 12 12 12 18 18 18 F=ones(1,3); G=ones(3,1); F*G ans 3 G*F ans= 1 1 1 1 1 1 1 1 1 （2）数组乘 数组乘的运算符为“.*”，运算符中的点号不能遗漏，也不能随意加空格符。参加数组乘运算的两数组的大小必须相等（即同维数组）。数组乘的结果是将两同维数组（矩阵）的对应元素逐一相乘，因此，A.*B和B.*A的计算结果是完全相同的，如： A=[1 1 1 1 1;2 2 2 2 2;3 3 3 3 3]; B=A; A.*B ans= 1 1 1 1 1 4 4 4 4 4 9 9 9 9 9 B.*A ans= 1 1 1 1 1 4 4 4 4 4 9 9 9 9 9 由于矩阵运算和数组运算的差异，能进行数组乘运算的两矩阵，不一定能进行矩阵乘运算。如 A=ones(1,3); B=A; A.*B ans= 1 1 1 A*A ???Error using= =>
佛山科学技术学院 实 验 报 告 课程名称 数值分析 实验项目 数值积分 专业班级 机械工程 姓 名 余红杰 学 号 2111505010 指导教师 陈剑 成 绩 日 期 月 日 一、实验目的 1、 掌握程序的录入和matlab 的使用和操作； 2、 了解影响线性方程组解的精度的因素——方法与问题的性态。 3、 学会Matlab 提供的“\”的求解线性方程组。 二、实验要求 1、按照题目要求完成实验内容； 2、写出相应的Matlab 程序； 3、给出实验结果(可以用表格展示实验结果)； 4、分析和讨论实验结果并提出可能的优化实验。 5、写出实验报告。 三、实验步骤 1、用LU 分解及列主元高斯消去法解线性方程组 a)??????? ??=??????? ????????? ??----15900001.582012151526099999.2310 7104321x x x x ， 输出b Ax =中系数LU A =分解的矩阵L 和U ，解向量x 和)det(A ;用列主元法的行交换次序解向量x 和求)det(A ；比较两种方法所得结果。 2、用列主高斯消元法解线性方程组b Ax =。 （1）、???? ? ??=????? ??????? ??--11134.981.4987.023.116.427 .199.103.601.3321x x x
（2）、???? ? ??=????? ??????? ??--11134.981.4990.023.116.427 .199.103.600.3321x x x 分别输出)det(,,A b A ，解向量x ，（1）中A 的条件数。分析比较（1）、（2）的计算结果 3、线性方程组b Ax =的A 和b 分别为 ??????? ??=1095791068565778710A ,?????? ? ??=31332332b 则解T x ),1,1,1,1(=. 用MATLAB 内部函数求)det(A 和A 的所有特征值和2)(A cond . 若令 ?????? ? ??=+98.99599.6989.998.585604.508.72.71.8710A A δ, 求解b x x A A =++))((δδ，输出向量x δ和2x δ，从理论结果和实际计算两方面分析线性 方程组b Ax =解的相对误差22/x x δ以及A 的相对误差 /A A δ的关系。 四、实验结果 1： %run311.m clc,clear; A = [10 -7 0 1;-3 2.099999 6 2;5 -1 5 -1;2 1 0 2]; b = [8;5.90001;5;1]; %L U 分解 format short %小数点后四位，不然会受到后面的影响 [L U] = lu(A) %解方程组，输出A ，det(A) y = L\b; format long %小数点后15位显示 x = U\y
小学数学总复习空间与 图形试题 Document number【SA80SAB-SAA9SYT-SAATC-SA6UT-SA18】
空间与图形试题精选 一、填空题。 1. 从直线外一点到这条直线可以画无数条线段，其中最短的是和这条直线（）的线段。 2. 下图中，∠1=（）度，∠2=（）度。 1 30 2 3. 一个三角形中，最小的角是46°，按角分类，这个三角形是（）三角 形。 4. 右图是三个半径相等的圆组成的图形，它有（）条对称轴。 5. 用百分数表示以下阴影部分是整个图形面积的百分之几。 6. 把一个底面直径2分米的圆柱体截去一个高1分米的圆柱体，原来的圆柱体表面积减少（）平方分米。 7. “”和“”的周长之比是（），面积之比是（）。 8. 右图是由棱长1厘米的小正方体木块搭成的，这个几何体的表面积是 （）平方厘米。至少还需要（）块这样的小正方体才能搭成一个 大正方体。 9. 画一个周长厘米的圆，圆规两脚间的距离是（）厘米，画成的圆的面 积是（）。 10. 下面的小方格边长为1厘米，估一估图①中“福娃”的面积，算一算图②中阴影部分的面积。
11. 一个梯形，上底长a 厘米，下底长b 厘米，高h 厘米。它的面积是（ ）平方厘米。如果a=b ，那么这个图形就是一个（ ）形。 12. 在一块边长是20厘米的正方形木板上锯下一个最大的圆，这个圆的面积是（ ）平方厘米，剩下的边料是（ ）平方厘米。 13. 将一个大正方体切成大小相同的8个小正方体，每个小正方体的表面积是18平方厘米，原正方体的表面积是（ ）平方厘米。 14. 5个棱长为30厘米的正方体木箱堆放在墙角（如右图），露在外面的表面积是（ ）平方厘米。 15. 如下左图，已知大正方形的边长是a 厘米，小正方形的边长是b 厘米。用字母表示阴影部分的面积是（ ）平方厘米。 二、选择题。 1. 小青坐在教室的第3行第4列，用（4，3）表示，小明坐在教室的第1行第3列应当表示为（ ）。 A. （1，3） B. （3，1） C. （1，1） D. （3，3） 2. 在同一平面内，画已知直线的垂线，可以画（ ）。 A. 1条 B. 4条 C. 2条 D. 无数条 3. 用100倍的放大镜看40°的角，这个角的度数是（ ）度。 A. 4 B. 40 C. 400 D. 4000 4. 下面图形是用木条钉成的支架，最不容易变形的是（ ）。 D C B A 5. 下列图形中，对称轴条数最多的是（ ）。
列主元法 function lianzhuyuan(A,b) n=input('请输入n:') %选择阶数A=zeros(n,n); %系数矩阵A b=zeros(n,1); %矩阵b X=zeros(n,1); %解X for i=1:n for j=1:n A(i,j)=(1/(i+j-1)); %生成hilbert矩阵A end b(i,1)=sum(A(i,:)); %生成矩阵b end for i=1:n-1 j=i; top=max(abs(A(i:n,j))); %列主元 k=j; while abs(A(k,j))~=top %列主元所在行 k=k+1; end for z=1:n %交换主元所在行a1=A(i,z); A(i,z)=A(k,z); A(k,z)=a1; end a2=b(i,1); b(i,1)=b(k,1); b(k,1)=a2; for s=i+1:n %消去算法开始m=A(s,j)/A(i,j); %化简为上三角矩阵 A(s,j)=0; for p=i+1:n A(s,p)=A(s,p)-m*A(i,p); end b(s,1)=b(s,1)-m*b(i,1); end end X(n,1)=b(n,1)/A(n,n); %回代开始 for i=n-1:-1:1 s=0; %初始化s for j=i+1:n s=s+A(i,j)*X(j,1);
end X(i,1)=(b(i,1)-s)/A(i,i); end X 欧拉法 clc clear % 欧拉法 p=10; %贝塔的取值 T=10; %t取值的上限 y1=1; %y1的初值 r1=1; %y2的初值 %输入步长h的值 h=input('欧拉法please input number(h=1 0.5 0.25 0.125 0.0625):h=') ; if h>1 or h<0 break end S1=0:T/h; S2=0:T/h; S3=0:T/h; S4=0:T/h; i=1; % 迭代过程 for t=0:h:T Y=(exp(-t)); R=(1/(p-1))*exp(-t)+((p-2)/(p-1))*exp(-p*t); y=y1+h*(-y1); y1=y; r=r1+h*(y1-p*r1); r1=r; S1(i)=Y; S2(i)=R; S3(i)=y; S4(i)=r; i=i+1; end t=[0:h:T]; % 红线为解析解，'x'为数值解 plot(t,S1,'r',t,S3,'x')
空间曲面绘图 （1） 直接绘图-ezmesh 和ezsurf ezmesh 绘制三维网格图，ezsurf 绘制三维表面图。 例1 绘制抛物柱面2x 2z ?=的图形。 指令：ezmesh('2-x^2',[-1,1,-1,1]) 图7.43 指令：ezsurf('2-x^2',[-1,1,-1,1]) 图7.44 例2 绘制)xy sin(z =的图形。 指令：ezmesh('sin(x*y)',[0,4,0,4])
图7.45 例3 绘制马鞍面2y 2x z 2 2?=的图形。 指令：ezmesh('x^2/2-y^2/2') 图7.46 例4 绘制椭圆抛物面22x y 2z +=和抛物柱面2x 2z ?=所围的图形。 ezmesh('2-x^2',[-1,1,-1,1] hold on % 在同一图形窗口中继续绘图 ezmesh('2*y^2+x^2',[-1,1,-1,1]) axis([-1,1,-1,1,0,4])
图7.47 （2） 自定义网格绘图-meshgrid 和mesh/surf 例5 绘制23y x z ?=的图形。 [x,y]=meshgrid(-2:0.2:2,-3:0.1:3) % 自定义网格数据 z=x.^3-y.^2 mesh(x,y,z) 图7.48 surf(x,y,z) 图7.49 contour3(x,y,z,50) %绘制50条三维等高线
图7.50 Contour(x,y,z,40) %绘制40条二维等高线 图7.51 例6 在同一坐标系中绘制23y x z ?=和0z =的图形。 [x,y]=meshgrid(-2:0.2:2,-3:0.1:3) z=x.^3-y.^2 mesh(x,y,z) zz=zeros(size(z)) hold on mesh(x,y,zz)
如何控制matlab中的数据输出格式？ 2009/04/02 22:33 [useful ] 如何控制matlab中的数据输出格式？ 针对数据的显示形式控制，matlab中有专用命令format 但该命令不影响数据的储存形式和计算精度 format 默认格式 format short 5字长定点数 format long 15字长定点数 format short e 5字长浮点数 format long e 15字长浮点数 format short g format long g format hex 16进制 format bank 定点货币形式 format rat 小数分数表示 format + ＋，－，空格 format compact 压缩空格 format loose 包括空格和空行 当然也可以不用format 命令，可以修改系统的默认设置格式，File->Preferences->Command Window->Text Display. 如果想得到分数可以用rats()函数 如果想得到根号，只有用符号表示了，当然用符号表示是万能的。例如： a=1/3 得到0.33333 rats(a) 得到1/3 a='sqrt(3)' 得到sqrt(3) 不过这是字符型，要转换成数值型可以用eval(a),或str2num(a). 但是无论Matlab中采取什么样的输出格式，在系统内核中的变量的精度总是保持精确的（尽可能）。除非你人为的改变它的计算精度。我们可以用一个简单的例子来说明这个问题：如a=1/3 显示为0.33333333 a=a*3 得到1 （不是0.99999999） 在很多其它的程序设计语言中是不可能得到1的。这就说明了Matlab在计算的过程中不会损失用户的计算信息（包括中间结果）。t ormat：设置输出格式 对浮点性变量，缺省为format short.
MATLAB 编程题库 1.下面的数据表近似地满足函数2 1cx b ax y ++=，请适当变换成为线性最小二乘问题，编程求最好的系数c b a ,,，并在同一个图上画出所有数据和函数图像. 625 .0718.0801.0823.0802.0687.0606.0356.0995 .0628.0544.0008.0213.0362.0586.0931.0i i y x ---- 解： x=[-0.931 -0.586 -0.362 -0.213 0.008 0.544 0.628 0.995]'; y=[0.356 0.606 0.687 0.802 0.823 0.801 0.718 0.625]'; A=[x ones(8,1) -x.^2.*y]; z=A\y; a=z(1); b=z(2); c=z(3); xh=-1:0.1:1; yh=(a.*xh+b)./(1+c.*xh.^2); plot(x,y,'r+',xh,yh,'b*')
2.若在Matlab工作目录下已经有如下两个函数文件，写一个割线法程序，求出这两个函数 10 的近似根，并写出调用方式： 精度为10 解： >> edit gexianfa.m function [x iter]=gexianfa(f,x0,x1,tol) iter=0; while(norm(x1-x0)>tol) iter=iter+1; x=x1-feval(f,x1).*(x1-x0)./(feval(f,x1)-feval(f,x0)); x0=x1;x1=x; end >> edit f.m function v=f(x) v=x.*log(x)-1; >> edit g.m function z=g(y) z=y.^5+y-1; >> [x1 iter1]=gexianfa('f',1,3,1e-10) x1 = 1.7632 iter1 = 6 >> [x2 iter2]=gexianfa('g',0,1,1e-10) x2 = 0.7549 iter2 = 8