最新Matlab第五章-数据和函数的可视化资料ppt课件
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
❖ >> y=[9.5600 8.1845 5.2616 2.7917 2.2611 1.7340 1.2370 1.0674 1.1171 0.7620];
❖ >> q=polyfit(x,y,2);w=polyfit(x,y,4); ❖ >> e=0:0.0001:3; ❖ >> r=polyval(q,e);r1=polyval(w,e); ❖ >> plot(e,r,'b',e,r1,'r',x,y,'hb'); ❖ legend('2阶','4阶','原始数据')
❖ >> roots(p) ❖ ans = ❖ -2.4008 ❖ 0.3254 + 0.4544i ❖ 0.3254 - 0.4544i ❖ 0.1429 + 0.3499i ❖ 0.1429 - 0.3499i ❖ >> polyval(p,ans(1,1)) ❖ ans = ❖ 7.2742e-013 ❖ >> polyval(p,-2.4008) ❖ ans = ❖ -0.0198
❖ >> plot(x,y,'r')
❖ >> hold on
❖ >> plot(x,zeros(size(x)))
❖ >> hold off
❖ >> zoom on
%获局部放大图
❖ >>[tt,yy]=ginput(5);zoom off
❖ >> tt
❖ tt =
x 106 4
❖ -0.9908
❖ -0.4378
❖ end
❖ >> xx
❖ xx =
❖ -2.4008 -2.4008 -2.4008 -2.4008
题目4：已知数据如下，请用2阶及4阶多项式对数据进行拟 合，并比较拟合效果。
❖ >> x=[0.0129 0.0247 0.0530 0.1550 0.3010 0.4710 0.8020 1.2700 1.4300 2.4600];
❖ >> y_c=inline('(4*x.^3+7*x.^2-5*x+3).*(7*x.^2-2*x+1)')
❖ y_c =
❖ Inline function:
❖ y_c(x) = (4*x.^3+7*x.^2-5*x+3).*(7*x.^2-2*x+1)
❖ >> x=-10:0.01:10;
❖ >> y=y_c(x);
上网查询极限的定义，并用matlab计算
❖ >> clear ❖ >> syms x ❖ >> limit(((1+tan(x))^0.5-(1+sin(x))^0.5)/(x*(1-cos(x))),x,0) ❖ ❖ ans = ❖ ❖ 1/2
第五章 数据和函数的可视化
离散函数和连续函数的可视化 二维绘图指令plot
Matlab第五章-数据和函数 的可视化资料
上节课作业
❖ 题目1：上网查询矩阵行列式的定义，随机生成一个3阶方 阵，用定义及matlab自带函数分别计算其行列式的值。
❖ A=rand(3) ❖ A= ❖ 0.9501 0.4860 0.4565 ❖ 0.2311 0.8913 0.0185 ❖ 0.6068 0.7621 0.8214 ❖ >> a=(A(1,1)*A(2,2)*A(3,3)+A(1,2)*A(2,3)*A(3,1)+
A(1,3)*A(2,1)*A(3,2))- (A(1,3)*A(2,2)*A(3,1)+A(1,1)*A(2,3)*A(3,2)+ A(1,2)*A(2,1)*A(3,3)) ❖ a= ❖ 0.4289 ❖ >> b=det(A) ❖ b= ❖ 0.4289
题目2：用matlab求解下列方程组
❖ >> A=[-1 2 0;2 1 1;4 5 7;1 1 5]; ❖ b=[2;3;0;-5]; ❖ >> x=A\b ❖ x= ❖ 1.1818 ❖ 1.3766 ❖ -1.5974 ❖ >> A*x-b ❖ ans = ❖ -0.4286 ❖ -0.8571 ❖ 0.4286 ❖ -0.4286
坐标控制和图形标识axis 多次叠绘和多子图hold，subplot
获取数据指令ginput
离散函数和连续函数的可视化
离散数据和离散函数的可视化
众所周知：一对实数标量 (x, y)可表示为平面上的一个点；进而， 一对实数“向量”(x, y) 可表现为平面上的一组点。MATLAB就是利用这 种几何比拟法实现了离散数据可视化。
❖ >> A=[-1 2 0;2 1 1;4 5 7];b=[ห้องสมุดไป่ตู้;3;0]; ❖ >> x=A\b
❖ x=
❖ 1.7273
❖ 1.8636
❖ -2.3182
❖ >> A*x-b
❖ ans =
❖ 1.0e-015 *
❖ 0.4441
❖
0
❖
0
题目3：已知f(x)=4x3+7x2-5x+3，g(x)=7x2-2x+1，用roots（上 网查询或使用matlab帮助）命令及图像法求出f(x)*g(x)=0的 解。
❖ >> f=[4 7 -5 3] ❖ f= ❖ 4 7 -5 3 ❖ >> g=[7 -2 1] ❖ g= ❖ 7 -2 1 ❖ >> p=conv(f,g) ❖ p= ❖ 28 41 -45 38 -11 3 ❖ >> help roots ❖ ROOTS Find polynomial roots. ❖ ROOTS(C) computes the roots of the polynomial whose coefficients ❖ are the elements of the vector C. If C has N+1 components, ❖ the polynomial is C(1)*X^N + ... + C(N)*X + C(N+1). ❖ See also POLY, RESIDUE, FZERO. ❖ Overloaded methods ❖ help gf/roots.m ❖ help localpoly/roots.m
3
❖ 0.5300 2
❖ 1.9124
❖ -2.7880
1
❖ >> yy
❖ yy =
0
❖ 1.0e+004 * -1
❖ 5.9942
❖ -4.2398
-2
❖ -4.2398
❖ -4.2398
-3
-10 -8 -6 -4 -2
0
2
4
6
8 10
❖ 5.9942
❖ >> for n=1:4
❖ [xx(n),yyv(n)]=fzero(y_c,tt(n));
❖ >> q=polyfit(x,y,2);w=polyfit(x,y,4); ❖ >> e=0:0.0001:3; ❖ >> r=polyval(q,e);r1=polyval(w,e); ❖ >> plot(e,r,'b',e,r1,'r',x,y,'hb'); ❖ legend('2阶','4阶','原始数据')
❖ >> roots(p) ❖ ans = ❖ -2.4008 ❖ 0.3254 + 0.4544i ❖ 0.3254 - 0.4544i ❖ 0.1429 + 0.3499i ❖ 0.1429 - 0.3499i ❖ >> polyval(p,ans(1,1)) ❖ ans = ❖ 7.2742e-013 ❖ >> polyval(p,-2.4008) ❖ ans = ❖ -0.0198
❖ >> plot(x,y,'r')
❖ >> hold on
❖ >> plot(x,zeros(size(x)))
❖ >> hold off
❖ >> zoom on
%获局部放大图
❖ >>[tt,yy]=ginput(5);zoom off
❖ >> tt
❖ tt =
x 106 4
❖ -0.9908
❖ -0.4378
❖ end
❖ >> xx
❖ xx =
❖ -2.4008 -2.4008 -2.4008 -2.4008
题目4：已知数据如下，请用2阶及4阶多项式对数据进行拟 合，并比较拟合效果。
❖ >> x=[0.0129 0.0247 0.0530 0.1550 0.3010 0.4710 0.8020 1.2700 1.4300 2.4600];
❖ >> y_c=inline('(4*x.^3+7*x.^2-5*x+3).*(7*x.^2-2*x+1)')
❖ y_c =
❖ Inline function:
❖ y_c(x) = (4*x.^3+7*x.^2-5*x+3).*(7*x.^2-2*x+1)
❖ >> x=-10:0.01:10;
❖ >> y=y_c(x);
上网查询极限的定义，并用matlab计算
❖ >> clear ❖ >> syms x ❖ >> limit(((1+tan(x))^0.5-(1+sin(x))^0.5)/(x*(1-cos(x))),x,0) ❖ ❖ ans = ❖ ❖ 1/2
第五章 数据和函数的可视化
离散函数和连续函数的可视化 二维绘图指令plot
Matlab第五章-数据和函数 的可视化资料
上节课作业
❖ 题目1：上网查询矩阵行列式的定义，随机生成一个3阶方 阵，用定义及matlab自带函数分别计算其行列式的值。
❖ A=rand(3) ❖ A= ❖ 0.9501 0.4860 0.4565 ❖ 0.2311 0.8913 0.0185 ❖ 0.6068 0.7621 0.8214 ❖ >> a=(A(1,1)*A(2,2)*A(3,3)+A(1,2)*A(2,3)*A(3,1)+
A(1,3)*A(2,1)*A(3,2))- (A(1,3)*A(2,2)*A(3,1)+A(1,1)*A(2,3)*A(3,2)+ A(1,2)*A(2,1)*A(3,3)) ❖ a= ❖ 0.4289 ❖ >> b=det(A) ❖ b= ❖ 0.4289
题目2：用matlab求解下列方程组
❖ >> A=[-1 2 0;2 1 1;4 5 7;1 1 5]; ❖ b=[2;3;0;-5]; ❖ >> x=A\b ❖ x= ❖ 1.1818 ❖ 1.3766 ❖ -1.5974 ❖ >> A*x-b ❖ ans = ❖ -0.4286 ❖ -0.8571 ❖ 0.4286 ❖ -0.4286
坐标控制和图形标识axis 多次叠绘和多子图hold，subplot
获取数据指令ginput
离散函数和连续函数的可视化
离散数据和离散函数的可视化
众所周知：一对实数标量 (x, y)可表示为平面上的一个点；进而， 一对实数“向量”(x, y) 可表现为平面上的一组点。MATLAB就是利用这 种几何比拟法实现了离散数据可视化。
❖ >> A=[-1 2 0;2 1 1;4 5 7];b=[ห้องสมุดไป่ตู้;3;0]; ❖ >> x=A\b
❖ x=
❖ 1.7273
❖ 1.8636
❖ -2.3182
❖ >> A*x-b
❖ ans =
❖ 1.0e-015 *
❖ 0.4441
❖
0
❖
0
题目3：已知f(x)=4x3+7x2-5x+3，g(x)=7x2-2x+1，用roots（上 网查询或使用matlab帮助）命令及图像法求出f(x)*g(x)=0的 解。
❖ >> f=[4 7 -5 3] ❖ f= ❖ 4 7 -5 3 ❖ >> g=[7 -2 1] ❖ g= ❖ 7 -2 1 ❖ >> p=conv(f,g) ❖ p= ❖ 28 41 -45 38 -11 3 ❖ >> help roots ❖ ROOTS Find polynomial roots. ❖ ROOTS(C) computes the roots of the polynomial whose coefficients ❖ are the elements of the vector C. If C has N+1 components, ❖ the polynomial is C(1)*X^N + ... + C(N)*X + C(N+1). ❖ See also POLY, RESIDUE, FZERO. ❖ Overloaded methods ❖ help gf/roots.m ❖ help localpoly/roots.m
3
❖ 0.5300 2
❖ 1.9124
❖ -2.7880
1
❖ >> yy
❖ yy =
0
❖ 1.0e+004 * -1
❖ 5.9942
❖ -4.2398
-2
❖ -4.2398
❖ -4.2398
-3
-10 -8 -6 -4 -2
0
2
4
6
8 10
❖ 5.9942
❖ >> for n=1:4
❖ [xx(n),yyv(n)]=fzero(y_c,tt(n));