數(shù)學(xué)實(shí)驗(yàn)作業(yè)匯總

千里之行，始于足下讓知識(shí)帶有溫度。第第2頁/共2頁精品文檔推薦數(shù)學(xué)實(shí)驗(yàn)作業(yè)匯總（1）產(chǎn)生一個(gè)5階魔方矩陣M：M=magic(5)

（2）將矩陣M的第3行4列元素賦值給變量t：t=M(3,4)

（3）將由矩陣M第2，3，4行第2,5列構(gòu)成的子矩陣賦給變N：N=M(2:4,2:3:5)

（4）將由矩陣M的前3行賦給變量N：N=M(1:3,:)

（5）將由矩陣M的后3列賦給變量N：N=M(:,end:-1:end-2)

（6）提取M的主對(duì)角線元素，并以這些對(duì)角線元素構(gòu)成對(duì)角矩陣N：N=diag(diag(M))或N=tril(triu(M))(7)隨機(jī)產(chǎn)生1000個(gè)100以內(nèi)的整數(shù)賦值給變量t：t=round(rand(1,1000)*100)

（8）隨機(jī)產(chǎn)生100*5個(gè)100以內(nèi)的實(shí)數(shù)賦值給變量M：M=rand(100,5)*100

（1）刪除矩陣M的第7個(gè)元素M(7)=[]

（2）將含有12個(gè)元素的向量t轉(zhuǎn)換成3*4的矩陣：reshape(t,3,4)

（3）產(chǎn)生和M同樣大小的單位矩陣：eye(size(M))

（4）尋覓向量t中非零元素的下標(biāo)：find(t)

（5）逆序顯示向量t中的元素：t(end:-1:1)

（6）顯示向量t偶數(shù)位置上的元素：t(2:2:end)

