2022年华数杯C题插层熔喷非织造材料的性能控制-完整代码.zip

Program 1 :Three spline interpolation method close all clear all clc x=[0 250 500 700]; y=[96.5465 148.44 129.214 115.6626]; cs=spline(x,[0 y 0]); xx=linspace(0,700,100); plot(x,y,'o',xx,ppval(cs,xx),'-'); Program 2 :Genetic algorithm main program clc clear all close all ts=AA(:,1); Hs=sqrt(AA(:,2).^2+AA(:,3).^2); xxx=AA(:,2); yyy=AA(:,3); [PM PN]=size(AA); for times=1:1 % for m=1:PN % A(m,PN)=0; % end n=200; ger=200; pc=0.9; pm=0.01; % [Ax,Ay]=size(A) % 生成初始种群 v=init_population(n,2); [N,L]=size(v); disp(sprintf('Number of generations:%d',ger)); disp(sprintf('Population size:%d',N)); disp(sprintf('Crossover probability:%.3f',pc)); disp(sprintf('Mutation probability:%.3f',pm)); % 计算适应度，并画出图形 [fit,xx,yy,QS]=Sun_Shadow_fun(v,ts,AA(:,2),AA(:,3)); % 初始化 vmfit=[]; it=1; vx=[]; %C=[]; % 开始进化 while it<=ger %Reproduction(Bi-classist Selection) vtemp=Sun_Shadow_roulette(v,fit); %Crossover v=Sun_Shadow_crossover(vtemp,pc); v=Sun_Shadow_mutation(vtemp,pm); %Results [fit,xx,yy,QS]=Sun_Shadow_fun(v,ts,AA(:,2),AA(:,3)); [sol,indb]=min(fit); v(1,:)=v(indb,:); media=mean(fit); vx=[vx sol]; vmfit=[vmfit media]; it=it+1; end %%%% 最后结果 disp(sprintf('\n')); %空一行 % 显示最优解及最优值 v(indb,: % 图形显示最优及平均函数值变化趋势 figure(1); plot(vx); %title('最优,平均函数值变化趋势'); xlabel('Generations'); ylabel('f(x)'); hold on; %plot(vmfit,'r'); hold off; fa(times)=180*v(1,1)-90; H(times)=2*v(1,2)+1; fit(1)/20; Fit(times)=fit(1)/sum(Hs)/2; end runtime=toc res=[fa H Fit]; [b,bint,r,rint,stats]=regress(Hs,QS(:,1)) rcoplot(r,rint) Program 3 :Polynomial fitting solution code clear; x=[0 250 500 700]; y=[96.5465 148.44 129.214 115.6261]; p=polyfit(x,y,3); px=poly2str(p,'x'); pv=polyval(p,x); plot(x,y,'*',x,pv,'k-') xlabel('x') ylabel('y') p px x_aver=sum (x)/length(x); y_aver=sum (y)/length(y); fori=1:length(x) Lxx_1(i)=(x(i)-x_aver)^2; Lyy_1(i)=(y(i)-y_aver)^2; Lxy_1(i)=x(i)*y(i)-x_aver*y_aver; end Lxx=sum(Lxx_1); Lyy=sum(Lyy_1); Lxy=sum(Lxy_1); b_guji=Lxy/Lxx; a_guji=y_aver-b_guji*x_aver; y_guji_x_200=a_guji+b_guji*200; sigamma_guji=1/(length(x)-2)*(Lyy-b_guji^2*Lxx); rhoxy_guji=Lxy/(sqrt(Lxx)*sqrt(Lyy)); disp(['R=',num2str(rhoxy_guji)]) Program 4 :Residual test code close all clear all clc x=[143 145 146 147 149 150 153 154 155 156 157 158 159 160 162 164]'; X=[ones(16,1) x]; Y=[88 87 88 91 92 93 93 95 96 98 97 96 98 99 100 102]'; [b,bint,r,rint,stats]=regress(Y,X); rcoplot(r,rint);