更多请见：http://bbs.dartmouth.edu/~fangq/MATH/FAQ/

40)..Matlab中如何作非线性回归？
:#FangQ(Qianqian.Fang@dartmouth.edu), 2002/6/22. BigGreen/en_Matlab#

请参考
http://www.mathworks.com/support/solutions/data/10652.shtml

matlab默认只提供了多项式拟合的函数polyfit，对于其他稍微简单
一点的拟合，如标准的指数、对数、高阶多项式拟合，都有解析公式，参见：
http://mathworld.wolfram.com/LeastSquaresFitting.html
对于更加复杂的非线性函数，建议使用Mathematica或者DataFit

Mathematica中提供了Fit[],以及
<< Statistics`NonlinearFit`
NonlinearFit[],NonlinearRegress[]
可以拟合任意复杂的表达式。

DataFit可以自定义拟合模型，适用于复杂系统的拟合。

=================================== - [返回]
41)..Matlab中为什么我对m文件、simulink模块，mat文件的修改不起作用呢？
:hyphone,2002/7/6,SMTH/MathTools#

检查Matlab路径中是否有与你的m文件、mdl文件或者mat文件同名的m文件、
mdl文件或者mat文件。Matlab执行搜索到的第一个文件。

=================================== - [返回]
42).Matlab中，函数里面怎样使用基本工作空间中的变量？
为什么inline函数不能使用外面的变量？
:hyphone,2002/7/6,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,2002/7/6,SMTH/MathTools#

函数m文件可以使用Function & Tables中的Fcn模块。如果有多个输入，
用Mux组合成一个向量，然后在Fcn模块的Expression填
MyFunction(u(1),u(2))"。如果有多个输出，用Demux分解成多个标量。

=================================== - [返回]
44).Matlab中怎样进行数制转换？
:hyphone,2002/7/6,SMTH/MathTools#

参见：hex2dec,oct2dec,bin2dec,dec2bin,dec2oct,dec2hex.

=================================== - [返回]
45).matlab中的*.p是什么文件？怎么用？
:#David Martin(dg.martin@home.deleteme.com), 2000/07/07,comp.soft-sys.matlab#

p文件是pre-parsed的缩写，即matlab在第一次运行某m文件时，matlab把该
文件先编译成一种matlab的pseudo-code，当你再次运行该m文件时，节省了
parse的时间。

p文件可以离开m文件单独运行

:#Mike Robbins (michael.robbins@us.cibc.com),2001/04/20,comp.soft-sys.matlab#

在matlab中用pcode('yourfilename')来生成p文件

=================================== - [返回]
46).在Matlab中有goto语句吗？
:#FangQ(Qianqian.Fang@dartmouth.edu), 2002/7/23. BigGreen/en_Matlab#

matlab中没有提供goto，因为结构化程序设计不推荐使用goto,但在matlab
中，goto的功能可以部分用结构化的异常处理机制来实现，比如从多重循环
中跳出，具体代码为：

try
for i=1:10
for j=1:10
for k=1:10
do_something;
if(jump_condition)
errorid=-2;
error('I want to get out!');
end
end
end
end
catch
fprintf(1,'catch error:%d',errorid);
end

=================================== - [返回]
47).请问matlab6.X的那个matlab server是做什么的？(matlab开机运行问题)
:#hairhead, 2002/4/3. SMTH/MathTools#

如果你在安装matlab时选择了web server的话，在每次重新启动后系统
进程中就会出现matlab/matlabserver两个进程，如果你不需要使用
matlab webserver服务的话，建议在安装时不要选择web server，
如果已经安装的话，可以在"控制面板->管理工具->服务"中找到
matlab webserver，然后把它disable掉。

=================================== - [返回]
48).Matlab中如何用鼠标取得坐标？
:#FangQ(Qianqian.Fang@dartmouth.edu), 2002/11/19. BigGreen/en_Matlab#

matlab中用ginput来取得图像上的一点。注意：image()/imagesc()对象
的y轴是和一般图的y轴反的。

如果复杂的话，可以在ButtonDownFcn中get(0,'PointerLocation')
或者get(gcf,'Position');

=================================== - [返回]
49).Matlab中有阶乘函数吗？
:#FangQ(Qianqian.Fang@dartmouth.edu), 2002/7/23. BigGreen/en_Matlab#

matlab没有提供直接的阶乘函数，但可以用prod(1:n)来求n!
用prod(1:2:2n-1)或者prod(2:2:2n)来求解n!!

=================================== - [返回]
50)..怎样才能把Maltab学精？
:#energy(长白山), 2002/9/23. SMTH/MathTools#

"带着问题学，活学活用，学用结合，急用先学，立竿见影，
在'用'字上狠下功夫。"

=================================== - [返回]
51)..Matlab如何计算大阶乘？
:#lll, 2002/9/18. SMTH/MathTools#

如果只需要大致的值，取log10,
计算出result=log10(1000!)=log10(1)+log10(2)+...log10(1000)
然后求10^result=10^result的小数部分*10^result的整数部分


=================================== - [返回]
52)..Matlab中怎样求变上限二重积分？
:#hyphone,2002/9/22, SMTH/MathTools #

对于解析函数，用两次int即可。
如求x+y在0[返回]
53)..用符号积分算出来Ei是什么意思，怎样求值？
:#hyphone,2002/9/22, 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,2002/10/21,comp.soft-sys.matlab#
:#Peter J. Acklam(pjacklam@online.no),2002/6/26,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 (slord@mathworks.com),2002-10-22,comp.soft-sys.matlab#

使用computer命令：
[C,MAXSIZE] = computer

=================================== - [返回]
56)..如何使用整型矩阵来节省内存？(Duane Hanselman)
:#Duane Hanselman (masteringmatlab@yahoo.com),2002/10/22,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),2002/11/10,SMTH/MathTools#

如果不想安装补丁，只能使用matlab -nojvm的形式启动
matlab，否则需要参照如下页面的解决方案来安装补丁：
http://groups.google.com/groups?hl=en&lr=&ie=UTF-8&oe=UTF-8&selm=
NQ8a9.6835%24ob2.611653%40newsread1.prod.itd.earthlink.net
http://www.mathworks.com/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.m

=================================== - [返回]
59).Matlab/Mathematica中如何中断当前运算？(FangQ)
:#FangQ(Qianqian.Fang@dartmouth),2002/11/19,BigGreen/MathTools#

Mathematica中使用: Alt+./Alt+,
Matlab中使用: Ctrl+C

=================================== - [返回]
60)).Matlab/Mathematica中如何检查括号匹配？(FangQ)
:#FangQ(Qianqian.Fang@dartmouth.edu),2002/11/11,SMTH/MathTools#

Mathematica中使用: Ctrl+.
Matlab editor中使用: Ctrl+B

=================================== - [返回]
61)..Matlab的GUI中为何无法使用uicontrol的句柄？(FangQ)
:#FangQ(Qianqian.Fang@dartmouth.edu),2002/11/11,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,2002/11/9,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),2002/11/19,refer to comp.soft-sys.matlab#
使用selectmoveresize函数，例如：

figure
h=uicontrol('style','pushbutton');
set(h,'ButtonDownFcn',...
'selectmoveresize;set(h,''selected'',''off'')',...
'Enable','inactive')

=================================== - [返回]
64)..Matlab如何求解维数巨大的稀疏矩阵方程？(FangQ)
:#FangQ(Qianqian.Fang@dartmouth),2002/11/19,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