（7）利用find函數(shù)，將向量t中小于10的整數(shù)置為0：t(find(t>M=[t;t+10;t+20;t+30;t+40]

M=

12345678910

11121314151617181920

21222324252627282930

31323334353637383940

41424344454647484950

>>N=M(2:4,2:9)

N=

1213141516171819

2223242526272829

3233343536373839

>>sum(sum(M))-sum(sum(n))

ans=

663

2）n取100、1000、10000，求序列1、1/2、1/3……1/n的和。

>>n=100;

>>t=[1:n];

>>formatrat

>>M=t.^-1;

>>S=sum(M)

S=

2630/507

>>n=1000;

>>t=[1:n];

>>formatrat

>>M=t.^-1;

>>S=sum(M)

S=

1804/241

>>n=10000;

>>t=[1:n];

>>formatrat

>>M=t.^-1;

>>S=sum(M)

S=

1106/113

1.在同一坐標(biāo)系下繪制y1=sin(t),y2=sin(2t),y3=sin(3t),其中y1的數(shù)據(jù)點(diǎn)用星號(hào)，線形為黑色虛線，

y2的數(shù)據(jù)點(diǎn)用方塊，線形為紅色實(shí)線，y3的數(shù)據(jù)點(diǎn)用小圓圈，線形為藍(lán)色點(diǎn)線。（要求采納一次繪出和逐次填加兩種方式完成繪圖）

>>t=linspace(0,2*pi,100);

>>y1=sin(t);

>>y2=sin(2*t);

>>y3=sin(3*t);

>>plot(t,y1,’*k:’,t,y2,’sr-’,t,y3,’ob-.’)

>>t=linspace(0,2*pi,100);

>>y1=sin(t);

>>plot(t,y1,’*k:’)

>>holdon

>>y2=sin(2*t);

>>plot(t,y2,’sr-’)

>>holdon

>>y3=sin(3*t);

>>plot(t,y3,’ob-.’)

>>holdoff

2.分離用plot和fplot函數(shù)繪制y=sin（1/x）的曲線，分析兩曲線的差別

>>x=linspace(0,1/(2*pi),100);

>>y=sin(x.^-1);

>>plot(x,y,’*-’)

>>fplot(’sin(x.^-1)’,[0,1/(2*pi)],’o-’)

兩曲線的差別:plot曲線在確定自變量x的取值間隔時(shí)采納平均間隔，圖像不是非常精確?????；fplot曲線自動(dòng)取值，在函數(shù)值變化平穩(wěn)時(shí)，它的數(shù)值點(diǎn)會(huì)自動(dòng)相對(duì)稀疏一點(diǎn)，在函數(shù)值變化強(qiáng)烈處，所取點(diǎn)會(huì)自動(dòng)密集一點(diǎn)，所以曲線越發(fā)光潔精確?????。

6.

已知曲面方程f（x,y）=

,x∈[-1.5π，1.5π]，y∈[-2.5π，2.5π]，用建立子窗

口的辦法在同一圖形窗口繪制出三維線圖，網(wǎng)線圖，曲面圖。>>x=-1.5*pi:pi/50:1.5*pi;>>y=-2.5*pi:pi/50:2.5*pi;>>[X,Y]=meshgrid(x,y);

>>Z=sin(sqrt(X.^2+Y.^2))./sqrt(1+X.^2+Y.^2);>>subplot(1,3,1);plot3(X,Y,Z);>>subplot(1,3,2);mesh(X,Y,Z);>>subplot(1,3,3);surf(X,Y,Z);

8.將peaks函數(shù)生成的最高峰削去，并用色圖矩陣“cool”修飾。>>[x,y,z]=peaks(30);>>x1=x(1,:);y1=y(:,1);>>i=find(y1>1>>j=find(x1>-1>>z(i,j)=NaN*z(i,j);>>surf(x,y,z)>>colormap(cool)

3.定義一個(gè)函數(shù)，函數(shù)的自變量為整數(shù)n，函數(shù)的功能是：隨機(jī)產(chǎn)生n個(gè)三位整數(shù)，將其中小于平均值

的數(shù)用0代替。function[mean,x]=ff(n)x=floor(100+899*rand(1,n));m=length(x);mean=sum(x)/m;x(x>t=rand(1,20);

>>disp('max=');disp(max(t))

max=

0.7942

>>disp('min=');disp(min(t))

min=

0.0503

函數(shù)文件

functionf3(n)

t=rand(1,n);

disp('max=');disp(max(t));disp('min=');disp(min(t));

end

3.求下列函數(shù)的一階和二階導(dǎo)數(shù)

()

()cos

2

y=

tan3x

x

+

>>symsx

>>diff(2/tan(x)+cos(x)/3,’x’,1)

ans=

-sin(x)/3-(2*(tan(x)^2+1))/tan(x)^2

>>symsx

diff(2/tan(x)+cos(x)/3,’x’,2)

4.求積分

>>symsx

int(sqrt(exp(x)+1),x)

ans=

2*(exp(x)+1)^(1/2)+2*atan((exp(x)+1)^(1/2)*i)*i

5.求下列級(jí)數(shù)的和

>>symsn

>>s=symsum((-1)^(n+1)*1/n,1,inf)

s=

log(2)

6.求函數(shù)在x=0處的泰勒綻開式

>>symsx

>>taylor((exp(x)+exp(-x))/2,x,5,0)

ans=

x^4/24+x^2/2+1

1.利用randn函數(shù)聲稱符合正態(tài)分布的10*5隨機(jī)矩陣A，舉行以下操作：

(1).A的各列元素的均值和標(biāo)準(zhǔn)方差

(2).A的最大元素及其所在位置

(3).A的每行元素的和以及所有元素之和

(4).分離對(duì)A的每行元素按升序排序

(5).將A中的每行元素的總和按從大到小的挨次存入line_sum中，相應(yīng)的行號(hào)存入line_num中

>>A=randn(10,5);

>>a1=mean(A)

>>a2=std(A)

>>AA=max(max(A))

>>[ij]=find(A==AA)

>>a3=sum(A,2)

>>a4=sum(sum(A))

>>a5=sort(A,2)

>>[line_sum,line_num]=sort(sum(A,2),'descend')

2、補(bǔ)充題：

利用導(dǎo)入向?qū)Вɑ蚪柚瘮?shù)imread）導(dǎo)入一幅單色圖片存入變量ima_data中，然后依次完成下列操作：（1）用imshow函數(shù)顯示圖片；（2）刪除圖片前若干行（例如前100行）再次顯示該圖片。

（3）將圖片上、下翻轉(zhuǎn)再次顯示圖片。

先找到一個(gè).bmp的文件，把它放入工作名目下，并修改名稱為‘1.bmp’,執(zhí)行下列操作。

ima_data=imread(’1.bmp’);

(1)imshow(ima_data);

(2)a=ima_data(101:end,:);imshow(a);

(3)imshow(flipud(ima_data));

3.下表所示是0~90度內(nèi)某些數(shù)的正弦近似值

利用線性、樣條差值求x=20、40、80度時(shí)正弦值，這兩種辦法哪個(gè)好？為什么

試驗(yàn)步驟：利用inerp1函數(shù)先分離求出線性插值和三次樣條插值所得到的y11和y12，再利用sin（x）函數(shù)得到精確?????的y1，比較y11和y1，y12和y12，不難得出結(jié)論。

所用語句clear;clc;

x=[0153045607590]./180.*pi;

y=sin(x);

x1=[204080]./180.*pi;

y11=interp1(x,y,x1,’linear’);

y12=interp1(x,y,x1,’spline’);

y1=sin(x1);

主要結(jié)果y11=0.33920.63810.9773;

y12=0.34200.64280.9849;

y1=0.34200.64280.9848;

4.已知某次試驗(yàn)測(cè)得數(shù)據(jù)如下：

（1）請(qǐng)用3次多項(xiàng)式舉行擬合，并給出擬合函數(shù)在0、0.5、1、1.5^9、9.5處的值

（2）估量用幾階多項(xiàng)式擬合的效果較好，并說明理由。

4.(1)clear;clc;

x=1:0.4:9.4;

y=[0.870.525.213.5114.2919.4314.1341.5313.9158.5614.99130.4744.8221.2543.15281.25200.09177.93344.53509.84531.07260.49];

x1=0:0.5:9.5;

p=polyfit(x,y,3);

y1=polyval(p,x1);

主要結(jié)果:y1=[50.5533.0318.918.381.61-1.230.055.6215.6530.3249.8074.28103.92

138.91179.41225.61277.67335.79400.12470.85]

(2)19階擬合效果最好。理由通過編寫差方和函數(shù)(基于最小二乘原理）f(n)

f(n)函數(shù)如下：

functiontz=f(n)

t=[];

x=1:0.4:9.4;

y=[0.870.525.213.5114.2919.4314.1341.5313.9158.5614.99130.4744.8221.2543.15281.25200.09177.93344.53509.84531.07260.49];

fori=1:n

p=polyfit(x,y,i);

y1=polyval(p,x);

c=sum((y-y1).^2,2);

t=[tc];

end

tz=find(t==min(t));

令n=22（一共22組數(shù)據(jù)）f函數(shù)值最小時(shí)是19階時(shí)

所以得出結(jié)論19階多項(xiàng)式擬合效果最好。

再用擬合圖像（p=polyfit(x,y,19)，plot(x,y,’:o’,x,polyval(p,x),’-*’)）也可以看出19階多項(xiàng)式擬合效果最好。

2、自行練習(xí)題。下列填空題是期中考試出錯(cuò)比較多的題目，請(qǐng)仔細(xì)考慮并上機(jī)調(diào)試。

（6）逆序顯示向量t中的元素：

（7）顯示向量t偶數(shù)位置上的元素：

（9）刪除向量t中最小的5個(gè)數(shù)：

（17）將1~50按列優(yōu)先存放到5*10的矩陣M中：

（18）求矩陣M最大值所在的位置：

（19）統(tǒng)計(jì)字符串S中小寫字母的個(gè)數(shù)：

（20）設(shè)A是n階0、1方陣，A邊界上1的個(gè)數(shù)：

(6).t(end:-1:1)

(7).t(2:2:end)

(9).M=sort(t)

a=find(t=’a’4,3,-1,2;8,3,-3,4;3,3,-2,-2];>>b=[4,6,12,6]’;>>inv(A)*b運(yùn)用左除運(yùn)算符

>>A=[2,2,-1,1;4,3,-1,2;8,3,-3,4;3,3,-2,-2];>>b=[4,6,12,6]’;>>x=A\b運(yùn)用矩陣分解

>>A=[2,2,-1,1;4,3,-1,2;8,3,-3,4;3,3,-2,-2];>>b=[4,6,12,6]’;>>[Q,R]=qr(A);>>x=R\(Q\b)

4.在區(qū)間[30,50]內(nèi)，求

3()5sin()2log1.8fxxx=-+的零點(diǎn)。

>>f=’5*sin(x)-2*(log(x)/log(3))+1.8’;>>ezplot(f,30,50)>>fzero(f,33)ans=32.5547>>fzero(f,34)ans=33.3960>>fzero(f,38)ans=39.0426

>>fzero(f,[39.4，39.5])ans=39.4785則方程有四個(gè)零點(diǎn)6.給出試驗(yàn)數(shù)據(jù)如下：

試分離用

bx

b

yaeyax

==+

和做擬合形式，求出a和b及擬合曲線，并畫圖舉行比較。>>x=[2:16];

>>y=[6.24,8.20,9.58,9.60,9.60,10.02,9.93,9.99,10.47,10.59,10.60,10.80,10.60,10.90,10.75];>>X=1./x;>>Y=log(y);>>P=polyfit(X,Y,1)P=

-1.15522.4629>>exp(2.4629)ans=11.7388

則a=11.7388b=-1.1552作圖：

>>Y1=polyval(P,X)>>y1=exp(Y1);>>plot(x,y,’:o’,x,y1,’-*’)>>x=[2:16];

>>y=[6.24,8.20,9.58,9.60,9.60,10.02,9.93,9.99,10.47,10.59,10.60,10.80,10.60,10.90,10.75];>>Y=1./y;>>X=1./x;

>>P=polyfit(X,Y,1)P=

0.13840.0815則a=0.0815b=0.1384作圖：

>>Y1=polyval(P,X);>>y1=1./Y1;

>>plot(x,y,’:o’,x,y1,’-*’)

3.求下列方程或方程的根在指定點(diǎn)的近似根

23

sin()ln703210

50yxyzxzxyz?++-=?+-+=??++-=?

，初值0001,1,1xyz===functionf=myFun(x)

f(1)=sin(x(1))+x(2)^2+log(x(3))-7;f(2)=3*x(1)+2^x(2)-x(3)^3+1;f(3)=x(1)+x(2)+x(3)-5;>>X=[1,1,1]’;

>>op=optimset(’display’,’off’);>>x=fsolve(@myfun,X,op)

x=0.59912.39592.00502.已知

2sincos2(02)yxxxπ=+≤≤，求y的單調(diào)增區(qū)間和y的極值

>>fplot(’2*sin(x)+cos(2*x)’,[0,pi/2])>>symsx

>>f=2*sin(x)+cos(2*x);>>s=diff(f)s=

2*cos(x)-2*sin(2*x)

>>fzero(’2*cos(x)-2*sin(2*x)’,0.5)ans=0.5236

由圖知單調(diào)遞增區(qū)間為[0,0.5236]；將ans的值代入原式中，得y的極值為1.5。3.求解線性約束最優(yōu)化問題

2212121212121212min()0.526222..230,0

fxxxxxxxxxxxstxxxx=++≤??-+≤??

+≤??≥≥?

functionf=fop(x)

f=0.5*x(1)^2+x(2)^2-x(1)*x(2)-2*x(1)-6*x(2);>>x0=[0.5;0.5];>>A=[1,1;-1,2;2,1];>>b=[2;2;3];>>lb=[0;0];

>>options=optimset(’display’,’off’);

>>[x,f]=fmincon(@fop,x0,A,b,[],[],lb,[],[],options)x=0.66671.3333f=-8.2222

1、請(qǐng)你構(gòu)造一個(gè)生成素?cái)?shù)的公式,并將你的工作與Euler的工作比較。采納素?cái)?shù)生成公式p=n^2-79*n+1601

