论坛中常见的Matlab问题：
1.      数据生成问题

问题描述：如何生成两列数字x,y，使得x+y<=1？

解答：首先生成一列随机数x，第二列在(-inf,1-x)区间内取值。如有其他限制可以依次类推。

2.      数据导入问题

问题描述：希望将excel中的daily数据导入到matlab中，要求使用xlsread?

解答：可以先保存为字符串读入，再使用日期函数将字符串转化过来。

3.      符号方程求解：

问题描述：求解三次方程s=x^3+(a+2b-c)x^2-(a^2+3c-4d)x+(b^2+5c^2-2d)=0

解答：syms a b c x;

s=x^3+(a+2b-c)*x^2-(a^2+3*c-4*d)*x+(b^2+5*c^2-2*d);

solve(s);

4.      动态规划中0-1变量的定义

问题：我在用matlab求解动态规划，需要定义一个0-1变量，在matlab里面怎么表示呢？

解答：使用dummyvar。

5.      GARCH(1,1)模型

问题：如何在Matlab中实现GARCH(1,1)模型，另外Eviews和Matlab中拟合的结果不同，为何？

解答：Matlab中可以用金融工具箱拟合GARCH模型，拟合结果的差距更多是计算上的微小波动，并不影响结果的分析。

6.      双变量求和问题

问题：如何实现双变量求和，公式为？

答案：

result=(CT-C0)./(log(CT)-log(C0)).*log(IT./I0);



finalResult=sum(result(:));

7.      ARMA模型诊断

问题：如何用MATLAB确定ARMA模型的阶数？、

解答：通过自相关和偏自相关图。

8.      画图问题

问题：画y=(a*(b*a+1-b)-sqrt(a*(1-a)*(1-b*a)*(b*a+1-b)))/(2*a-1)的图像，定义域为[0.1]

解答：

b=0.5;



x=linspace(0,1);



y=(x.*(b.*x+1-b)-sqrt(x.*(1-x).*(1-b.*x).*(b.*x+1-b)))./(2*x-1);



plot(x,y);

设定不同的b就可以实现不同参数下的图像了。

9.      面板数据的单位根检验

问题：如何在Matlab中进行单位根检验？

解答：面板数据的单位根检验LLC、JPS-W检验、ADF-Fisher和PP-Fisher检验，需要根据检验原理编制程序。

不过LLC检验可以参考

http://onlyzhusy.spaces.live.com ... 504687420!606.entry。

10.  CSV文件读取

问题：What are the commands to manipulate with .csv files? Do the arrays start with 0? And how to play around with text format columns?

解答：Please try csvread().csvread() only supports numerial format,so if you need to use text format columns, try text or excel.

11.  portopt中矩阵的非正定问题

问题：如何解决portopt（）中的矩阵不是正的半定矩阵的问题啊(矩阵有随机数生成三角矩阵后转置填充得到)

解答：该矩阵必须是正定的，否则会导致拟合结果不稳定；参考的解决方案是使用AA’的形式。

12.  GARCH残差分布推广

问题：GARCH模型的残差性一般假定服从正态分布、t分布或广义误差分布，这些假定与金融市场的典型事实并不是很相符，请高手指点残差为稳定分布或双曲分布的GARCH模型的估计

解答：残差分布的推广后依然需要使用极大似然法进行估计，这是可以考虑蒙特卡洛模拟法来估计参数。

13.  如何将几列很长的数据在matlab中排列成矩阵

问题：如上

解答：reshape

14.  数据转换

问题：对于一个只有0，1两个元素的矩阵，在给定一个比例0.02的条件下，将0换成1，1换成0.

解答：

idx=rand(size(A));



A(idx<.02)=1-A(idx<.02);

15.  matlab求解规划问题

问题：如上

解答：linprog

16.  求解带参数矩阵的特征值 特征向量

问题：如上

解答：

syms x

a=(1-x)*eye(4)+x*ones(4,1)*ones(4,1)'

[V,D]=eig(a)

17. AR模型拟合

              问题：如上

              解答：

x=randn(100);

y=eye(100,1);

for i=2:100

    y(i)=1+0.9*y(i-1)+x(i);

end

plot(y)

18 参数估计问题

问题：有50维的数据，想分别用t、高斯、clayton、frank、gumble这几个copula函数来对其联合分布进行描述，怎么估计参数和与其对应的AIC值呢？

              解答：蒙特卡洛模拟

19 加权地理回归

              问题：如何运用Matlab进行地理加权回归分析（GWR）

              解答：按照原理编程。

20 包安装问题

              问题：如何将下载到的扩展包安装在matlab中

              解答：安装在WORK中

21 逻辑问题

              问题：在matlab中作两种条件当其中之一满足时，则执行某种操作。什么命令是表示或者关系的？

              解答：if else语句

