《matlab中遇到的一些问题》由会员分享，可在线阅读，更多相关《matlab中遇到的一些问题（19页珍藏版）》请在金锄头文库上搜索。
1、1什么是Gaussian拟合？:#FangQ(Qianqian.FangDartmouth.Edu),2002/9/23, BigGreen/MathTools # 高斯拟合(Gaussian Fitting)即使用形如： Gi(x)=Ai*exp(x-Bi)2/Ci2) 的高斯函数对数据点集进行函数逼近的拟合方法。 其实可以跟多项式拟合类比起来，不同的是多项式拟合是用幂函数系， 而高斯拟合是用高斯函数系。 使用高斯函数来进行拟合，优点在于计算积分十分简单快捷。这一点 在很多领域都有应用，特别是计算化学。著名的化学软件Gaussian98 就是建立在高斯基函数拟合的数学基础上的。2如何在给定句
2、柄的axis里绘图?:#FangQ(Qianqian.FangDartmouth.Edu),2002/6/12, SMTH/MathTools # plot(data,parent,haxis); 或者 hbar=bar(data); set(hbar,parent,haxis);3)由Matlab符号运算得到的公式怎么才能将数据代进去运算？:#ramjet (德芙)2002/3/3, SMTH/MathTools # 使用subs(),或先将值赋予一个符号变量，然后用eval()5)在Matlab中如何求最值点？如何求一维数组的极值？:#FangQ(Qianqian.FangDartmout
3、h.Edu),2002/6/18, SMTH/MathTools# 最值： 一维或多维数组最值用max(data(:) 如果想返回最值所在的位置，用Y,I=max(data):#FangQ(Qianqian.FangDartmouth.Edu), 2001/4/21,UESTC/Math# 极值： data是你的数据， find(diff(sign(diff(data)=-2)+1 找到极大值的位置 find(diff(sign(diff(data)=2)+1 找到极小值的位置 data(find(diff(sign(diff(data)=-2)+1)和 data(find(diff(sign
4、(diff(data)=2)+1) 返回的是极大值和极小值6)Matlab中如何作线性拟合/线性回归/多元线性回归？:#FangQ(Qianqian.FangDartmouth.Edu),2002/6/21, BigGreen/MathTools # 即用y=a*x+b来拟合一组数据x1,y1,x2,y2xn,yn matlab中使用polyfit x=data(:,1); y=data(:,2); p=polyfit(x,y,1); p(1)为斜率a，p(2)为截距b 多元线性回归即用y=a1*x1+a2*x2+.+am*xm来拟合数据点x1i,x2i,xmi,yi (i=1n) |x11,
5、x21,xm1| A=|x12,x22,xm2| | | |x1n,x2n,xmn| Y=y1,y2,y3,yn 则系数a1,a2,am=pinv(A)*Y 在matlab中使用 coeff=AY 则可以得到最小二乘意义上的拟合系数7)Matlab中如何作圆回归？:#Peter Boettcher (boettcherll.mit.edu),2002/5/16, comp.soft-sys.matlab# Q5.5: How can I fit a circle to a set of XY data? = An elegant chunk of code to perform least-s
6、quares circle fitting was written by Bucher Izhak and has been floating around the newgroup for some time. The first reference to it that I can find is in: function xc,yc,R,a = circfit(x,y) %CIRCFIT Fits a circle in x,y plane % % XC, YC, R, A = CIRCFIT(X,Y) % Result is center point (yc,xc) and radiu
7、s R.A is an % optional output describing the circles equation: % % x2+y2+a(1)*x+a(2)*y+a(3)=0 % by Bucher izhak 25/oct/1991 n=length(x); xx=x.*x; yy=y.*y; xy=x.*y; A=sum(x) sum(y) n;sum(xy) sum(yy). sum(y);sum(xx) sum(xy) sum(x); B=-sum(xx+yy) ; -sum(xx.*y+yy.*y) ; -sum(xx.*x+xy.*y); a=AB; xc = -.5*
8、a(1); yc = -.5*a(2); R = sqrt(a(1)2+a(2)2)/4-a(3); Tom Davis provided a more sophisticated approach that works for more cases in and Code included.8)Matlab中如何绘制箭头？:#FangQ(Qianqian.FangDartmouth.Edu),2002/6/21, SMTH/MathTools # 到http:/ 2-D Plotting and Graphics中查找arrow.m，或者 http:/ http:/www.math.umd.
9、edu/jec/matcomp/matcompmfiles/mfiles.html9)Matlab中如何作二维数据的插值?:#FangQ(Qianqian.FangDartmouth.Edu),2002/6/21, BigGreen/MathTools # 对于一维、二维、三维规则数据点阵使用interp1/interp2/interp3， 二维、三维非规则数据用griddata/griddata310)Matlab中如何绘制三维数据阵？:#FangQ(Qianqian.FangDartmouth.Edu),2002/6/21, BigGreen/MathTools # 如果使用matlab，
10、打开帮助窗口，在目录树上找到 MATLABUsing Matlab 3-D Visualization: Volume Visualization Techniques 如果图形复杂，建议使用Tecplot，参见Tecplot手册中数据格式，将你 的三维数据读入Tecplot，双击zone，可以设置mesh/contour/surface transparency等。 在Field菜单中有3D Iso-surface Details和3D Slice Details，可以绘制等值 面和任意平面的截面图。11)Matlab中如何注解一大段代码？:#hyphone,2002/7/6, SMTH/M
11、athTools # 注释大段代码选中代码，Ctrl+R；取消注释，选中代码，Ctrl+T。 或者用Edit菜单或者右键弹出中的注释。:#misc,2002/6/21, SMTH/MathTools # if(0) 大段的代码 end12)Matlab中如何计算程序运行的时间？:#misc,2002/6/21, SMTH/MathTools # tic your_code; toc 或者使用 t=cputime; your_operation; cputime-t13)Matlab中如何改变默认的工作路径？:#SindyGong, 2002/4/7, SMTH/MathTools # 编辑一个startup.m文件，其中cd yourpath 或者在X:matlabtoolboxlocalmatlabrc.m的最后添加cd yourpath 参见： http:/
