院系計算機學院
班級2004級二學位
指導教師張興永
姓名田琳
實驗目錄
程式內容:源程式:\
程式清單:
function s(x)
s=1;
for(i=1:8);
s=s*cos(2^i*x);
ends
實驗結果:s(1)= 3.4160e-004;s(2)= 8.1827e-004;s(pi)=1; s(pi/2)=-1;s(pi/3)= 0.0039
程式內容: fourier、laplace、ztrans變換
源程式:\
程式清單:
syms t;
f=5*sin(2*t)-3*cos(2*t)
fourier_f=fourier(f)
laplace_f=laplace(f)
ztrans_f=ztrans(f)
實驗結果:
f =5*sin(2*t)-3*cos(2*t)
fourier_f =pi*(5*i*dirac(w+2)-3*dirac(w+2)-5*i*dirac(w-2)-3*dirac(w-2))
laplace_f =10/(s^2+4)-3*s/(s^2+4)
ztrans_f =
10*z*cos(1)*sin(1)/(-4*z*cos(1)^2+z^2+2*z+1)-3*(z+1-2*cos(1)^2)*z/(-4*z*cos(1)^2+z^2+2*z+1)
對矩陣進行高斯消去法變換
求三次初等變換的總的等價乘子,並用matlab求原來a的行列式、秩和跡
源程式 :\
程式清單:
a=[1 0 7;4 1 5;2 -1 9]; a0=a; %輸入a,並保留乙個備份
a(2,: ) = -4*a(1,: )+a(2,: )
a1=a,
b1=a1/a0消去a(2,1),求b1
a(3,: ) = -2*a(1,: )+a(3,: )
a2=a,
b2=a2/a1消去a(3,1)
a(3,: ) = -a(3,2)/a(2,2)*a(2,: )+a(3,: )
a3=a,
b3=a3/a2消去a(3,2)
b0 = a3/a0求三次初等變換的總的等價乘子
det_a=det(a0),
rank_a=rank(a0),
tr_a=trace(a0求原來a的行列式、秩和跡
實驗結果:
, det_a = -28 rank_a =3tr_a =11
心形線繪製:
clear
t=0:0.001:2*pi;
subplot(2,2,1);
polar(a, 1+cos(t))
subplot(2,2,2);
plot(cos(t).^3,sin(t).^3)
subplot(2,2,3);
polar(t,abs(sin(t).*cos(t)))
subplot(2,2,4);
polar(t,(cos(2*t)).^0.5)
影象:圖表 1四曲線
圖表 2 fib 曲線
實驗結果:100以內的fibonacci 數:f =
1 1 2 3 5 8 13 21 34 55 89
實驗結果:1000以內的fibonacci數:f=
1 1 2 3 5 8 13 21 34 55 89
144 233 377 610 987
資料:12 34 56 78 99 123 165 198 243 277 353 345 303 288 275的三次擬合
程式:\
圖表 3 三次擬合曲線
實驗內容: [1]p11頁試驗課題1
clear;
clc;
syms x;
resualt=zeros(4);
n_r=1;
a=input('輸入a的值:');
方案1i1=log(a+1)-log(a);
for n=1:10
i1=-1*a*i1+1/n;
f=x^n/(a+x);
i0=int(f,'x',0,1);
i0=vpa(i0,500);
i0=vpa(i0,8);
i=vpa(i1,6);
wucha1=abs((i0-i1)/i0);
wucha1=vpa(wucha1,8) ;
db=[n i0 i wucha1];
if n==1
resualt=db;
else
resualt=[resualt;db];
endendresualt
方案2resualt=zeros(4);
n=13;
if a>=n/(n+1)
i2=(2*a+1)/(2*a*(a+1)*(n+1));
else
i2=0.5*(1/((a+1)*(n+1))+1/n);
endfor n=n:-1:1
f=x^n/(a+x);
i0=int(f,'x',0,1);
i0=vpa(i0,500);
i0=vpa(i0,8);
i2=(-1*i2+1/n)/a;
i=vpa(i2,6);
wucha2=abs((i0-i2)/i0);
wucha2=vpa(wucha2,8);
db=[n i0 i wucha2];
if n==n
resualt=db;
else
resualt=[resualt;db];
endendresualt
輸入a的值:0.05
次數精確值迭代值相對誤差
1, .84777388, .847774, .22248517e-8]
2, .45761131, .457611, .85349539e-8]
3, .31045277, .310453, .63500227e-8]
4, .23447736, .234477, .68175840e-8]
5, .18827613, .188276, .10198168e-7]
6, .15725286, .157253, .44935959e-9]
7, .13499450, .134994, .10844168e-8]
8, .11825028, .118250, .42221299e-7]
9, .10519860, .105199, .25088308e-7]
10, .94740070e-1, .947401e-1, .13928927e-8]
次數精確值迭代值相對誤差
13, .72971840e-1, .889587e-1, .21908205]
12, .79024729e-1, -.112507, 2.4236878]
11, .86172087e-1, 4.06831, 46.211490]
數學實驗實驗報告
天水師範學院數學與統計學院 實驗報告 實驗專案名稱 所屬課程名稱 實驗型別線性代數實驗 實驗日期 班級學號 姓名成績 附錄1 源程式 附錄2 實驗報告填寫說明 1 實驗專案名稱 要求與實驗教學大綱一致。2 實驗目的 目的要明確,要抓住重點,符合實驗教學大綱要求。3 實驗原理 簡要說明本實驗專案所涉及...
數學實驗報告
實驗名稱 matlabe的認識與應用 學院土木與環境工程學院 專業班級土木0804 姓名張偉康 學號40812104 2010年 10月 一 實驗目的 對matlabe軟體做初步的認識與應用,學會解決一些基本的實際問題 二 實驗任務 解決老師布置的教材應用問題,在訓練中進一步加強對該軟體的應用能力 ...
數學實驗報告
計算02 1班 王莉 實驗八曲柄滑塊機構的運動規律 實驗目的 本實驗主要涉及微積分中對函式特性的研究。通過實驗複習函式求導法,taylor公式和其他有關知識。著重介紹運用建立近視似模型並進行數值計算來研究討論函式的方法。問題重述 曲柄滑塊機構是一種常用的機械結構,它將曲柄的轉動轉化為滑塊在直線上的往...