温馨提示！由于新浪微博认证机制调整，您的新浪微博帐号绑定已过期，请重新绑定！&&|&&
乐于交友,乐于交往够意思的朋友,乐于交往所有乐于交往的朋友。
LOFTER精选
网易考拉推荐
用微信&&“扫一扫”
将文章分享到朋友圈。
用易信&&“扫一扫”
将文章分享到朋友圈。
阅读(5878)|
用微信&&“扫一扫”
将文章分享到朋友圈。
用易信&&“扫一扫”
将文章分享到朋友圈。
历史上的今天
在LOFTER的更多文章
loftPermalink:'',
id:'fks_081067',
blogTitle:'MATLAB中的点乘和点除问题',
blogAbstract:'\r\n&\r\n如：A=[1 2 3;4 5 6;7 8 9]\r\n&&&&&& B=[2 4 6;8 10 12;14 16 18]\r\nA*B= 60&&& 72&&& 84&&&&& &132&& 162&& 192&&&&& &204&& 252&& 300\r\nB/A=\r\nWarning: Matrix is singular to working precision. \r\nans =\r\n&& NaN&& NaN&& NaN&& NaN&& NaN&& NaN&& NaN&& NaN&& NaN\r\n%矩阵这样做是不行的！！！',
blogTag:'',
blogUrl:'blog/static/',
isPublished:1,
istop:false,
modifyTime:2,
publishTime:0,
permalink:'blog/static/',
commentCount:0,
mainCommentCount:0,
recommendCount:1,
bsrk:-100,
publisherId:0,
recomBlogHome:false,
currentRecomBlog:false,
attachmentsFileIds:[],
groupInfo:{},
friendstatus:'none',
followstatus:'unFollow',
pubSucc:'',
visitorProvince:'',
visitorCity:'',
visitorNewUser:false,
postAddInfo:{},
mset:'000',
remindgoodnightblog:false,
isBlackVisitor:false,
isShowYodaoAd:false,
hostIntro:'乐于交友,乐于交往够意思的朋友,乐于交往所有乐于交往的朋友。',
hmcon:'1',
selfRecomBlogCount:'0',
lofter_single:''
{list a as x}
{if x.moveFrom=='wap'}
{elseif x.moveFrom=='iphone'}
{elseif x.moveFrom=='android'}
{elseif x.moveFrom=='mobile'}
${a.selfIntro|escape}{if great260}${suplement}{/if}
{list a as x}
推荐过这篇日志的人：
{list a as x}
{if !!b&&b.length>0}
他们还推荐了：
{list b as y}
转载记录：
{list d as x}
{list a as x}
{list a as x}
{list a as x}
{list a as x}
{if x_index>4}{break}{/if}
${fn2(x.publishTime,'yyyy-MM-dd HH:mm:ss')}
{list a as x}
{if !!(blogDetail.preBlogPermalink)}
{if !!(blogDetail.nextBlogPermalink)}
{list a as x}
{if defined('newslist')&&newslist.length>0}
{list newslist as x}
{if x_index>7}{break}{/if}
{list a as x}
{var first_option =}
{list x.voteDetailList as voteToOption}
{if voteToOption==1}
{if first_option==false},{/if}&&“${b[voteToOption_index]}”&&
{if (x.role!="-1") },“我是${c[x.role]}”&&{/if}
&&&&&&&&${fn1(x.voteTime)}
{if x.userName==''}{/if}
网易公司版权所有&&
{list x.l as y}
{if defined('wl')}
{list wl as x}{/list}小木虫 --- 600万学术达人喜爱的学术科研平台
&&查看话题
matlab程序问题求解答
function=bisect(f,a,b,delta)
% Input&&- f is the function input as a string 'f'
%& && &&&- a and b are the left and right endpoints
%& && &&&- delta is the tolerance
% Output - c is the zero
%& && &&&- yc = f(c)
%& && &&&- err is the error estimate for c
ya = feval(f,a);
yb = feval(f,b);
if ya * yb & 0, end
max1 = 1+ round((log(b-a)-log(delta))/log(2));
for k =1:max1
& & c =(a+b)/2;
& & yc = feval(f,c);
& & if yc==0;
& && &&&a=c;
& && &&&b=c;
& & elseif yb * yc & 0
& && &&&b=c;
& && &&&yb=
& && &&&a=c;
& && &&&ya=
& & if b - a & delta, break, end
c=(a+b)/2;
err=abs(b-a);
yc=feval(f,c);& &&&% 求函数在c处的函数值
for T=200:10:300
& &K = exp(-9./T+0.835931*log(T)...
& && & +2.*T-1.*T^2+5.*T^3)
& &f=inline('(4*K-1)*x^2-8*K*x+4*K');
& &bisect(f,0,1,10^-5)
主程序调用函数出错
展开阅读全文
学术科研必备，90%的学术科研者都在使用
关于matlab程序问题求解答的相关话题在小木虫APP已经有273位虫友给出了详细回复。
赶快查看回复吧！
给了一个与温度相关的平衡方程，求200-300每隔10的平衡常数，以及平衡转化率。。。要求用二分法求解~
主要是程序出现的问题是什么呢
& &26.5468
Error using inline/feval (line 23)
Not enough inputs to inline function.
Error in bisect (line 8)
ya = feval(f,a);
Error in main (line 5)
& &bisect(f,0,1,10^-5)
就是我声明的函数有K，被认为是变量，所以就少参数了。。。但是前面的K是可以求出来的~这个就不知道怎么弄了，它为什么不代进去啊~我刚开始的时候也像你一样迷茫，同学给我推荐了，上面很多牛人分享的科研经验，对新手很有帮助，你也可以下载来试试！
多谢~能顺便解释下为什么嘛~
关于matlab程序问题求解答的相关话题在小木虫APP已经有273位虫友给出了详细回复。
赶快查看回复吧！
学术必备与600万学术达人在线互动！
扫描下载送金币
地址：北京市海淀区北三环甲18号中鼎大厦A座1层102室 电话：010-
浏览器进程
打开微信扫一扫
随时随地聊科研[转载]matlab中遇到的一些问题
已有 14452 次阅读
|个人分类:|系统分类:|文章来源:转载
1什么是Gaussian拟合？:#FangQ(Qianqian.Fang@Dartmouth.Edu),, BigGreen/MathTools #&&&&&&&& 高斯拟合(Gaussian Fitting)即使用形如：&&&&&&&&&&&&&&& Gi(x)=Ai*exp((x-Bi)^2/Ci^2)&&&&&&&& 的高斯函数对数据点集进行函数逼近的拟合方法。&&&&&&&& 其实可以跟多项式拟合类比起来，不同的是多项式拟合是用幂函数系，&&&&&&&& 而高斯拟合是用高斯函数系。&&&&&&&& 使用高斯函数来进行拟合，优点在于计算积分十分简单快捷。这一点&&&&&&&& 在很多领域都有应用，特别是计算化学。著名的化学软件Gaussian98&&&&&&&& 就是建立在高斯基函数拟合的数学基础上的。2如何在给定句柄的axis里绘图?:#FangQ(Qianqian.Fang@Dartmouth.Edu),, SMTH/MathTools #&&&&&&&& plot(data,'parent',haxis);&&&&&&&& 或者&&&&&&&& hbar=bar(data);&&&&&&&& set(hbar,'parent',haxis);3)由Matlab符号运算得到的公式怎么才能将数据代进去运算？:#ramjet (德芙), SMTH/MathTools #&&&&&&&& 使用subs(),或先将值赋予一个符号变量，然后用eval()5)在Matlab中如何求最值点？如何求一维数组的极值？:#FangQ(Qianqian.Fang@Dartmouth.Edu),, SMTH/MathTools# &&&&&&&& 最值：&&&&&&&& 一维或多维数组最值用max(data(:))&&&&&&&& 如果想返回最值所在的位置，用[Y,I]=max(data):#FangQ(Qianqian.Fang@Dartmouth.Edu), ,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(diff(data)))==2)+1)&&&&&&&& 返回的是极大值和极小值6)Matlab中如何作线性拟合/线性回归/多元线性回归？:#FangQ(Qianqian.Fang@Dartmouth.Edu),, BigGreen/MathTools #&&&&&&&& 即用y=a*x+b来拟合一组数据{{x1,y1},{x2,y2}&{xn,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=1~n)&&&&&&&&&& |x11,x21,&xm1|&&&&&&&& A=|x12,x22,&xm2|&&&&&&&&&& |&&&&&&&& |&&&&&&&&&& |x1n,x2n,&xmn|&&&&&&&& Y={y1,y2,y3,&,yn}'&&&&&&&& 则系数{a1,a2,&,am}'=pinv(A)*Y&&&&&&&& 在matlab中使用&&&&&&&& coeff=A\Y&&&&&&&& 则可以得到最小二乘意义上的拟合系数7)Matlab中如何作圆回归？:#Peter Boettcher (boettcher@ll.mit.edu),, 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-squares 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 radius R.A is an &&&&&&&&&&&&&&&&& % optional output describing the circle's equation:&&&&&&&&&&&&&&&&& %&&&&&&&&&&&&&&&&& %&&& x^2+y^2+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)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=A\B;&&&&&&&&&&&&&&&&& xc = -.5*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.Fang@Dartmouth.Edu),, SMTH/MathTools #&&&&&&&& 到/matlabcentral/fileexchange/index.jsp&&&&&&&& 2-D Plotting and Graphics中查找arrow.m，或者&&&&&&&& /matlabcentral/spotlight/arrows.shtml&&&&&&&& http://www.math.umd.edu/~jec/matcomp/matcompmfiles/mfiles.html9)Matlab中如何作二维数据的插值?:#FangQ(Qianqian.Fang@Dartmouth.Edu),, BigGreen/MathTools #&&&&&&&& 对于一维、二维、三维规则数据点阵使用interp1/interp2/interp3，&&&&&&&& 二维、三维非规则数据用griddata/griddata310)Matlab中如何绘制三维数据阵？:#FangQ(Qianqian.Fang@Dartmouth.Edu),, BigGreen/MathTools #&&&&&&&& 如果使用matlab，打开帮助窗口，在目录树上找到&&&&&&&& MATLAB\Using 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,, SMTH/MathTools #&&&&&&&& 注释大段代码选中代码，Ctrl+R；取消注释，选中代码，Ctrl+T。&&&&&&&& 或者用Edit菜单或者右键弹出中的注释。:#misc,, SMTH/MathTools #&&&&&&&& if(0)&&&&&&&& 大段的代码&&&&&&&& end12)Matlab中如何计算程序运行的时间？:#misc,, SMTH/MathTools #&&&&&&&& tic&&&&&&&&&& your_&&&&&&&& toc&&&&&&&& 或者使用&&&&&&&& t= &&&&&&&&&& your_ &&&&&&&& cputime-t13)Matlab中如何改变默认的工作路径？:#SindyGong, , SMTH/MathTools #&&&&&&&& 编辑一个startup.m文件，其中cd yourpath&&&&&&&& 或者在X:\matlab\toolbox\local\matlabrc.m的最后添加cd yourpath&&&&&&&& 参见：&&&&&&&& /support/solutions/data/25164.shtml14)Matlab如何改变默认的图形字体?:#comp.soft-sys.matlab FAQ#&&&&&&&& 编辑一个startup.m文件，其中&&&&&&&& set(0,'DefaultObjectnamePropertyName',Value)&&&&&&&& 或者在X:\matlab\toolbox\local\matlabrc.m的最后添加&&&&&&&& set(0,'DefaultObjectnamePropertyName',Value)15)如何在Matlab中实现交互操作？:#FangQ(Qianqian.Fang@Dartmouth.Edu),,BigGreen/MathTools #&&&&&&&& 如果只在命令窗口进行交互操作，请参见demo中的例子，主要是&&&&&&&& 通过input命令和pause/clear/disp等实现的，还有一些窗口资源可以使&&&&&&&& 用：&&&&&&&& uigetfile,uiputfile,uiwait,uisetcolor,uisetfont, uiopen,uisave&&&&&&&& inputdlg,msgbox,helpdlg,questdlg,warndlg,errordlg16)Matlab中为什么只能在小数点后显示四位?:#FangQ(Qianqian.Fang@Dartmouth.Edu),,BigGreen/MathTools #&&&&&&&& 用format命令来改变命令窗口数字的显示格式和精度，但不会影&&&&&&&& 响matlab的计算精度，matlab的矩阵运算默认都是双精度浮点型运算。17)Matlab如何在命令窗口按照格式输出？:#FangQ(Qianqian.Fang@Dartmouth.Edu),,SMTHTools #&&&&&&&& fprintf(1,&your_format_string&,var1,var2,&);18)如何在Matlab中画隐函数曲线？:#FangQ(Qianqian.Fang@Dartmouth.Edu),,BigGreen/MathTools #&&&&&&&& 在/matlabcentral/fileexchange/index.jsp&&&&&&&& 查找implicit，会找到一个Arthur Jutan写的implot.m&&&&&&&& Mathematica中绘制隐函数用ImplicitPlot[]&&&&&&&& 或者ImplicitPlot3D[]&&&&&&&& Maple中为implicitplot(),implicitplot3d()&&&&&&&& 参见&&&&&&&& http://engineering.dartmouth.edu/~fangq/MATH/download/source/&&&&&&&&& ImplicitPlot3D.htm19)Matlab中什么函数可以删除矩阵的某一行或列？:#FangQ(Qianqian.Fang@Dartmouth.Edu),,BigGreen/MathTools #&&&&&&&& A(j,:)=[];&& %删除A的第j行&&&&&&&& A(:,i)=[];&& %删除A的第i列20)Matlab中能开的最大数组是由什么决定的？:# chenft (mike),, SMTH/MathTools #&&&&&& I have had similar problems. Below is an explanation I received from&&&&&&&& Ian Boyd &&&&&& from Mathworks (just giving credit where credit is due) that explains &&&&&& what's happening. You solution is to run matlab with the -nojvm mode.&&&&&& &The heap memory system in JAVA consists of data and handle elements. &&&&&& When you allocate a variable you get a handle and data. As long as&&&&&&&& data has an associated handle, the JVM considers it valid and &&&&&& will not clean it up.&&&&&& However, when you call the clear function in MATLAB, all handles are &&&&&& destroyed, and the data associated is now invalid. This means that&&&&&&&& the JAVA engine can free up that data (garbage collection), but does &&&&&& not mean that it will clean it up at that moment.&&&&&& Calling the PACK command encourages JAVA to run the garbage collector &&&&&& and de-fragment the memory. But it does not force it to (This is part&&&&&&&& of the JAVA design). Even though the memory is 'freed' on the heap,&&&&&&&& it is not actually free to the OS, it is only free to the JVM. Here&&&&&&&& is one way to think of it:&&&&&& [MATLAB]&&&&&& [JAVA]&&&&&&&& [OS]&&&&&& MATLAB runs on JAVA (virtual machine), and Java runs on the OS &&&&&& (physical machine). So when MATLAB is running in JAVA mode memory &&&&&& allocations&& are requested from the JRE, not the OS.&&&&&& One problem you may be running into is that the default maximum &&&&&& JAVA heap size is relatively low ( &= 64 MB), so that is all the &&&&&& memory one session of MATLAB will ever get on your system.&&&&&& The good news is that you can increase this value. You will need &&&&&& to create a java.opts file in $MATLAB/bin/$ARCH (or in the current &&&&&& directory when&& you start MATLAB) and put the following command:&&&&&& %%%BEGIN CODE%%%&&&&&& maxHeapSize = &&&&&& %%%END CODE%%%&&&&&& This will give you 256MB of JVM memory and you can adjust the &&&&&& parameter as needed.&&&&&& Note: $MATLAB is the root directory and $ARCH is your system&&&&&& architecture. This solution works on Windows as well as Solaris, &&&&&& Linux,Alpha, and SGI. A similar operation is possible on IBM and &&&&&& HPUX, but with a different syntax.&&&&&& For the 1.1.8 JVM (Windows, Linux, Solaris, Alpha, SGI) our &&&&&& defaults are:&&&&&&&&&&&&&&& minHeapSize = &&&&&&&&& maxHeapSize = &&&&&& These are the structure field names in&& that correspond to &&&&&& -ms and -mx, and the settings above are roughly 16MB and 64MB.&&&&&& To investigate the Java heap a bit, ask via the following:&&&&&&&&& && java.lang.Runtime.getRuntime.totalMemory&&&&&&&&& && java.lang.Runtime.getRuntime.freeMemory&&&&&& When the free memory hits zero, Java will double the heap size &&&&&& (up to the maximum setting).&&&&&& If you choose to run without Java, you will remove the overhead &&&&&& of the middle man, but you will also lose some MATLAB functionality &&&&&& (mostly graphics and the Editor). You will still have most &&&&&& of the computational power though.&&&&&& Without JAVA, memory management will come directly from the OS, &&&&&& and a CLEAR operation will result in memory being freed back to &&&&&& the OS.
21)如何在Matlab中添加新的工具箱？:#FangQ(Qianqian.Fang@Dartmouth.Edu),,BigGreen/MathTools #&&&&&&&& 如果是Matlab安装光盘上的工具箱，重新执行安装程序，选中即可。&&&&&&&& 如果是单独下载的工具箱，一般情况下仅需要把新的工具箱解压到某&&&&&&&& 个目录，然后用addpath(对于多个目录的使用genpath())或者pathtool添&&&&&&&& 加工具箱的路径，然后用which newtoolbox_command.m来检验是否可&&&&&&&& 以访问。如果能够显示新设置的路径，则表明该工具箱可以使用了。&&&&&&&& 具体请看工具箱自己代的README文件。22)如何读写Matlab的.mat文件？:#FangQ(Qianqian.Fang@Dartmouth.Edu),,BigGreen/MathTools #&&&&&&&& 文件结构参见：&&&&&&&& http://www.mathworks.de/access/helpdesk/help/pdf_doc/matlab/&&&&&&&&&& matfile_format.pdf&&&&&&&& /support/solutions/data/8757.shtml&&&&&&&& /pub/tech-support/solutions/s8757/&&&&&&&&&& readmemat.txt&&&&&&&& 建议使用matlab自己提供的函数来读写简单安全，或者参考：&&&&&&&& http://engineering.dartmouth.edu/~fangq/MATH/download/&&&&&&&&&&& source/mat_file.txt&&&&&&&& 来自matlab的c math library23)如何得到contour线上的坐标点？:#FangQ(Qianqian.Fang@Dartmouth.Edu),,BigGreen/MathTools #&&&&&&&& lcount=5;&&&&&&&& [c,h]=contour(peaks,lcount);&&&&&&&& x=get(h,'xdata');&&&&&&&& y=get(h,'ydata');&&&&&&&& 这里得到的x和y都是cell数组，用x{1}/y{1}来得到每条线上的坐标对，&&&&&&&& 注意，每条线的最后一个数据是NaN24)如何将Matlab绘制的三维网格图帖到word里？:#FangQ(Qianqian.Fang@Dartmouth.Edu),,BigGreen/MathTools #&&&&&&&& 如果需要位图，好处是所见即所得，坏处是图像精度差，不能放缩：&&&&&&&&&&&&&&&& 1.用拷屏 Alt+PrintScreen&&&&&&&&&&&&&&&& 2.在图形窗口菜单Edit\Copy Options&.\选择Bitmap，可以&&&&&&&&&&&&&&&&&& 选择透明背景，然后Edit\Copy Figure&&&&&&&& 如果需要拷贝矢量图：&&&&&&&&&&&&&&&& 在图形窗口菜单Edit\Copy Options&.\选择Metafile,然后&&&&&&&&&&&&&&&& Edit\Copy Figure，在Word中粘贴&&&&&&&& 经常地，按照Metafile方式粘贴的图片曲线会出现锯齿，最好的方式是&&&&&&&& 使用eps文件：&&&&&&&&&&&&&&&& 1.将需要拷贝的图作为当前窗口&&&&&&&&&&&&&&&& 2.再转换到matlab命令窗口，print -deps filename.eps&&&&&&&&&&&&&&&& 3.-deps还可以用depsc,deps2,depsc2&&&&&&&&&&&&&&&& 4.在word中插入图片，选中该eps，如果是word 2000以前版本&&&&&&&&&&&&&&&&&& ，不会显示图片内容，但可以打印，word XP即可显示，又可打印。&&&&&&&&&&&&&&&& 5.如果不满意，可以在word中双击编辑，如果安装有Adobe &&&&&&&&&&&&&&&&&& Illustrator等矢量图像编辑软件，也可以进行编辑。25)请问可以查看Matlab中函数的源代码吗?:#FangQ(Qianqian.Fang@Dartmouth.Edu),,BigGreen/MathTools #&&&&&&&& Matlab除了buildin函数和mex/dll文件看不到原码，其他如工具箱等都可&&&&&&&& 以直接看到代码，首先确认该文件安装在matlab中，即which &&&&&&&& filename.m存在，然后可以edit filename.m26).Matlab有没有求矩阵行数/列数/维数的函数？:#FangQ(Qianqian.Fang@dartmouth.edu),,BigGreen/MathTools #&&&&&&&&& ndims(A)返回A的维数&&&&&&&&& size(A)返回A各个维的最大元素个数&&&&&&&&& length(A)返回max(size(A))&&&&&&&&& [m,n]=size(A)如果A是二维数组，返回行数和列数&&&&&&&&& nnz(A)返回A中非0元素的个数27).Matlab中如何中断运算？:#FangQ(Qianqian.Fang@dartmouth.edu),,BigGreen/MathTools #&&&&&&&&& 在命令窗口按Ctrl+C，在UNIX/LINUX会立即中断运算，在Windows可&&&&&&&&& 能由于操作系统的原因，有时会出现死机和等待的情况。28).Matlab中有没有画圆或椭圆的函数？:#FangQ(Qianqian.Fang@dartmouth.edu),,BigGreen/MathTools #&&&&&&&&& 没有，Matlab没有提供直接绘圆的图元函数，需要自己写代码，其实&&&&&&&&& 就两句：&&&&&&&&& sita=0:pi/20:2*&&&&&&&&& plot(r*cos(sita),r*sin(sita)); %半径为r的圆&&&&&&&&& plot(a*cos(sita+fi),b *sin(sita+fi)); %椭圆&&&&&&&&& 如果是单位圆，可以使用rectangle('Curvature', [1 1])29).Matlab下如何定义整形:#修改：fhorse (马不停蹄),,SMTH/MathTools #&&&&&&&&& Matlab默认的矩阵数据结构都是双精度浮点型，即64位来表示一个数&&&&&&&&& 字，大多数的函数和操作都定义在double数据结构，如果你需要&&&&&&&&& 把double的数据转换为整形，然后再参与运算，需要使用&&&&&&&&& double(int32(x))或者floor/round/ceil等函数&&&&&&&&& 如果为了节省内存，只进行赋值、打印等简单操作，可以参&&&&&&&&& 见uint8/uint16/uint32命令的帮助30).Matlab如何产生均匀分布的白噪声? :#misc,,SMTH/MathTools #&&&&&&&&& help rand&&& 均匀分布百噪声&&&&&&&&& help randn高斯分布百噪声31).在Matlab中debug的时候能否跟踪变量的?:#FangQ(Qianqian.Fang@dartmouth.edu),,BigGreen/MathTools #&&&&&&&&& 可以，如果使用medit，设置断点后可以用鼠标移到所看的变量上，显&&&&&&&&& 示当前的值，或者在命令窗口打该变量名直接回车。如果在代码中实&&&&&&&&& 现调试断点等功能，参&&&&&&&&& 见dbstop,dbcont,dbstep,dbclear,dbtype,dbstack,dbup,dbdown,dbstatus,&&&&&&&&& dbquit32).请问在Matlab中怎样输入特殊符号啊或者上标、下标？:#FangQ(Qianqian.Fang@dartmouth.edu),southerner(笑着),,SMTH/MathTools#&&&&&&&&& matlab的text/title/xlabel/ylabel对象支持简单的TeX排版语法，如希腊字&&&&&&&&& 母，上下标等例如&&&&&&&&& text(0.5,0.5,'\alpha^\beta_2');33).Matlab中如何后台运行一个DOS程序?:#FangQ(Qianqian.Fang@dartmouth.edu), . BigGreen/en_Matlab#&&&&&&&&& 这里是一个后台执行一个需要外部输入的DOS命令的例子,需要的输&&&&&&&&& 入实事先都写在同目录下的input.txt文件中：&&&&&&&&& dos('myexe & input.txt &')34).Matlab如何加载输入文件(批处理模式) ？:#翻译自：comp.sys-soft.Matlab FAQ. BigGreen/en_Matlab#&&&&&&&&& PC上可以使用matlab /r参数来在matlab启动的时候直接加载运行m文件&&&&&&&&& ，在UNIX上，使用&&&&&&&&& matlab & MyMFile & MyOutputFile&&&&&&&&& 来外部执行MyMFile，&&&&&&&&& 以上执行方式都可以通过脚本文件实现批处理35).Matlab如何启动时执行规定的文件？:#FangQ(Qianqian.Fang@dartmouth.edu), .BigGreen/en_Matlab#&&&&&&&&&&&&&&&&& 参见上一个问题的回答36).如何在Matlab GUI中使用图形背景？:#FangQ(Qianqian.Fang@dartmouth.edu), .BigGreen/en_Matlab#&&&&&&&&& 这是一个简单的例子：&&&&&&&&& [A,map]=imread('yourimg.gif');&&&&&&&&& imagesc(A)&&&&&&&&& colormap(map)&&&&&&&&& set(gca,'position',[0 0 1 1])&&&&&&&&& axis off&&&&&&&&& ax2=axes('position',[0.2,0.2,0.6,0.6]);&&&&&&&&& plot(rand(1,10),'parent',ax2);&&&&&&&&&&&& set(ax2,'color','none')&&&&&&&37).大量数据点Matlab绘图为什么很慢？:#FangQ(Qianqian.Fang@dartmouth.edu), .BigGreen/en_Matlab#&&&&&&&&& 1.首先看能否用已有函数对整个矩阵绘图，比&&&&&&&&&&& 如mesh/plot3/trimesh等&&&&&&&&& 2.如果必须一点一点/或者一条线一条线的添加，最好作如下&&&&&&&&&&& 设置：&&&&&&&&&&& doublebuffer=on&&&&&&&&&&& erasemode=none&&&&&&&&&&& backingstore=off&&&&&&&&&&& renderer＝opengl&&&&&&&&&&& 以及参考MathWorks对于高速绘图的tips:&&&&&&&&&&& /support/tech-notes/v5/.shtml,38).Matlab中如何求解广义积分？即积分限到有无穷的或者有奇异点的积分(瑕积分)？:#FangQ(Qianqian.Fang@dartmouth.edu), . BigGreen/en_Matlab#&&&&&&&&& Matlab的quad/quad8只能作定积分，广义积分需要自己来写程序逼近，&&&&&&&&& 流程大概如下：&&&&&&&&&&&&&&&&&&&&&&&&& 1.设定收敛限epsi&&&&&&&&&&&&&&&&& 2.把为inf/-inf或者歧义点的积分限设置为一个初始值,k=1&&&&&&&&&&&&&&&&& 3.计算定积分Q(k)&&&&&&&&&&&&&&&&& 4.然后朝着inf/-inf或者歧义点移动一个步长，然后计算定积分&&&&&&&&&&&&&&&&&&& Q(k+1)&&&&&&&&&&&&&&&&& 5.判断(abs(Q(k+1)-Q(k)) &&&&&&&&&&&&&&&&&&&&&&&&& Mathematica中可以使用NIntegrate[]，对于无穷振荡的函数，可以使用Method-&&&&&&&&&& QuasiMonteCarlo或者Oscillatory]39).为什么我的Matlab程序这么慢？:#FangQ(Qianqian.Fang@dartmouth.edu), . BigGreen/en_Matlab#&&&&&&&&& 我们工学院的收发室的门上贴着一张小纸条，写的是&&&&&&&&& &Our policy is always blaming the computer&&&&&&&&&& 大多数的人在遇到问题的时候，总是责备计算机如何如何，别人如何&&&&&&&&& 如何，其实，最最主要的因素是在于自己。&&&&&&&&& 一个程序运行快慢，有很多因素决定，最主要的是算法，简炼而优美&&&&&&&&& 的的数学公式胜过100遍的优化。能从算法上改进，才能比别人有根&&&&&&&&& 本的优势。计算机也很重要，以前我总把自己用的PC看成万能的加&&&&&&&&& 以崇拜，对UNIX嗤之以鼻，结果当自己真正开始算起来，才知道差&&&&&&&&& 别有多大。搞大型数值计算的，没有好的工作站或者并行系统，就输&&&&&&&&& 在了起跑线上了。然后是程序的优化，看看变量是否占用太多内存，&&&&&&&&& 看看是否有功能重复的模块或者计算，经常的是用牺牲内存来换取速&&&&&&&&& 度，具体取舍，具体需要来决定。用profile看看哪些语句占用时间最&&&&&&&&& 多，然后把核心部分进行优化。&&&&&&&&& 如果是使用Matlab，使用vectorization和矩阵整体操作的代码要比大量&&&&&&&&& 的for循环快很多，eval/inline函数如果出现在核心循环，也会让速度下&&&&&&&&& 降几时倍的。40)..Matlab中如何作非线性回归？:#FangQ(Qianqian.Fang@dartmouth.edu), . BigGreen/en_Matlab#&&&&&&&&& 请参考&&&&&&&&& /support/solutions/data/10652.shtml&&&&&&&&&&&&&&&&& matlab默认只提供了多项式拟合的函数polyfit，对于其他稍微简单&&&&&&&&& 一点的拟合，如标准的指数、对数、高阶多项式拟合，都有解析公式，参见：&&&&&&&&& /LeastSquaresFitting.html&&&&&&&&& 对于更加复杂的非线性函数，建议使用Mathematica或者DataFit&&&&&&&&&&&&&&&&&&&&&&&&& Mathematica中提供了Fit[],以及&&&&&&&&& && Statistics`NonlinearFit`&&&&&&&&& NonlinearFit[],NonlinearRegress[]&&&&&&&&& 可以拟合任意复杂的表达式。&&&&&&&&&&&&&&&&&&&&&&&&& DataFit可以自定义拟合模型，适用于复杂系统的拟合。41)..Matlab中为什么我对m文件、simulink模块，mat文件的修改不起作用呢？:hyphone,,SMTH/MathTools#&&&&&&&&& 检查Matlab路径中是否有与你的m文件、mdl文件或者mat文件同名的m文件、&&&&&&&&& mdl文件或者mat文件。Matlab执行搜索到的第一个文件。
42).Matlab中，函数里面怎样使用基本工作空间中的变量？为什么inline函数不能使用外面的变量？:hyphone,,SMTH/MathTools#&&&&&&&&& 函数只能存取它自己的工作空间中的变量。要在函数之间，或者函数与基本&&&&&&&&& 工作空间之间传递数值，尝试以下方法：&&&&&&&&& 1、使用全局变量，用global定义全局变量&&&&&&&&& 2、使用evalin：&&&&&&&&&&&&&&&&& evalin('base','v=1;');&&& %在基本工作空间中执行命令：v=1;&&&&&&&&&&&&&&&&& evalin('caller','v=1;');%在调用该函数的函数的工作空间中&&&&&&&&&&&&&&&&&&&&&&&&& 执行命令：v=1;&&&&&&&&& 3、参数传递&&&&&&&&& 4、使用assignin：&&&&&&&&&&&&&&&&& assignin('base','v',v); %将v赋给基本工作空间中的变量v;&&&&&&&&&&&&&&&&& assignin('caller','v',v);&&& %将v赋给调用者工作空间中的变量v;&&&&&&&&& inline函数里面只能出现函数和参数，要传递一个可变系数，需要用evalin。&&&&&&&&& 如：fzero('f(evalin(''base'',''x1(i)''),x3)',求解区间）43).怎样在Simulink中调用m文件？:hyphone,,SMTH/MathTools#&&&&&&&&& 函数m文件可以使用Function & Tables中的Fcn模块。如果有多个输入，&&&&&&&&& 用Mux组合成一个向量，然后在Fcn模块的Expression填&&&&&&&&& MyFunction(u(1),u(2))&。如果有多个输出，用Demux分解成多个标量。44).Matlab中怎样进行数制转换？:hyphone,,SMTH/MathTools#&&&&&&&&&&&&&&&&& 参见：hex2dec,oct2dec,bin2dec,dec2bin,dec2oct,dec2hex.45).matlab中的*.p是什么文件？怎么用？ :#David Martin(dg.martin@), ,comp.soft-sys.matlab#&&&&&&&&& p文件是pre-parsed的缩写，即matlab在第一次运行某m文件时，matlab把该&&&&&&&&& 文件先编译成一种matlab的pseudo-code，当你再次运行该m文件时，节省了&&&&&&&&& parse的时间。&&&&&&&&& p文件可以离开m文件单独运行:#Mike Robbins (michael.robbins@),,comp.soft-sys.matlab#&&&&&&&&& 在matlab中用pcode('yourfilename')来生成p文件46).在Matlab中有goto语句吗？ :#FangQ(Qianqian.Fang@dartmouth.edu), . BigGreen/en_Matlab#&&&&&&&&& matlab中没有提供goto，因为结构化程序设计不推荐使用goto,但在matlab&&&&&&&&& 中，goto的功能可以部分用结构化的异常处理机制来实现，比如从多重循环&&&&&&&&& 中跳出，具体代码为：&&&&&&&&& try&&&&&&&&&&&& for i=1:10&&&&&&&&&&&&&& for j=1:10&&&&&&&&&&&&&&&&&&&& for k=1:10&&&&&&&&&&&&&&&&&&&&&&& do_&&&&&&&&&&&&&&&&&&&&&&&&&&& if(jump_condition)&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& errorid=-2;&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& error('I want to get out!');&&&&&&&&&&&&&&&&&&&&&&&&&&& end&&&&&&&&&&&&&&&&&&&&& end&&&&&&&&&&&&&&&&&&& end&&&&&&&&&&&&&&&&& end&&&&&&&&& catch&&&&&&&&&&&&& fprintf(1,'catch error:%d',errorid);&&&&&&&&& end47).请问matlab6.X的那个matlab server是做什么的？(matlab开机运行问题):#hairhead, . SMTH/MathTools#&&&&&&& 如果你在安装matlab时选择了web server的话，在每次重新启动后系统&&&&&&& 进程中就会出现matlab/matlabserver两个进程，如果你不需要使用&&&&&&& matlab webserver服务的话，建议在安装时不要选择web server，&&&&&&& 如果已经安装的话，可以在&控制面板-&管理工具-&服务&中找到&&&&&&& matlab webserver，然后把它disable掉。48).Matlab中如何用鼠标取得坐标？:#FangQ(Qianqian.Fang@dartmouth.edu), . BigGreen/en_Matlab#&&&&&&& matlab中用ginput来取得图像上的一点。注意：image()/imagesc()对象&&&&&&& 的y轴是和一般图的y轴反的。&&&&&&& 如果复杂的话，可以在ButtonDownFcn中get(0,'PointerLocation')&&&&&&& 或者get(gcf,'Position');49).Matlab中有阶乘函数吗？:#FangQ(Qianqian.Fang@dartmouth.edu), . BigGreen/en_Matlab#&&&&&&&&&&& matlab没有提供直接的阶乘函数，但可以用prod(1:n)来求n!&&&&&&& 用prod(1:2:2n-1)或者prod(2:2:2n)来求解n!!51)..Matlab如何计算大阶乘？:#lll, . SMTH/MathTools#&&&&&&&&&&& 如果只需要大致的值，取log10,&&&&&&& 计算出result=log10(1000!)=log10(1)+log10(2)+...log10(1000)&&&&&&& 然后求10^result=10^result的小数部分*10^result的整数部分52)..Matlab中怎样求变上限二重积分？:#hyphone,, SMTH/MathTools #&&&&&&&&& 对于解析函数，用两次int即可。&&&&&&&&& 如求x+y在0[返回]53)..用符号积分算出来Ei是什么意思，怎样求值？:#hyphone,, SMTH/MathTools #&&&&&&&&& Ei是maple中的指数积分函数，表示exp(-x*t)/t^n对t从1到正无穷大的积分。&&&&&&&&& Ei(n,x) = int(exp(-x*t)/t^n, t=1..infinity)&&&&&&&&& 用命令&mhelp Ei&可以查看详细说明；用符号计算得到的不明白的函数都可以通过&&&&&&&&& mhelp命令得到帮助。&&&&&&&&& 求Ei(1,2)可以用maple命令得到：str2num(maple('evalf(Ei(1,2))'))。
54)..Maltab中使用\n换行在notepad中显示为小黑块，为什么？(Randy Poe):#Randy Poe,,comp.soft-sys.matlab#:#Peter J. Acklam(pjacklam@online.no),,comp.soft-sys.matlab#&&&&&&&&& 换行和回车是不同的，而且在不同的操作系统，解释也不相同。&&&&&&&&& \n一般会操作系统被翻译成&行的结束&，即LF(Line-Feed)&&&&&&&&& \r会被翻译成&回车&，即CR(Cariage-Return)&&&&&&&&& 对于文本文件的新行，在UNIX上，一般用\n(LF)来表示，Mac上用\r(CR)来表示，&&&&&&&&& Windows上是用\n\r(CR-LF)来表示。&&&&&&&&&&&&&&&&& 所以在matlab中使用\n来写回车，在windows上打开会出现小方块。如果想&&&&&&&&& 避免这种情况，打开文件时使用t参数：&&&&&&&&&&&&&&&&&&&&& fid = fopen('myfile.dat','wt');55)..Matlab中能开多大数组？(Steven Lord):#Steven Lord (),,comp.soft-sys.matlab#&&&&&&&&& 使用computer命令：&&&&&&&&& [C,MAXSIZE] = computer56)..如何使用整型矩阵来节省内存？(Duane Hanselman):#Duane Hanselman (),,comp.soft-sys.matlab#&&&&&&&&& 如下方法无需先生成一个double的数组，然后转换为int8&&&&&&&&&&&&&&&&& && rc=[3 4] % row and column sizes&&&&&&&&& && a(prod(rc))=int8(0) % example using int8&&&&&&&&& && class(a)&&& % they are all int8s&&&&&&&&& && reshape(a,rc) % make it the size you want&&&&&&&&&&&&&&&&& 同时也可以使用repmat来实现上述功能：&&&&&&&&&&&&&&&&&&&&&&& && repmat(int8(0),rc)57).Matlab在P4芯片上无法启动的解决方案:#FangQ(Qianqian.Fang@dartmouth),,SMTH/MathTools#&&&&&&&&& 如果不想安装补丁，只能使用matlab -nojvm的形式启动&&&&&&&&& matlab，否则需要参照如下页面的解决方案来安装补丁：&&&&&&&&& /groups?hl=en&lr=&ie=UTF-8&oe=UTF-8&selm=&&&&&&&&&&&&&&&& NQ8a9..newsread1.prod.itd.earthlink.net&&&&&&&&& /support/solutions/data/27293.shtml=================================== 58).如何求解对离散点的最优椭圆拟合？(Andrew Fitzgibbon, et al):#Authors: Andrew Fitzgibbon, Maurizio Pilu, Bob Fisher&Direct Least Squares Fitting of Ellipses&, IEEE T-PAMI, 1999#&&&&&&&&& http://bbs.dartmouth.edu/~fangq/MATH/Source/fitellipse.m59).Matlab/Mathematica中如何中断当前运算？(FangQ):#FangQ(Qianqian.Fang@dartmouth),,BigGreen/MathTools#&&&&&&&&& Mathematica中使用:&&&& Alt+./Alt+,&&&&&&&&& Matlab中使用:&&&&&&&&& Ctrl+C60)).Matlab/Mathematica中如何检查括号匹配？(FangQ):#FangQ(Qianqian.Fang@dartmouth.edu),,SMTH/MathTools#&&&&&&&&& Mathematica中使用:&&&& Ctrl+.&&&&&&&&& Matlab editor中使用: Ctrl+B61)..Matlab的GUI中为何无法使用uicontrol的句柄？(FangQ):#FangQ(Qianqian.Fang@dartmouth.edu),,SMTH/MathTools#&&&&&&&&& 如果你在执行GUI时出现&handles not defined&的错误时，&&&&&&&&& 打开GUIDE,把figure的HandleVisibility设置为on或者callback,&&&&&&&&& 则你可以在各个callback中直接使用handles，而不用声明global&&&&&&&&&&&&&&&&& 例如:&&&&&&&&&&&&&&&&& dat=get(handles.figure1,'userdata');&&&&&&&&& dat2=fliplr(dat);&&&&&&&&& se(handles.edit1,'userdata');62)..Matlab中如何把向量拓展成矩阵？(Zealous/FangQ):#Zealous/FangQ,,SMTH/MathTools#&&&&&&&&& 可以使用repmat()，例如：&&& repmat([1,2,3]',1,5)&&&&&&&&& 或者使用kron(),例如:&&&&&&& kron([1 2 3]',ones(1,5))63)..Matlab的GUI中的按钮如何在运行时移动及改变大小？(FangQ) :#FangQ(Qianqian.Fang@dartmouth),,refer to comp.soft-sys.matlab#&&&&&&&&& 使用selectmoveresize函数，例如：&&&&&&&&&&&&&&&&& figure&&&&&&&&& h=uicontrol('style','pushbutton');&&&&&&&&& set(h,'ButtonDownFcn',...&&&&&&&&&&&&& 'set(h,''selected'',''off'')',...&&&&&&&&&&&&& 'Enable','inactive')64)..Matlab如何求解维数巨大的稀疏矩阵方程？(FangQ):#FangQ(Qianqian.Fang@dartmouth),,BigGreen/MathTools#&&&&&&&&& Matlab提供了非常丰富的迭代型矩阵求解器，方法包括CG,BiCG,BiCGSTAB,&&&&&&&&& CGS,GMRES,LSQR,MINRES,P-CG,QMR,SYMMLQ等，为了加快矩阵求解速度，还&&&&&&&&& 提供了两个Preconditioner函数：luinc(不完全LU分解)，和cholinc(不完全&&&&&&&&& Cholesky分解)。下面是一个使用GMRES求解方程Amat*x=rhs的一个简单例子：&&&&&&&&&&&&&&&&& [L2,U2] = luinc(Amat,1e-3);&&&&&&&&& tic&&&&&&&&&&&& x=gmres(Amat,rhs,[],1e-6,100,L2,U2);&&&&&&&&& toc
在matlab中如何控制数据的输出格式？（完整总结）sunooy 发表于
format 缺省时为默认短格式方式与format short相同 format short 短格式方式，显示5位定点十进制数。 format long 长格式方式，显示15位定点十进制数。
format short g当数据大于1000或小于1时便会以科学记数法显示（-e），若想坚持用整数部分加小数部分的格式来显示，就要再后边加 g
format long gformat hex 十六进制格式方式。 format bank 银行格式。按元、角、分（小数点后具有两位）的固定格式。 format + +格式，以＋，&和空格分别表示中的正数，负数和零元素 format short&& e 短格式e方式，显示5位浮点十进制数 format long&&& e 长格式e方式，显示15位浮点十进制数。 format rat 分数格式形式。用有理数逼近显示数据。如pi显示为355/113。 format loose 松散格式。数据之间有空行。 format compact 紧凑格式。数据之间无空行。
vpa(date,n)将数据date以n位有效数字显示。
在图形窗口，依次进入 edit /aexs properties / （X，Y，Z，style，lengedn，选择你要改变字体的对象）/font size 选择你要的字体大小就可以了。
转载网址：
1我有一组x,y,z值，非规则点阵，如何绘制曲面或者等高线或者插值？:#FangQ(Qianqian.Fang@Dartmouth.Edu),, SMTH/MathTools #&&&&&&&& 1.如果你使用matlab，有如下四种解决办法&&&&&&&&&& 1.1: &&&&&&&&&&&&&&&&&&&&&&&&& tri=delaunay(x,y);&&&&&&&&&&&&&&&&&&&&&&&&& trimesh(tri,x,y,z); or trisurf(tri,x,y,z);&&&&&&&&&& 1.2: 如果你想得到更加精细和光滑的曲线，而且想绘制coutour&&&&&&&&&&&&&&&&&&& 则需要用griddata()插值：&&&&&&&&&&&&&&&&&&&&&&&&& [xi,yi]=meshgrid(min(x):dx:max(x),min(y):dy:max(y));&&&&&&&&&&&&&&&&&&&&&&&&& [zi,xi,yi]=griddata(x,y,z,xi,yi,'cubic');&&&&&&&&&&&&&&&&&&&&&&&&& contourf(xi,yi,zi,20);&&&&&&&&&& 1.3: 如果你的matlab安装有PDE Toolbox，可以使用pdeplot()，参见帮助&&&&&&&&&& 1.4: 如果你可以连接国外，可以使用MIT的Kirill K. Pankratov写的SaGA&&&&&&&&&&&&&&& 工具箱，下载地址在:&&&&&&&& http://www.ai-geostats.org/software/Geostats_software/SAGA.htm&&&&&&&& http://globec.whoi.edu/software/saga/&&&&&&&& 2.如果你使用Tecplot&&&&&&&&&& 1.先将三列数据按照顺序保存在一个文本文件中&&&&&&&&&&&&&&&& data.dat文件格式为&&&&&&&&&&&&&&&&& x1 y1 z1&&&&&&&&&&&&&&&&& x2 y2 z2&&&&&&&&&&&&&&&&& ...&&&&&&&&&& 2. Tecplot菜单File\Load Data file,在左上角选择2D&&&&&&&&&& 3. 选中zone,菜单Data\Triangulate&&&&&&&&&& 4. 双击该zone,将mesh掩藏掉，把contour打开，然后在&&&&&&&&&&&&&&&&& 左上角选择3D&&&&&&&&&& 5. 如果想插值，先创建一个更细的网格，然后用data菜单中&&&&&&&&&&&&&&&&&&&&& Interpolation&&&&&&&& 3.如果你使用Mathematica&&&&&&&&&& data={{x1,y1,z1},{x2,y2,z2},......};&&&&&&&&&& fun=Interpolation[data];&&&&&&&&&& Plot3D[fun[x,y],{x,xmin,xmax},{y,ymin,ymax}];&&&&&&&&&& ContourPlot[fun[x,y],{x,xmin,xmax},{y,ymin,ymax}];&&&&&&&& 当然，其他还有很多绘图工具，有待补充2请问哪里能找到绘制等高线的算法？:#FangQ(Qianqian.Fang@Dartmouth.Edu),, SMTH/MathTools #&&&&&&&&& 在一般的计算机图形学的教材中都有常见算法的描述，下面的网址&&&&&&&&&&&&&&&&& 有免费的C代码&&&&&&&&&&&&&&&&& http://members.bellatlantic.net/~vze2vrva/thesis.html3特征值与奇异值有什么区别？:#dimond,, SMTH/MathTools #&&&&&&&& 特征值就不用说了，大家都知道。&&&&&&&& 奇异值的定义：对于任一个复(实)矩阵A(m*n),设A^为A的共轭转置矩阵，&&&&&&&& 即A^是n*m阶,则 A^*A （n*n阶）的特征值的非负平方根为 A 的奇异值，&&&&&&&& 也就是A共有n个奇异值，且全部&=0.&&&&&&&& A^*A 是一个半正定矩阵，所以它的特征值&=0&&&&&&&& 对于半正定的Hermit方阵而言，特征值与奇异值是一一对应的
4条件数很大(病态)的矩阵如何求逆？:#granule,, SMTH/MathTools #&&&&&&&& 有许多解法！实际上这作为数学的一个分支已经数学家&&&&&&&& 们抽象出来研究了很多年(ill-posed problems).文献&&&&&&&& 多的不得了。&&&&&&&& 简单的解法就是用regularization tools.（典型的是&&&&&&&& Tikhonov regularization).还有许多迭代解法。&&&&&&&& 基本思想是要加约束。&&&&&&&& for Ax=b, its standard Tikhonov solution is&&&&&&&& x=inv(trans(A)*A+alfa*I)*Trans(A)*b&&&&&&&& trans() inv()&&&&&&&& alfa is a positive sclar. The quality of solution&&&&&&&& depends on alfa directly. There is a optimal alfa&&&&&&&& You can refer to the relevant topics such as:&&&&&&&& inverse problems, ill-conditioned problems.
转载本文请联系原作者获取授权，同时请注明本文来自孙亭科学网博客。链接地址：
上一篇：下一篇：
当前推荐数：1
评论 ( 个评论)
扫一扫，分享此博文
作者的其他最新博文
热门博文导读
Powered by
Copyright &