Matlab程序設(shè)計作業(yè)

Matlab程序設(shè)計作業(yè)

姓名:

學(xué)號：

專業(yè):

《MATLAB程序設(shè)計》作業(yè)

1、考慮如下x-y一組實驗數(shù)據(jù)：

x=[l,2,3,4,5,6,7,8,9,10]

y=[1.2,3,4,4,5,4.7,5,526,7.2]

分別繪出plot的原始數(shù)據(jù)、一次擬合曲線和三次擬合曲線，給出

MATLAB代碼和運行結(jié)果。

7

6

5

4

3

2

1_____________________________________________________________

12345678910

原始曲線

7

6.5

6

5.5

5

4.5

4

3.5

3

2.5

2一

12345678910

一次擬合

8i-------------------------------------------，----------------------------1---------------------------i---------------------------i---------------------------.---------------------------

7

6-

5

4

3

2

1------------------------------------

12345678910

三次擬合

x=[l,2,3,4,5,6,7,8,9,10];

y=[l.2,3,4,4Z5,4.7,5,5.2,6,7.2];

figure;

plot(x,y)

pl=polyfit(x,y,1);

yl=polyval(pl,x);

figure;

plot(x,yl)

p2=polyfit(x,y,3);

y2=polyval(p2,x);

figure;

plot(x,y2)

2、在［0,3捫區(qū)間，繪制y=sin(x)曲線(要求消去負半波,即(兀，2兀)區(qū)

間內(nèi)的函數(shù)值置零)，求出曲線y的平均值，以及y的最大值及其最大

值的位置。給出執(zhí)行代碼和運行結(jié)果。

x=0:pi/1000:3*pi;

y=sin(x);

yl=(y>=0).*y;當(dāng)消去負半波figure(1);

plot(x,yl,fbT);

a=mean(yl)省求出yl的平均值

b=max(yl)%求出yl的最大值b,以及最大值在矩陣中的位置；

d=x(find(yl==b))

?exl

a=

0.4243

b=

1

d=

1.57087.8540

?

1

0.9

0.8

0.7?

0.6

0.5

/\

I1\

0.43-/1\-

0.2

1

0.1

?

0

012345678910

3、給出函數(shù)z=xA2+yA2的三維曲面的繪制結(jié)果和MATLAB代碼,

x=-3:0.1:3;

y=-4:0.1:4;

[x,y]=meshgrid(x,y);

z=x.A2+y.A2;

mesh(x,y,z);

4、將一幅彩色圖像文件轉(zhuǎn)換為灰度圖像，存儲并顯示轉(zhuǎn)換后的灰度圖

像，給出運行結(jié)果和執(zhí)行代碼。

x=imread('hehua.jpg*);

x=rgb2gray(x);

figure(1);

imshow(x);

h=figure(1);

saveas(h,*gray.jpg');

原圖

處理后圖像

5、結(jié)合自身研究方向，自擬一題目，采用MATLAB語言實現(xiàn)一具有

完整功能的算法，給出運行結(jié)果和執(zhí)行代碼。

題目：用遺傳算法找函數(shù)的最大值.

functionmain()

clear

clc

popsize=100;%種群大小

chromlength=10;帶二進制編碼長度

pc=0.6;%交叉概率

pm=0.001;告變異概率

pop=initpop(popsizezchromlength);%初始種群

fori=l:100

[objvalue]=cal_objvalue(pop);3計算適應(yīng)度值(函數(shù)值)

fitvalue=objvalue;

[newpop]=selection(pop,fitvalue);%選擇操作

[newpop]=crossover(newpopzpc);會交叉操作

[newpop]=mutation(newpop,pm);%變異操作

pop=newpop;為更新種群

[bestindividual,bestfit]=best(pop,fitvalue);%尋?找最優(yōu)解

x2=binary2decimal(bestindividual);

xl=binary2decimal(newpop);

[yl]=cal_objvalue(newpop);

ifmod(ir10)==0

figure;

fplot(*10*sin(5*x)+7*abs(x-5)+10*,[010]);

holdon;

title(L迭代次數(shù)為n=1num2str(i)]);

ff

plot(xlzylz);

end

end

fprintf('thebestXis-->>%5.2f\n*Ax2);

fprintf('thebestYis-->>%5.2f\n1,bestfit);

迭代次數(shù)為n=20

60

50

”卜、T\

40

30

20

10

0

012345678910

迭代次數(shù)為n=30

60

50

40

30

20

10

0

012345678910

迭代次數(shù)為n=40

60

六

50.

卡\*

:

40

30

20

10

0

012345678910

迭代次數(shù)為n=50

60

50

40

30

20

10

0

0123456789