2016全国大学生数学建模竞赛培训(MATLAB).ppt

数学建模 算例1 contourf(X,Y,V,cv,k-) %axis(square) title(\fontname{隶书}\fontsize{22}偶极子的场),hold on quiver(X,Y,Ex,Ey,0.7) plot(a,b,wo,a,b,w+) plot(-a,-b,wo,-a,-b,w-) xlabel(x);ylabel(y),hold off 数学建模 算例1 数学建模 算例2 hist指令的使用示例。 randn(state,1),rand(state,31) x=randn(100,1);y=rand(100,1); subplot(1,2,1),hist(x,7) subplot(1,2,2),histfit(x,20) 数学建模 算例2 n_y1=min(y):0.1:max(y);n_y2=min(y):0.05:max(y); subplot(1,2,1),hist(y,n_y1) subplot(1,2,2),hist(y,n_y2) 数学建模 算例3 取发生信号的原始模型为 。 在[0, 4]中取值； 受到噪声0.3*(rand(n,1)-0.5) 的污染。 本例演示：如何以 为模型，通过fminsearch从受污染数据中，估计出参数 。 数学建模 算例3 （1）创建M文件xydata.m [XYDATA.M] function [x,y,STDY]=xydata(k_noise) %xydata.m x=[0:0.2:4]; yo=3*exp(-0.4*x)+12*exp(-3.2*x); rand(seed,234) y_noise=k_noise*(rand(size(x))-0.5); y=yo+y_noise; STDY=std(y_noise); 数学建模 算例3 创建M文件twoexps.m [TWOEXPS.M] function E=twoexps(a,x,y) %twoexps.m x=x(:); y=y(:); Y=a(1)*exp(-a(3)*x)+a(2)*exp(-a(4)*x); E=sum((y-Y).^2); 数学建模 算例3 （2）编写M脚本文件作为主文件 [EXM.M] %exm041022_1.m k_noise=0.3; [x,y,STDY]=xydata(k_noise); a0=[1 1 1 1]; options=optimset(fminsearch); options.TolX=0.01; options.Display=off; a=fminsearch(@twoexps,a0,options,x,y); chi_est=twoexps(a,x,y)/STDY^2; freedom=length(x)-length(a0); Q=1-chi2cdf(chi_est,freedo 数学建模 算例3 y_est=a(1)*exp(-a(3)*x)+a(2)*exp(-a(4)*x); ych=y_e_s_t=; a1=num2str(a(1)); a2=num2str(a(2)); a3=num2str(a(3)); a4=num2str(a(4)); char_y_est=[ych,a1,*exp(-,a3,*x)+,a2,*exp(-,a4,*x)]; plot(x,y,b+);holdon, plot(x,y_est,r); hold off, axis([0,4,0,16]) 数学建模 算例3 text(0.4,14,y=3*exp(-0.4*x)+12*exp(-3.2*x)) text(0.4,12,char_y_est),text(2.5,9,[chi2= , num2str(chi_est)]) text(2.5,7,[freedom= , num2str(freedom)]) text(2.5,5,[Q= , num2str(Q)]) 数学建模 算例3 （3） exm 作业 完成实验报告（一） * * 一行的数值是变量的个数么？？？？？ * 数学建模 一、函数 子函数 M文件中可以包含一个以上的函数，文件中除主函数以外的其他函数称为子函数，并且它们只对主函数或同一文件中的其他函数可见。 主函数放在最上面，子函数放在下面，各子函数排名不分先后。 即使在同一个M文件中，子函数也不能获取用于主函数或其他子函数的参数，除非它们声明为全局变量，或者作为参数传递。 从M文件内部调用函数时，MATLAB首先检查文件，看函数是否为子