(1)編寫函數(shù)f(x),用來計(jì)算素?cái)?shù)多項(xiàng)式生成公式，在100以內(nèi)和1000以內(nèi)，產(chǎn)生素?cái)?shù)的百分比,程序如下：functiontz=f(x)n=0:x(1,3);

t=n.^2+x(1,1)*n+x(1,2);t1=find(isprime(t));

tz=length(t1)/length(n);

end

(2)代入Euler公式系數(shù)x1=[141100],x2=[1411000]與p=n^2-79*n+1601系數(shù)y1=[-791601100],y2=[-7916011000]比較

得到結(jié)果

f(x1)=0.8614;f(x2)=0.5814;

f(y1)=0.9505;f(y2)=0.6014;

所以可得結(jié)論該公式比Eluer的公式生成素?cái)?shù)的概率要高；

2、討論百萬以內(nèi)素?cái)?shù)的間隔邏輯。

a=primes(1000000);

b=a;b(1)=[];a(length(a))=[];

t=b-a;

plot(a,t,’.’);

t1=unique(t)%求相鄰素?cái)?shù)間的間隔值

t1=

Columns1through14

12468101214161820222426

Columns15through28

2830323436384042444648505254

Columns29through42

5658606264666870727476788082

Columns43through52

84868890929698100112114

