可視化計(jì)算離線作業(yè)

1、浙江大學(xué)遠(yuǎn)程教育學(xué)院可視化計(jì)算課程作業(yè)2016年（夏）姓名：賀鵬學(xué) 號(hào)：716128012033年級(jí)：16秋專本計(jì)算機(jī)學(xué)習(xí)中心：浙大校內(nèi)直屬學(xué)習(xí)中心（華家池）注意：所有圖像的標(biāo)題必須呈現(xiàn)足夠你本人信息1. （10分）求解下列線性方程組的解：2. （10分）信號(hào)y = 5*sin(pi*20*t)+3*cos(2*pi*50*t)幅度為1的白噪聲的干擾，請(qǐng)畫出此信號(hào)，然后進(jìn)行傅立葉變換，畫出變換后的頻域信號(hào)。clear;title('賀鵬716128012033')t=0:0.001:0.6;y=5*sin(pi*20*t)+3*cos(2*pi*50*t);y=y+randn(

2、1,length(t);subplot(1,2,1)plot(t,y)title('賀鵬')xlabel('變換前信號(hào)')Y=fft(y,512);subplot(1,2,2)Y=real(Y);plot(Y(1:512)title('716128012033')xlabel('變換后的頻域信號(hào)')3. （10分）在空間有一個(gè)帶正電的點(diǎn)電荷，請(qǐng)畫出此點(diǎn)電荷的空間電位分布與電場(chǎng)強(qiáng)度的空間分布圖。clear;subplot(1,2,1)x,y=meshgrid(-4:0.2:4,-4:0.2:4);z=1./sqrt(x.2+y.2+

3、0.01);mesh(x,y,z)title('賀鵬')subplot(1,2,2)x,y=meshgrid(-2:0.1:2,-2:0.1:2);z=1./sqrt(x).2+y.2+0.01);dx,dy=gradient(z);contour(x,y,z,-12,-8,-5,-3,-1,-0.5,-0.1,0.1,0.5,1,3,5,8,12);hold onquiver(x,y,dx,dy)title('716128012033')4. （10分）仿照課本第11章的太陽地球月亮衛(wèi)星，繞轉(zhuǎn)動(dòng)畫實(shí)例，呈現(xiàn)地球繞太陽運(yùn)轉(zhuǎn)的動(dòng)畫。Clear;figure(

4、9;name','地球繞太陽轉(zhuǎn)動(dòng)');title('賀鵬 716128012033')s1=0:.01:2*pi;hold on;axis equal;axis off r1=10;w1=1;t=0;pausetime=.002;sita1=0;sita2=0;set(gcf,'doublebuffer','on') plot(-20,18,'color','r','marker','.','markersize',40);text(-17,1

5、8,'太陽');p1=plot(-20,16,'color','b','marker','.','markersize',20);text(-17,16,'地球');plot(0,0,'color','r','marker','.','markersize',60);plot(r1*cos(s1),r1*sin(s1);set(gca,'xlim',-20 20,'ylim'

6、;,-20 20);p1=plot(r1*cos(sita1),r1*sin(sita1),'color','b','marker','.','markersize',30);while 1set(p1,'xdata',r1*cos(sita1),'ydata',r1*sin(sita1);sita1=sita1+w1*pausetime;pause(pausetime);drawnowEnd5. （10分）設(shè)計(jì)一個(gè)低通濾波器，從混合信號(hào)：x(t)=sin(2*pi*10*t)+co

7、s(2*pi*100*t)+0.2*randn(size(t)中獲取10Hz的信號(hào)（10分）。clear;figure('name','賀鵬716128012033');ws=1000;t=0:1/ws:0.4;x=3*sin(2*pi*10*t)+2*cos(2*pi*100*t)+0.8*randn(size(t);wn=ws/2;B,A=butter(10,30/wn);y=filter(B,A,x);plot(t,x,'b-')hold onplot(t,y,'r','MarkerSize',10)lege

8、nd('input','Output')6. （20分）設(shè)計(jì)一個(gè)程序，應(yīng)用函數(shù)subplot(1,2,1)、subplot(1,2,2)分別顯示您本人的二張照片，然后對(duì)二張照片分別進(jìn)行傅立葉變換，并分別畫出變換后的頻域信號(hào)。再把2個(gè)頻域信號(hào)相加，經(jīng)傅立葉逆變換后，顯示時(shí)域信號(hào)的圖像。clear;A=imread('C:1.jpg');subplot(1,2,1)imshow(A);title('賀鵬716128012033')B=imread('C:2.jpg');subplot(1,2,2)imshow(B);t

9、itle('賀鵬716128012033')A1=fft2(A);B1=fft2(B);subplot(1,2,1);imshow(A1)title('賀鵬716128012033')subplot(1,2,2);imshow(B1) title('賀鵬716128012033')頻域信號(hào)相加后經(jīng)傅立葉逆變換C1=A1+B1;C=ifft2(C1);imshow(C)xlabel('716128012033');ylabel('賀鵬');7. （30分）小論文 根據(jù)工作中的實(shí)際需要，請(qǐng)?jiān)O(shè)計(jì)一個(gè)實(shí)際工程問題的可視化。

10、可以選擇以下之一：（1）工程動(dòng)畫的可視化；（2）大數(shù)據(jù)處理中的可視化；（3）算法與模型計(jì)算的可視化；（4）實(shí)際生產(chǎn)流程的可視化；（5）或其它有創(chuàng)新意義的可視化科學(xué)計(jì)算。要求：（1）題目有實(shí)際意義。（2）有分析、算法描述（3）程序源代碼設(shè)計(jì)。（4）問題結(jié)果有可視化顯示。（4）題目的問題有一定的新意。小論文的字?jǐn)?shù)不能少于2500字，格式由下列各部分組成：Matlab在公司中如何用來計(jì)算工資摘要： 隨著社會(huì)的不斷發(fā)展，許多大型的企業(yè)也在迅速崛起，企業(yè)員工也從以前的幾十個(gè)人發(fā)展到成百上千，乃至數(shù)萬人?？上攵?，單純的靠人力去計(jì)算每個(gè)人都每月的實(shí)發(fā)工資肯定是不現(xiàn)實(shí)的。最近學(xué)習(xí)了Matlab以后，發(fā)現(xiàn)它不

11、僅僅能編寫一些簡(jiǎn)單的程序，更可以把許多數(shù)據(jù)以圖像的方式呈獻(xiàn)給人們，這對(duì)于企業(yè)的管理以及工資的統(tǒng)計(jì)工作有著巨大的幫助。因?yàn)閷W(xué)習(xí)時(shí)間也不是很久，所以我在根據(jù)已學(xué)的知識(shí)基礎(chǔ)上，參考了一些其他書籍，集合書本上他人的設(shè)計(jì)經(jīng)驗(yàn)，自己寫了一個(gè)小程序，可以很快的計(jì)算出一個(gè)人每月的實(shí)際工資。中文關(guān)鍵詞： 圖像 統(tǒng)計(jì) 工資英文題目：Matlab in the company how to calculate salary英文摘要： With the continuous development of society, also is in rapid rise, many large enterprises fr

12、om 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 monthly net salary is not realistic. Recent learning Matlab, found that it not only can write some simple program,

13、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 am to learn knowledge, based on reference to some other books, a collection of books on the design experience o

14、f 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é)軟件。它的應(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ù)面前，如

