《大数据可视化部分问题详解》由会员分享,可在线阅读,更多相关《大数据可视化部分问题详解(7页珍藏版)》请在金锄头文库上搜索。
1、数据可视化 (即画图 )部分作业x acost1 已知椭圆的长、 短轴 a 4,b 2,用“小红点线” 画椭圆。(参见图 p5-1)y bsint程序:t=2*pi*(0:)/;a=4;b=2;x=a*cos(t);y=b*sin(t);plot(x,y,r.,MarkerSize,15)axis equal图 p5-12 根据表达式1 cos 绘制如图 p5-2 的心脏线。 (提示: 采用极坐标绘线指令 polar) clft = 0:.01:2*pi;P=1-cos(t); pline=polar(t,P,r), set(pline,LineWidth,5)title(P=1-costhe
2、ta)图 p5-23 A,B,C 三个城市上半年每个月的国民生产总值见表p5.1。试画出如图 p5-3 所示的三城市上半年每月生产总值的累计直方图。(提示:使用指令 bar)表 p5.1 各城市生产总值数据(单位:亿元)城市1月2月3月4月5月6月A170120180200190220B120100110180170180C70508010095120clear x=1:6;Y=170 120180200190220;120100110180170180;70508010095120bar(x,Y,grouped);bar(x,Y,stacked);colormap(cool);% legen
3、d(A,B,C,Location,NorthWest) legend(A,B,C,2)图 p5-34 二 阶 线 性 系 统 的 归 一化(即令n 1 ) 冲 激 响 应 可 表 示 为 :1 e t sin( t)01y(t)te t1 ,其中 |12 | , 为阻尼1e ( )t e ()t12系数。试在同一图上,画出0.2 :0.2:1.4 不同取值时的各条曲线,时间区间为t 0,18 。此外,要求( 1)1的各条曲线用蓝色,1的用红色,1 的用黑要求图形参见图 p5-4 )% exmp504.m clc,clf,clear; t=(0:0.05:18); N=length(t);zet
4、a=0.2:0.2:1.4;L=length(zeta); y=zeros(N,L); hold on for k=1:Lzk=zeta(k); beta=sqrt(abs(1-zk2);if zk-1-2*epsy=1/beta*exp(-zk*t).*sin(beta*t); plot(t,y,b) if zk0.4text(2.2,0.63,zeta = 0.2) endelseif (zk-1)1.2text(0.3,0.14,zeta = 1.4) endendend text(10,0.7,Deltazeta=0.2) axis(0,18,-0.4,0.8) hold off bo
5、x on grid on图 p5-45用绿实线绘制 x sin(t) ,y cos(t),z t的三维曲线, 曲线如图 p5-5所示。(提示:使用 plot3 指令)clf t=4*pi*(0:100)/100; x=sin(t);y=cos(t);z=t; plot3(x,y,z,g,LineWidth,3),box on图 p5-56 在区域 x,y 3,3 ,绘制 z24xe x的如图 p5-6 的三维(透视)网格曲面。不得使用 ezmesh)clfx=-3:0.1:3;y=x;X,Y=meshgrid(x,y);mesh(X,Y,Z) hidden off% colormap(cool
6、),% shading interp,syms x y zz=4*x*exp(-x2-y2);ezmesh(z,-3,3)hidden off图 p5-67在x,y 4 ,4 区间里,根据表达式 z sin(x y) ,绘制如图 p5-7所示的曲面。 xyclear allx=4*pi*(-50:50)/50;y=x;X,Y=meshgrid(x,y);Z=sin(X+Y)./(X+Y+(X+Y=0)*eps);surf(X,Y,Z)view(21,32)shading interp图 p5-78 试用图解法回答:1)方程组(1 x2 y2) 0.1 有多少个实数解?sin(x cos( y)
7、 02)求出离 x 0,y0最近、且满足该方程组的一个近似解。1) ezplot(y/(1+x2+y2)-0.1,-2*pi,2*pi,-pi/3,3.5*pi) hold on ezplot(sin(x+cos(y),-2*pi,2*pi,-pi/3,3.5*pi)图 p5-8可看到 6个交点,即方程组有 6 个实数解 (2)edit sushu function F=myfun8(x,y) F=y/(1+x2+y2)-0.1;sin(x+cos(y) end zoom onxy= ginput(1) f=fsolve(myfun8,xy) X,Y = ginput(1) f=fsolve(myfun8,X,Y) clear allsyms x ys=solve(y/(1+x2+y2)-0.1,sin(x+cos(y)离 x 0, y 0 最近的一个近似解为-0.9801y00 =0.2005