s=zeros(2,length(t1));

fori=1:length(t1)

s(1,i)=t1(i);s(2,i)=length(find(t==t1(i)));

end

disp(s)%統(tǒng)計(jì)間隔重復(fù)的次數(shù)

Columns1through7

124681012

18169814313549556970798005

Columns8through14

14161820222426

4233288149092401217226821175

Columns15through21

28303234363840

12341914550557767330424

Columns22through28

42444648505254

47620215519610677140

Columns29through35

56586062646668

53549616244813

Columns36through42

70727476788082

22131261335

Columns43through49

84868890929698

6414121

Columns50through52

100112114

211

max(t1)%求最大間隔值

ans=114

間隔邏輯:百萬以內(nèi)相鄰素?cái)?shù)間隔值有52個(gè),其中間隔值2,4,6,8,10，12重復(fù)的次數(shù)較多,最大間隔值為114;另外10000以內(nèi)最大間隔值為36,100000以內(nèi)最大間隔值為72，所以隨著整數(shù)范圍的擴(kuò)大,最大間隔值也隨著擴(kuò)大。

1、若在構(gòu)造Koch曲線的過程中將向量CE繞點(diǎn)C逆時(shí)針旋轉(zhuǎn)90度，并作出迭代三次的分形圖。functionq=koch(p)

