lec01(survey)-資料下載頁
2025-08-04 09:47本頁面
【正文】 而是 初步認(rèn)識(shí) 科學(xué)計(jì)算 Matlab, history ? Moler 是 Linpack 和 Eispack的開發(fā)者之一 ? 70年代末， Moler在教數(shù)值分析課程時(shí)，要求 學(xué)生使用EISPACK和 LINPACK程序庫 ? 發(fā)現(xiàn)很費(fèi)時(shí)間 ? 他決定自己寫接口 ? 1983年在 Stanford講學(xué) ? 吸引了工程師 John Little ? 1984年成立公司 ? Cleve Moler ? The founder of Matlab ? The chief scientist of Mathworks 通過實(shí)例了解 Matlab Fibonacci numbers Magic squares 例 1： Fibonacci 數(shù)列 斐波那契數(shù)列 (Fibonacci)”的發(fā)明者，是 意大利 數(shù)學(xué)家 列昂納多 斐波那契 （ Leonardo Fibonacci，生于公元 1170年， 卒于 1240年，籍貫大概是 比薩 ）。 1202年，他撰寫了 《 珠算原理 》 (Liber Abaci)一書。 他是第一個(gè)研究了 印度 和 阿拉伯 數(shù)學(xué)理論的 歐洲 人。 Fibonacci 數(shù)列 ? 以兔子繁殖為例子 ? 一般而言，兔子在出生兩個(gè)月后，就有繁殖能力，一對兔子每個(gè)月能生出一對小兔子來。如果所有兔都不死，那么一年以后可以繁殖多少對兔子？ ? 拿新出生的一對小兔子分析一下： ? 第一個(gè)月小兔子沒有繁殖能力，所以還是一對； ? 兩個(gè)月后，生下一對小兔民數(shù)共有兩對； ? 三個(gè)月以后，老兔子又生下一對，因?yàn)樾⊥米舆€沒有繁殖能力，所以一共是三對； ? …… ? 這個(gè)數(shù)列前兩項(xiàng)為 1，從第三項(xiàng)開始，每一項(xiàng)都等于前兩項(xiàng)之和。 Fibonacci 數(shù)列的 Matlab程序 ? N=20。 ? f=zeros(n,1)。 ? f(1)=1。 f(2)=1。 ? for k=3:n ? f(k)=f(k1)+f(k2)。 ? end Example 2 幻方 個(gè)數(shù)的和相等線所含以及兩條對角階方陣，其各行、各列個(gè)自然數(shù)組成的一個(gè)階幻方是由前 2nnnnMatlab 命令 ? N=4。 A=magic(n)。 ? sum(A) ? sum(A,2) ? sum(diag(A)) ? det(A) ? if det(A)~=0 B=inv(A)。 end 算法的差別 (向量表示 ) clear。 n=50000。 tic。 for i=1:n b(i)=sqrt(i)。 end t=toc。 disp([39。Time for loop method:39。,num2str(t)]) tic。 a=1:n。 b=sqrt(a)。 t=toc。 disp([39。Time for vector method:39。,num2str(t)]) 三維繪圖 ? x=::。 y=x39。 ? u=ones(size(y))*x。 ? v=y*ones(size(x))。 ? r=sqrt(u.^2+v.^2)。 %+eps。 ? z=sin(r)./r。 ? mesh(z)。 ? xlabel(39。x39。)。 ylabel(39。y39。)。 ? zlabel(39。z39。)。 聲音與動(dòng)畫 ? A=rand(1,10000)。 sound(A)。 ? load train。 sound(y)。 ? load chirp。 sound(y)。 ? load mddata。 ? movie(F1,3)。 movie(F2) 測驗(yàn)（ 30分鐘）
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號(hào)-1