常用算法matlab源程序.docVIP

参考目录： 微分方程的数值解 见F盘 微分方程的数值解.ppt 优化问题 111优化问题matlab27.doc 优化工具箱简洁.ppt 辅助优化计算与设计.pdf 优化建模与LINDO/LINGO软件/谢金星, 薛毅编著 51.81/3486 2005- 南中文借书处 可借 在馆 ###################################################### 曲线拟合及拟合后相应点的值 clear x=0:0.1:1; y=[-.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2]; n=2; p=polyfit(x , y , n); xi=linspace(0 , 1 , 100); z=polyval(p , xi); plot(x , y , -.ro , xi , z , : ) xlabel(x) ; ylabel(y=f(x)) title(Second Order Curve Fitting) legend(points,line)%图例 ###################################################### %plot画图的例子 x=0:pi/20:2*pi; y1=x.*sin(x); plot(x,y1,-.r*)%显示每一点的用*，并用虚线连起来 hold on %还要在这个加其他的图 y2=exp(x/100).*sin(x-pi/2); plot(x,y2,g) %用实线条画图 plot(x,y2,-mo) %用实线条画图 plot(x,y2,-mo) y3=sin(x-pi); plot(x,y3,gs) title(y1=..;y2=.;y3=..;)%加标题 legend(y1,y2,y3)%加上图利 xlabel(x axis biaoshi....);ylabel(y axis biaoshi..) axis([0 2*pi -5 2])%调整坐标轴 grid on %加上网格 %要想加上箭头可以手动 ##################################################################### 1.三维立体网格图mesh clear x=-5:0.01:5; y=-5:0.01:5; z=x.^2+y.^2; [X,Y]=meshgrid(x,y); Z=X.^2+Y.^2; figure,mesh(X,Y,Z) title(san wei surface) xlabel(x axis) ylabel(y axis) zlabel(z axis) axis([-5 5 -5 5 0 100]) grid on %画坐标方格 colorbar%加上彩色比照图 %%%%%%%%%%%% 其中将mesh改为meshz （后加colormap([1 0 0]) axis tight)和 waterfall(后家colormap axis tight) 分别可以得到 设立一个账莫平台 其中将mesh改为meshc画立体图提取登高线 contour提取登高线 contour3(Z)立体的登高线 ####################################################################################### 2.三维线和点 clear t=1:0.01:15; x=2*(exp(-0.2.*t).*cos(8.*t)); y=2*(exp(-0.4.*t).*cos(8.*t)); plot3(t,x,y) %标注 xlabel(x axis) ylabel(y axis) zlabel(z axis) axis([-5 5 -5 5 0 100]) ######################################################%灰色预测的matlab程序 %x0为原始值（ 输入六个原始数据） %f为预测值 得到8个数据（后两个为预测数据）精确度很高 %e为残差 function [f,e]=GM(x0) x1 = cumsum(x0); B = zeros(5,1); for i =1:5 B(i) = -(x1(i) + x1(i+1)) /2; end for i =1:5 y(i)