精品文档---下载后可任意编辑姓 名：学号：专业：《MATLAB 程序设计》作业1、考虑如下 x-y 一组实验数据：x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2]分别绘出 plot 的原始数据、一次拟合曲线和三次拟合曲线，给出 MATLAB 代码和运行结果。原始曲线一次拟合三次拟合12345678910123456781234567891022.533.544.555.566.571234567891012345678精品文档---下载后可任意编辑x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10];y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2];figure;plot(x,y)p1=polyfit(x,y,1);y1=polyval(p1,x);figure;plot(x,y1)p2=polyfit(x,y,3);y2=polyval(p2,x);figure;plot(x,y2)2、在[0，3π]区间，绘制 y=sin(x)曲线（要求消去负半波，即(π，2π)区间内的函数值置零），求出曲线 y 的平均值，以及 y 的最大值及其最大值的位置。给出执行代码和运行结果。x=0:pi/1000:3*pi;y=sin(x);y1=(y>=0).*y; %消去负半波 figure(1);plot(x,y1,'b');a=mean(y1) %求出 y1 的平均值b=max(y1) %求出 y1 的最大值 b,以及最大值在矩阵中的位置;d=x(find(y1==b))>> ex1a =b = 1d =>>3、给出函数 z=x^2+y^2 的三维曲面的绘制结果和 MATLAB 代码，x=-3:0.1:3;y=-4:0.1:4;[x,y]=meshgrid(x,y);z=x.^2+y.^2;mesh(x,y,z);01234567891000.10.20.30.40.50.60.70.80.91-4-2024-4-20240510152025精品文档---下载后可任意编辑4、将一幅彩色图像文件转换为灰度图像，存储并显示转换后的灰度图像，给出运行结果和执行代码。x=imread('hehua.jpg');x=rgb2gray(x);figure(1);imshow(x);h=figure(1);saveas(h,'gray.jpg');原图处理后图像5、结合自身讨论方向，自拟一题目，采纳 MATLAB 语言实现一具有完整功能的算法，给出运行结果和执行代码。题目:用遗传算法找函数的最大值.function main()clearclcpopsize = 100; %种群大小chromlength = 10; %二进制编码长度pc = 0.6; %交叉概率pm = 0.001; %变异概率pop = initpop(popsize,chromlength); %初始种群for i=1:100 [objvalue] = cal_objvalue(pop); %计算适应度值（函数值） fitvalue = objvalue; [newpop] = selection(pop,fitvalue); %选择操作 [newpop] = crossover(newpop,pc); %交叉操作 [newpop] = mutation(newpop,pm); %变异操作 pop = newpop; %更新种群 [bestindividual,bestfit]=best(pop,fitvalue);%寻找最优解 x2 = binary2decimal(bestindividual); x1 = binary2decimal(newpop);精品文档---下载后可任意编辑 [y1] = cal_objvalue(newpop); if mod(i,10)==0 figure; fplot('10*sin(5*x)+7*abs(x-5)+10',[0 10]); hold on; title(['迭代次数为 n=' num2str(i)]); plot(x1,y1,'*');endend fprintf('the best X is --->>%5.2f\n',x2); fprintf('the best Y is --->>%5.2f\n',bestfit);
