浙江大学远程教育学院《可视化计算》课程作业2０16年（夏)姓名:贺鹏学号:７1612８012０33年级:16秋专本计算机学习中心：浙大校内直属学习中心(华家池）————————————————————————————注意:所有图像的标题必须呈现足够你本人信息1.(10分）求解下列线性方程组的解：2.(10分）信号y=５*sin(pｉ*２０*t)＋３*cos（2＊pｉ＊5０*t)幅度为1的白噪声的干扰，请画出此信号,然后进行傅立叶变换，画出变换后的频域信号。ｃleaｒ;ｔitle(＇贺鹏71612８012０３3')t=0:0.０01:０.6;ｙ=5＊sin（pi＊20*ｔ)+３＊ｃoｓ（2＊pｉ＊5０*t);y=y+randn（1，lｅｎgth（t)）;subpｌot(１,２，1)plot（t,ｙ)tｉtｌｅ('贺鹏'）xlabel（'变换前信号＇)Y＝ｆft（y，５12);sｕbｐｌot(1,２,2）Ｙ=rｅａl(Y);pｌot(Y(1：5１2）)tｉtｌｅ(＇7１61２8０12033＇）xｌabel(＇变换后的频域信号')3.(10分）在空间有一个带正电的点电荷,请画出此点电荷的空间电位分布与电场强度的空间分布图。clear；ｓuｂpｌot(１，2,１)［x，y］＝meshgｒｉd(-４：0．２:4,－4:０.2:4）;ｚ=１./ｓqrt(ｘ.＾２+y.^2+0．0１);meｓh(x,y,z）titlｅ('贺鹏')subｐloｔ(１,2,２)[x,y]=meｓhｇrｉd(-２:０.１:２，-2:0.1：2）；z＝1./sqrt((x）．^２+y.＾2＋０.01);［ｄｘ,dｙ]=gｒaｄｉent(z）;conｔｏuｒ(ｘ,y,z，[-12,－8,-5,-３,-1,-0.５,-０.1,0.1,０.5,１,3，5,８，１２])；holｄonquiｖer(x,ｙ，ｄx，dy)tiｔle（'71612801203３')4.(10分)仿照课本第１1章的太阳｜地球｜月亮｜卫星,绕转动画实例,呈现地球绕太阳运转的动画。Ｃleａr；ｆigure('name','地球绕太阳转动'）;tiｔle('贺鹏716１２80１20３３'）ｓ１=[０:.0１:2*ｐi]；holdon;aｘiｓequal;ａxisｏffr1=10；w１＝1;t=０;pａusｅｔime=.002;siｔａ1=0;sita2＝0；ｓet(gcｆ，'doublebuffｅr','ｏｎ')pｌot(-20,１8，＇colｏr',＇r＇,＇ｍａrｋer'，'.','maｒｋerｓiｚｅ',4０)；ｔext（－１7，1８,＇太阳');p1=pｌot（-20，16,'ｃｏlor'，＇b','markｅr＇,'.'，'mａrkerｓize',20）；tｅxt(-17,１６,＇地球');ｐlot(０,0,＇coloｒ＇,'r','maｒｋer','.',＇maｒｋｅrsize',60）；ploｔ(r１*cｏｓ(ｓ１)，ｒ1*siｎ(s1））;sｅt(gcａ,'ｘliｍ',[-2020]，'ylim'，[-2０2０］）;p1=plｏｔ（r１*cｏs（sita1),r1*sｉn（sita１)，'color'，＇b'，'markｅr','.＇,＇ｍarkeｒsｉｚe＇，30);while１ｓet（p１,'xｄａtａ',r１*cos(sｉta１）,'yｄaｔa'，ｒ1＊ｓiｎ(sｉｔa1）);ｓiｔa1=ｓiｔa1＋w1*pａuｓetiｍe;ｐause(paｕsetimｅ);drawｎowＥnd5.(10分)设计一个低通滤波器,从混合信号:x(ｔ)=ｓｉn(2＊pｉ*10*t）+cｏs（2*pi*100*t)+0．2*ｒanｄｎ（size(ｔ))中获取１0Hz的信号（１0分）。clear；fiｇuｒe('name','贺鹏７1６1２80１203３＇）;wｓ＝１000;ｔ=0：1／ws:0.4;x=３*sin(2*ｐi*10＊t)＋２*ｃoｓ(2*pｉ*100*t)+０.8＊randn(sｉｚe(t)）;wn＝wｓ/2;[B,A］=buｔｔer（１0，30／wn);y＝ｆｉｌtｅr(B,A,x);plｏｔ（t,x,＇b-')holｄonｐlot(t,ｙ,'r','MａrkｅｒSiｚe',10)ｌegend（'iｎput',＇Outpuｔ'）6.(20分)设计一个程序,应用函数subploｔ(1,2,1）、subplot(１,２,2)分别显示您本人的二张照片,然后对二张照片分别进行傅立叶变换,并分别画出变换后的频域信号。再把２个频域信号相加,经傅立叶逆变换后，显示时域信号的图像。clear;A=iｍｒead('C:\1．jpｇ')；sｕbplｏt(1,２,1)imshow(A);tiｔle（'贺鹏71612801２０33')B=iｍreaｄ('Ｃ：\2.ｊpg');ｓubplｏt（1，2,2)imshow(B);tｉtle(＇贺鹏71612８0１2０３3')A1＝fft2（A);B1=ｆft２(B）;suｂplｏｔ(1,2,1);imshoｗ(A1）titｌｅ（'贺鹏71612８0１203３')ｓubｐｌｏｔ(1,2,2);ｉmsｈow(B1）tｉtlｅ('贺鹏７1612801２033')频域信号相加后经傅立叶逆变换C1=A１+B１;C=ifft2(C1)；imshｏw(C）xlａbel('71612８0１20３3＇);ylabel(＇贺鹏');7.(30分)小论文根据工作中的实际需要,请设计一个实际工程问题的可视化。可以选择以下之一:（1）工程动画的可视化；(２)大数据处理中的可视化;（3）算法与模型计算的可...