15、：每個(gè)月的工資統(tǒng)計(jì)。鑰匙還是使用人力計(jì)算，就需要大量的人和時(shí)間，還存在著一定的錯(cuò)誤率，這明顯是十分沒有效率的；但是用Matlab設(shè)計(jì)一個(gè)簡(jiǎn)單的程序，就能解決所有的問題，更能使正確率達(dá)到100%。由簡(jiǎn)至繁，我們還可以根據(jù)我們不同的需求來制定不同的程序，而要做的也僅僅只是改動(dòng)一下代碼，這樣既靈活多變，也能解決許多日常工作中繁瑣的問題。2. 算法基礎(chǔ)條件語句if else；elseif 例：某IT公司員工的工資計(jì)算方法如下： （1） 工作時(shí)數(shù)超過170小時(shí)者，超過部分加發(fā)15%。 （2） 工作時(shí)數(shù)低于150小時(shí)者，扣發(fā)700元。 （3） 其

16、余按每小時(shí)50元計(jì)發(fā)。 按輸入的工號(hào)和該號(hào)員工的工時(shí)數(shù)，計(jì)算應(yīng)發(fā)工資。3. 程序代碼clear;num=input('請(qǐng)輸入工號(hào):'); time=input('請(qǐng)輸入工時(shí)數(shù):');  if time<150 salary=time*50-700; elseif time>170 salary=84*170+(time-120)*84*(1+0.15); else  salary=time*50; end salary4. 結(jié)果分析（必須有可視化圖）當(dāng)工時(shí)小于150小時(shí)時(shí)：當(dāng)工時(shí)屬于150-170小時(shí)之間時(shí)：當(dāng)工時(shí)大于170小時(shí)時(shí)：結(jié)論由上所述，若某公司的工資按照每個(gè)員工每月所做工