22 相关系数问题

              问题：在matlab中：a=corrcoef（b）；那返回的a矩阵中a（i，j）是b中i行与j行的相关系数，还是i列与i列的相关系数？

              解答：行是观测 列是变量 所以这里的相关系数是指列之间的相关性

23 数据操纵问题

              问题：如何产生一个只有0和1的20*10的矩阵，且0和1的比例为1：1

              解答：

A=zeros(20,10);

idx=randperm(200);

A(idx(1:end/2))=1

24 拟合优度检验

              问题：MATLAB如何做检验拟合优度？

              解答：残差为res,数据为x，R=sqrt(1-res.^2/x.^2)

25 ARIMA模型定阶问题

              问题：请问用matlab做arima时，有没有什么命令可以自动定阶的？

              解答：自动定阶的也有，基本原理是尝试一定阶数范围内的ARIMA模型，根据一定的准则(AIC ,BIC等)选择某几个返还给统计员，由统计员最终决定使用哪个？

26 crosscorr问题

              问题：crosscorr函数的形式crosscorr(series1,series2,maxlag)

画出的图给出的lag区间为[-maxlag,maxlag]

但这个给出的相关系数为series1(t)与series2(t-lag)还是series1(t-lag)与series2(t)?

              解答：给出的是series1(t)与series2(t-lag)

27 规划求解

              问题：matlab目标方程线性,约束方程非线性,请问用optimization toolbox 哪个函数?

              解答：fmincon

28 GARCH定阶

              问题：如题

              解答：极大似然方法

29 密度函数绘制

              问题：如题

              解答：以正态分布为例

x=-10:0.01:10

y=normpdf(x)

plot(x,y)

30 数据操纵问题

              问题：我有一个m行2列的矩阵。所以，对于每一行，都有两个数。我只要那些行，它们的两个数都在-0.01到0.01之间的。其它不满足这条件的行，都不要，建立一个新的数据矩阵。

              解答：如下

m=10;A=randn(m,2)

A(all(A>-0.01&A<0.01,2),2)

31 二维正态分布问题

              问题：问高手用matlab语言如何生成二维标准正态分布数据！有现成的matlab语言么？X=normrnd(0,1,2,T)正确么？

解答：参考mvnrnd

32 绝对值拟合问题

              问题：请问matlab如何实现均值绝对偏差（mean absolute deviation）最小（含约束条件的）。这个应该属于线性规划，可不知道用哪个函数。请指教！

              解答：|x|=max(x,0)-min(x,0)，这样就可以写函数进行规划了，但是需要注意这里是非线性规划。

33 相同元素提取

              问题：有关matlab的问题---怎样从两个矩阵中取出相同元素

比如说有两个矩阵A=[1 2 3

                                   4 5 6

                                   7 8 9]

                        B=[ 4 10 12 7

                               7 5 25 6]

现在要把A和B中相同的元素提取出来组成一个新的矩阵，可以是任意行列的矩阵，请问有没有库函数可用？要是没有该怎么做？

              解答:

A=A(: )

B=B(: )

C=intersect(A,B)



34 matlab求解非线性方程组

问题：我想解出C11,C12,C21,C22,用以下方法



>> eq1='c11*c11+c12*c12=0.000449';

>> eq2='c11*c21+c12*c22=0.000181';

>> eq3='c21*c21+c22*c22=0.000753';

>> eq4='-411.693*c11+62.1194*c21=0';

>> solve(eq1,eq2,eq3,eq4)



得出的结果却是



ans =



c11: [4x1 sym]

c12: [4x1 sym]

c21: [4x1 sym]

c22: [4x1 sym]



好象不行啊!怎么才能解出C11,C12,C21,C22.

解答：1.fsolve : solves systems of nonlinear equations of several variables.



35 Step1:write an m-file



function f = weilai(x)



f(1) = x(1)*x(1)+x(2)*x(2)-0.000449;



f(2) = x(1)*x(3)+x(2)*x(4)-0.000181;



f(3) = x(3)*x(3)+x(4)*x(4)-0.000753;



f(4) = -411.693*x(1)+62.1194*x(3);



% save as weilai.m in a directory on the MATLAB path.





Step2: at the MATLAB command prompt. Enter



x0=[1,1,1,1]; %make a start guess at the solution



x=fsolve('weilai',x0)





x =





0.0041 0.0208 0.0273 0.0033





Step3.答案会随起始点不同而不同,请设定为你认为最适当的值.

36 数据导入

问题：matlab中如何调用Excel数据？

解答：xlsread

37 Moran's I指数

问题：              如何用MATLAB求moran'I 指数,有高手指点一、二吗?多谢.

解答：moran's I指数只是简单的最小二乘法，要给出目标值、空间权重，带入公式即可。

38 数据切割

问题：如何取矩阵中的某一行某一列

解答：A(i,:)取得第i行，A(:,i)取得第i列

39 多元正态随机数生成

问题：如何产生多元正态随机数？

解答：mvnrnd

做的好！

这个很好 可以多学点东西 呵呵