1、浙江大學(xué)遠(yuǎn)程教育學(xué)院可視化計(jì)算課程作業(yè)2015年（冬）姓名：學(xué) 號(hào)：年級(jí)：學(xué)習(xí)中心：注意：所有圖像的標(biāo)題必須呈現(xiàn)足夠你本人信息1. （10分）求解下列線性方程組的解：解：A=1 0 7 3 3;1.5 3.6 10 7 8;2 0.5 -3 8.5 0;9 -4 22 21 90; 7 -4 33 6 -20B=3;-4;20;5;16X=B'/AA = 1.0000 0 7.0000 3.0000 3.0000 1.5000 3.6000 10.0000 7.0000 8.0000 2.0000 0.5000 -3.0000 8.5000 0 9.0000 -4.0000 22.00

2、00 21.0000 90.0000 7.0000 -4.0000 33.0000 6.0000 -20.0000B = 3 -4 20 5 16X = 3.5056 -0.8979 -0.2745 0.1438 0.01372. （10分）信號(hào)y = 5*sin(pi*20*t)+3*cos(2*pi*50*t)幅度為1的白噪聲的干擾，請(qǐng)畫(huà)出此信號(hào)，然后進(jìn)行傅立葉變換，畫(huà)出變換后的頻域信號(hào)。 cleart=(0:0.001:1)'y = 5*sin(2*pi*10*t)+3*cos(2*pi*50*t);plot(y(1:150)title('');050100150