q=[];

t=90*pi/180;

M=[cos(t),-sin(t);sin(t),cos(t)];

fori=1:length(p)-1

A=p(:,i);B=p(:,i+1);

C=A/3*2+B/3;

E=A/3+B/3*2;

D=C+M*(E-C);

q=[q,A,C,D,E,B];

end

p=[0,1;0,0];

q=koch(koch(koch(p)));

plot(q(1,:),q(2,:))

axis([0100.6])

title(’迭代三次的koch曲線’)

2、修改Sierpinski三角形的生成元，使其不使用中點(diǎn)而用一個(gè)三等份點(diǎn)，黑色的三角形調(diào)節(jié)為隨機(jī)色彩的三角形，并作出迭代四次的分形圖。

functionq=sierpinsk(p)

q=[];

fori=1:3:length(p)

A=p(:,i);B=p(:,i+1);C=p(:,i+2);

D=A/3*2+B/3;E=B/3*2+C/3;F=C/3*2+A/3;

q=[q,A,D,F,B,E,D,C,F,E];

end

functionviewsierpinsk(p)

holdon

fori=1:3:length(p)

fill(p(1,i:i+2),p(2,i:i+2),rand());

end

holdoff

clf

pol=[-1,1,0;0,0,sqrt(3)];

q=sierpinsk(sierpinsk(sierpinsk(sierpinsk(pol))));

viewsierpinsk(q)

3、參考圖10-4，分析Minkowwski“香腸”的生成元，并作出迭代五次的分形圖。functionq=minkowwsk(p)

q=[];

t=90*pi/180;

M=[cos(t),-sin(t);sin(t),cos(t)];

N=[cos(-t),-sin(-t);sin(-t),cos(-t)];

fori=1:length(p)-1

A=p(:,i);B=p(:,i+1);

C=A/4*3+B/4;

E=(A+B)/2;

G=A/4+B/4*3;

D=C+M*(E-C);

F=E+N*(G-E);

H=E+N*(C-E);

J=G+M*(E-G);

q=[q,A,C,D,H,E,F,J,G,B];

end

p=[0,1;0,0];

q=minkowwsk(minkowwsk(minkowwsk(minkowwsk(minkowwsk(p)))));

plot(q(1,:),q(2,:))

2.對(duì)于logistic映射，選取適當(dāng)?shù)腶，使迭代序列進(jìn)入3,4,5,6周期，并給出周期軌道所用函數(shù)：

functiony=logistic(a,x0,n)

f=@(x)a*x*(1-x);

y=[];

fori=1:n

y=[y,x0];

x0=f(x0);

end

x=[];y=[];

fora=0:0.02:4

x0=0.2;f=@(x)a*x*(1-x);

fori=1:50

x0=f(x0);

end

fori=1:50

x0=f(x0);

end

fori=1:100

x0=f(x0);x=[x,a];y=[y,x0];

end

end

plot(x,y,’.’)

所用辦法：首先用logistic函數(shù)來生成迭代序列，第二構(gòu)造函數(shù)生成feigenbaum圖，然后通過調(diào)節(jié)a的取值范圍來觀看圖中周期分布并取近似值并一一試行。

所得結(jié)果：logistic(3.84,0.02,100)（即a=3.84可使迭代序列進(jìn)入3周期）

周期軌道：0.48800.95950.1494

logistic(3.46,0.02,100)（即a=3.46可使迭代序列進(jìn)入4周期）

周期軌道：0.83890.46750.86130.4132

logistic(3.74,0.02,100)（即a=3.74可使迭代序列進(jìn)入5周期）

周期軌道：0.65720.84250.49620.93490.2275

logistic(3.628,0.02,100)（即a=3.628可使迭代序列進(jìn)入6周期）

周期軌道：0.77050.64150.83440.50140.90700.3060

2、對(duì)于1000之內(nèi)的n，求Mersenne數(shù)Mn=2n-1是素?cái)?shù)的最大的n及對(duì)應(yīng)的Mersenne素?cái)?shù)的位數(shù)。只給出結(jié)果

對(duì)于1000之內(nèi)的n，Mersenne數(shù)Mn=2n-1