可视化计算2013年夏作业.docVIP

浙江大学远程教育学院 2013年夏《可视化科学计算》上机操作 姓名： 学 号： 年级： 学习中心： ————————————————————————————— 调试下列程序： h0=figure(toolbar,none,position,[450 280 370 230],... name,图形对象属性的获取和设置实例); t=0:pi/40:4*pi; hline1=plot(t,sin(t)); xlabel(t); ylabel(sin(t)); grid; ht=title(正弦曲线); hline2=line(t+0.2,sin(t),Color,r); haxes=gca; h0_color=get(h0,Color) haxes_color=get(haxes,Color) hline1_color=get(hline1,Color) hline2_color=get(hline2,Color) set(haxes,GridLineStyle,--); set(hline2,LineWidth,4,Color,[0.8 0.8 0.8]); set(ht,FontAngle,italic,FontWeight,bold); set(gca,Children,[hline1 hline2]); 写出程序的运行结果； h0_color=0.8000 0.8000 0.8000 haxes_color=1 1 1 hline1_color=0 0 1 hline2_color=1 0 0 修改程序，程序的标题为你的姓名，把正弦线由你学号顺序构成的各个位数的折线。（本题10分） h0=figure(toolbar,none,position,‘[7 1 3 0 0 9 0 1 2 0 0 1]’, 张德军,图形对象属性的获取和设置实例); t=0:pi:4*pi; hline1=plot(t,sin(t)); xlabel(t); ylabel(sin(t)); grid; ht=hline(713009012001); haxes=gca; 在同一画面中应用函数subplot(m, n, p)按上、下分别画出y=sin(2x),y=sin(x)*cos(2x)的图形，图形的标题是你单位的全称，横坐标中显示的是你的姓名，纵坐标是你的学号。（本题10分）。　 %zhdj_2.m% subplot(2,1,1) ezplot(sin(2*x))title(兴达数码快印)xlabel(张德军);ylabel(713009012001)subplot(2,1,2)ezplot(sin(x)*cos(2*x))title(兴达数码快印)xlabel(张德军);ylabel(713009012001) 3、有一质点，初速度为v0,加速度为a，经过t时间后，质点所运动的位移。请写出matlab源代码。 （本题10分）。 %zhdj_3.m% Clc;clear all;close all ; v0=input(请输入初始速度：); a=input(请输入加速度：); t=0:0.1:5; s=v0*t+0.5*a*t.^2; plot(t,s) xlabel(t/s) ylabel(s/m) 4、有下列实验数据： X 1 2 3 4 5 6 7 8 9 Y 10.2 14.1 13.9 12.3 10.5 11 12.4 13.5 16 用三次多项式拟合，求出表达式的值，画出原数值图与拟合示意图，图的标题设置为你的学号，并测算x为7.5、12时的y值（本题10分）。 %zhdj_4.m% Clc;clear all;close all; x=1:9;y=[10.2 14.1 13.9 12.3 10.5 11 12.4 13.5 16];p=polyfit(x,y,3);plot(x,y,r)y1=polyval(p,x);hold onplot(x,y1,y)y2=polyval(p,[7.5 12]) %当x为7.5,12时的y值 Clc;clear all;close all; a1=imread(d:\zhdj.jpg); imshow(a1) a2=rgb2gray(a1); a3=edge(a2); imshow(a3) ８、本题按以下步骤进行： （1）对信号y=sin(314t)+sin(628t)加上白噪声，画出图像 （2）进行离散傅立叶变换，并画出傅立叶变换后的频谱图、功率谱图