3、-8-6-4-202468?傅立葉轉(zhuǎn)變后： clear t=0:0.001:0.6; y=5*sin(2*pi*10*t)+3*cos(2*pi*50*t); plot(y(1:150);y=fft(y,512);Y=real(y);figure,plot(Y(1:512)title(' ') 0100200300400500600-200-1000100200300400500?3. （10分）在空間有一個(gè)帶正電的點(diǎn)電荷，請(qǐng)畫(huà)出此點(diǎn)電荷的空間電位分布與電場(chǎng)強(qiáng)度的空間分布圖。hh=figure('numbertitle','off','n

4、ame','點(diǎn)電荷的空間電位分布與電場(chǎng)強(qiáng)度的空間分布圖<2015年秋 >');r0=0.12;th=linspace(0,2*pi,13);x,y=pol2cart(th,r0)x=x;0.1*xy=y;0.1*yplot(x,y)grid onhold onplot(0,0,'o','MarkerSize',12)k=9e9;q=1e-9;r0=0.12;u0=k*q/r0;u=linspace(1,3,7)*u0x=linspace(-r0,r0,100)X,Y=meshgrid(x)r=sqrt(X.2+Y.2);U=k

5、*q./r;contour(X,Y,U,u)xlabel('x','fontsize',16)ylabel('y','fontsize',16)title(' ')-0.2-0.15-0.1-0.0500.050.10.15-0.2-0.15-0.1-0.0500.050.10.15xy?4. （10分）仿照課本第11章的太陽(yáng)地球月亮衛(wèi)星，繞轉(zhuǎn)動(dòng)畫(huà)實(shí)例，呈現(xiàn)地球繞太陽(yáng)運(yùn)轉(zhuǎn)的動(dòng)畫(huà)。clear;clc;close all;%相對(duì)應(yīng)的圓點(diǎn)坐標(biāo) 半徑等x0=0;y0=0;r0=80;Lmin0=0;Lmax0=0;T0=2

6、160;w0=0*pi/T0;q0=0;x1=0;y1=0;r1=40;Lmin1=25;Lmax1=30;T1=1080;w1=pi/T1;q1=0;%初始化hh=figure('numbertitle','off','name','地球 | 太陽(yáng)轉(zhuǎn)演示動(dòng)畫(huà) 余合哲');%設(shè)置擦除方式sun = line(0,0,'color','r','linestyle','.','erasemode','xor','markersize

7、',r0);%太陽(yáng)earth = line(x0,y0,'color','k','linestyle','.','erasemode','xor','markersize',r1);%地球%添加標(biāo)注axis offtitle('太陽(yáng) | 地球 ','fontname','宋體','fontsize',9,'FontWeight','demi','color',

8、9;black');text(-50,50,'太陽(yáng)');%對(duì)太陽(yáng)進(jìn)行標(biāo)識(shí)line(-55,50,'color','r','marker','.','markersize',80);text(-50,40,'地球');%對(duì)地球進(jìn)行標(biāo)識(shí)line(-55,40,'color','k','marker','.','markersize',40);%繪制軌道s1=0:.01:2*pi;line(Lmax1*

9、cos(s1),Lmin1*sin(s1),'linestyle',':');%地球軌跡axis(-60,60,-60,60);%開(kāi)始畫(huà)圖t=0;while 1 ifishandle(hh) return; end q0=t*w0;q1=t*w1;t=t+1;%設(shè)置運(yùn)動(dòng)規(guī)律 if t>=4320; t=0; end%得到一周 就重置 x0=Lmax0*cos(q1); y0 = Lmin0*sin(q1);%設(shè)置太陽(yáng) x1 = x0+Lmax1*cos(q1); y1=y0+Lmin1*sin(q1);%設(shè)置地球坐標(biāo) set(sun,'xdata&

10、#39;,x0,'ydata',y0); set(earth,'xdata',x1,'ydata',y1); drawnow;end太陽(yáng) | 地球 ?5. （10分）設(shè)計(jì)一個(gè)低通濾波器，從混合信號(hào)：x(t)=sin(2*pi*10*t)+cos(2*pi*100*t)+0.2*randn(size(t)中獲取10Hz的信號(hào)（10分）。clearhh=figure('numbertitle','off','name','低通濾波器<2015年秋 >');randn('

11、;state',1)ws=1000;%采樣頻率t=0:1/ws:0.4x=sin(2*pi*10*t)+cos(2*pi*100*t)+0.2*randn(size(t);%生成帶噪聲的多頻信號(hào)wn=ws/2;%Nyquest 頻率B,A=butter(10,30/wn);y=filter(B,A,x);plot(t,x,'b-',t,y,'r.','MarkerSize',10);legend('Input','Output',0);title(' ')00.050.10.150.20.2

12、50.30.350.4-2.5-2-1.5-1-0.500.511.522.5?InputOutput6. （20分）設(shè)計(jì)一個(gè)程序，應(yīng)用函數(shù)subplot(1,2,1)、subplot(1,2,2)分別顯示您本人的二張照片，然后對(duì)二張照片分別進(jìn)行傅立葉變換，并分別畫(huà)出變換后的頻域信號(hào)。再把2個(gè)頻域信號(hào)相加，經(jīng)傅立葉逆變換后，顯示時(shí)域信號(hào)的圖像。原圖：clearfigure(1);I=imread('d:a.jpg');K= imread('d:b.jpg');subplot(1,2,1),imshow(I)title('zhangqi')xlab

13、el('照片a');subplot(1,2,2),imshow(K)title('zhangqi')xlabel('照片b');figure(2)X1=fft(I)subplot(2,2,1),plot(X1(1:10)title('zhangqi')xlabel('照片a變換后');X2=fft(K);subplot(2,2,2),plot(X2(1:512)title('zhangqi')xlabel('照片b變換后');X3=X1+X2;X4=ifft(X3)subplot(2

14、,2,3),plot(X4(1:512)title('wangji')xlabel('相加后逆變換');7. （30分）小論文 根據(jù)工作中的實(shí)際需要，請(qǐng)?jiān)O(shè)計(jì)一個(gè)實(shí)際工程問(wèn)題的可視化?？梢赃x擇以下之一：（1）工程動(dòng)畫(huà)的可視化；（2）大數(shù)據(jù)處理中的可視化；（3）算法與模型計(jì)算的可視化；（4）實(shí)際生產(chǎn)流程的可視化；（5）或其它有創(chuàng)新意義的可視化科學(xué)計(jì)算。要求：（1）題目有實(shí)際意義。（2）有分析、算法描述（3）程序源代碼設(shè)計(jì)。（4）問(wèn)題結(jié)果有可視化顯示。（4）題目的問(wèn)題有一定的新意。小論文的字?jǐn)?shù)不能少于2500字，格式由下列各部分組成：中文題目 Matlab在公司中如何

15、用來(lái)計(jì)算工資摘要： 隨著社會(huì)的不斷發(fā)展，許多大型的企業(yè)也在迅速崛起，企業(yè)員工也從以前的幾十個(gè)人發(fā)展到成百上千，乃至數(shù)萬(wàn)人。可想而知，單純的靠人力去計(jì)算每個(gè)人都每月的實(shí)發(fā)工資肯定是不現(xiàn)實(shí)的。最近學(xué)習(xí)了Matlab以后，發(fā)現(xiàn)它不僅僅能編寫(xiě)一些簡(jiǎn)單的程序，更可以把許多數(shù)據(jù)以圖像的方式呈獻(xiàn)給人們，這對(duì)于企業(yè)的管理以及工資的統(tǒng)計(jì)工作有著巨大的幫助。因?yàn)閷W(xué)習(xí)時(shí)間也不是很久，所以我在根據(jù)已學(xué)的知識(shí)基礎(chǔ)上，參考了一些其他書(shū)籍，集合書(shū)本上他人的設(shè)計(jì)經(jīng)驗(yàn)，自己寫(xiě)了一個(gè)小程序，可以很快的計(jì)算出一個(gè)人每月的實(shí)際工資。中文關(guān)鍵詞： 圖像 統(tǒng)計(jì) 工資英文題目：Matlab in the company how to ca

16、lculate salary英文摘要： With the continuous development of society, also is in rapid rise, many large enterprises from the dozens of personal development to enterprise employees also have had hundreds of thousands, even tens of thousands of people. Predictably, the pure by manpower to calculate each mon

17、thly net salary is not realistic. Recent learning Matlab, found that it not only can write some simple program, can put a lot more people present in the form of image data, for the enterprise management with the help of the statistical work and wages are huge. Because learning time is not long, so I

18、 am to learn knowledge, based on reference to some other books, a collection of books on the design experience of others, oneself write a small program, can quickly calculate the one real wages each month.英文關(guān)鍵詞：mage statistics wages1. 引言Matlab是一套高性能的數(shù)值計(jì)算和可視化的工具軟件。也是一套十分強(qiáng)大的計(jì)算機(jī)科學(xué)計(jì)算、可視化輔助設(shè)計(jì)和可視化的教學(xué)軟件。它的

19、應(yīng)用十分的廣泛，包括科學(xué)計(jì)算、信號(hào)和圖像處理、通訊、控制系統(tǒng)計(jì)算、測(cè)試和測(cè)量、財(cái)務(wù)建模和分析，以及計(jì)算生物學(xué)等眾多科學(xué)領(lǐng)域。日常，我們計(jì)算一些數(shù)據(jù)量小，數(shù)據(jù)個(gè)數(shù)不多的文件，可能人工計(jì)算也能很快的得出結(jié)果，但是在大數(shù)據(jù)面前，如：每個(gè)月的工資統(tǒng)計(jì)。鑰匙還是使用人力計(jì)算，就需要大量的人和時(shí)間，還存在著一定的錯(cuò)誤率，這明顯是十分沒(méi)有效率的；但是用Matlab設(shè)計(jì)一個(gè)簡(jiǎn)單的程序，就能解決所有的問(wèn)題，更能使正確率達(dá)到100%。由簡(jiǎn)至繁，我們還可以根據(jù)我們不同的需求來(lái)制定不同的程序，而要做的也僅僅只是改動(dòng)一下代碼，這樣既靈活多變，也能解決許多日常工作中繁瑣的問(wèn)題。2. 算法基礎(chǔ)條件語(yǔ)句if else；elseif 例：某IT公司員工的工資計(jì)算方法如下： （1） 工作時(shí)數(shù)超過(guò)170小時(shí)者，超過(guò)部分加發(fā)15%。 （2） 工作時(shí)數(shù)低于150小時(shí)者，扣發(fā)700元。 （3）