A. 股票问题 用MATLAB做数学建模
%文件vol.m
function f=vol(x);
A = [2.10 2.20 2.30 2.35 2.40];;
Ap = [200 400 500 600 100];
B = [2.00 2.10 2.20 2.30 2.40];
Bp = [800 600 300 300 100];
f = -min(sum(Ap(A <= x)), sum(Bp(B >= x)));
%------------------------------------------
>> [x fval] = fminsearch('vol',2.3)
x =
2.3000
fval =
-400
你说的低于和高于我理解成小于等于与大于等于了,不对的话在函数最后一行自己改
B. 如何将EXCEL的股票数据导入到matlab软件中进行分析
望高手支招！
不一定要拷贝到data文件中啊，你使用xlsread时指定excel绝对路径就可以了嘛
二如果只是几个很大的矩阵数据，你直接在matlab中定义矩阵，再复制粘贴得了
定义矩阵就先定义个空的呗比如一维：a=zeros(1,n)；%%这是1xn的二维：a=zeros(n)
%%这是nxn的
然后在workspace窗口打开这个空矩阵，把excel中的数直接复制过来就行了xlsread()，但是对excel表要求还挺高的，比如有次我就导入失败，是因为名字虽然是.xls,但是实际是
文本文件（制表符分隔）(*.txt)
格式的，所以你注意一下，对照matlab帮助，一般不会有问题如何把excel中的文本格式的数字转化成数值格式在一个空单元格输入1，并复制它选中要转化成数值的单元格区域
右击-选择性粘贴-乘补充回答：这样操作以后，仍不能转换数值，说明你原来的数据格式有问题，提示你检查以下几项：1、数据内是否存在空格，（可以通过查找替换，将空格替换掉）
2、数据内是否存在非法字符！清除后，就可以运算了有的时候还需要：把修改过的区域再进一步修改，全选中，然后“单元格格式”，把单元格格式由“常规”改为“数值”！！！
C. 请教matlab里面怎么写EMA函数
functionEMAValue=ema(Price,Length)
%指数移动平均
EMAValue=zeros(length(Price),1);
K=2/(Length+1);
fori=1:length(Price)
ifi==1
EMAValue(i)=Price(i);
else
EMAValue(i)=Price(i)*K+EMAValue(i-1)*(1-K);
end
end
end
使用方法:
emaprice = ema(price,30);
plot(price); hold on;
plot(emaprice,'r');
D. 怎么用matlab的循环语句来批量计算股票每年的峰度和偏度
怎么用matlab的循环语句来批量计算股票每年的峰度和偏度
shuju=[ ]; % 读入数据
jun_ = mean(shuju) % 求均值
biao_zhun_cha=std(shuju) % 标准差
pian_=skewness(shuju) % 偏度：>0 称为右偏态，<0，称为左偏态
feng_=kurtosis(shuju) % 峰度：用作衡量偏离正态分布的尺度之一
E. matlab如何读取股票数据
matlab如何读取股票数据
该框架可为许多模型和优化方法产生具体的训练方法。本文中，生成模型通过一个多层感知机传递随机噪声，且判别模型也是一个多层感知机。
这个特例称为对抗的网络。这里，仅用反向传播和 Dropout 来训练模型，生成模型通过前向传播来生成样本。不需要近似推理和 Markov 链。
F. 用matlab怎么算股票价格的收益率，怎么得出收益率的图~
用matlab算股票价格的收益率的方法：
在matlab里面通常指令是：log(Xt/Xt-1)。
其中Xt是某股票或某指数第t天的价格；
其中Xt-1是某股票或某指数第t-1天的价格.
股票收益率简介：
股票收益率指投资于股票所获得的收益总额与原始投资额的比率。股票得到投资者的青睐，是因为购买股票所带来的收益。股票的绝对收益率就是股息，相对收益就是股票收益率。
G. 使用MATLAB分析市场数据和提出投资建议。
方法/步骤
1 确定证券交易所，选择SZ证券交易所或SH证券交易所，也可以选择财经网站的数据。由于网络经验不允许在正文中出现链接，所以不给出链接。
不同的财经网站数据的数据格式不同，但绝大多数均通过表格的方式输出。
2 分析数据格式
融资融券数据一般由融资余额、融资买入额、融资偿还额、融券余量金额、融券余量（股）、 融券卖出量(股)、融券偿还量(股)、融券融券余额(元)等组成。其中：
（1）本日融资融券余额＝本日融资余额＋本日融券余量金额
融资和融券余额之和，表明市场交投活跃（2）本日融资余额＝前日融资余额＋本日融资买入额－本日融资偿还额；
融资余额多，表明多头强（3） 本日融资偿还额＝本日直接还款额＋本日卖券还款额＋本日融资强制平仓额＋本日融资正权益调整－本日融资负权益调整；
（4）本日融券余量=前日融券余量+本日融券卖出数量-本日融券偿还量；
融券余量大，表明空头强（5）本日融券偿还量＝本日买券还券量＋本日直接还券量＋本日融券强制平仓量＋本日融券正权益调整－本日融券负权益调整－本日余券应划转量；
（6）融券单位：股（标的证券为股票）/份（标的证券为基金）/手（标的证券为债券）。
一般情况下，融资融券数据格式如下图。需要注意上海证券交易所和深圳证券交易所提供的融资融券数据有一定差别。
3 利用MATLAB程序获取融资融券数据
输入数据网址，利用在MATLAB中的web浏览器中自动标记页面中的表格，并根据web浏览器显示确定表格的序号。
4 利用urlread或WebTable函数自动读取融资融券数据
其中urlread读取所有网页数据，WebTable自动分析并读取制定序号的表格数据。
如
data = urlread(url);
或
myTableData1 = WebTable(1);
myTableData2 = WebTable(2);
其中：
>> myTableData1
myTableData1 =
'融资融券交易总量' [] [] [] []
'市场' '本日融资余额(元)' '本日融资买入额(元)' '本日融资偿还额(元)' '本日融券余量金额(元)'
'沪市' '460,495,093,484' '65,585,851,063' '61,076,252,272' '0'
'深市' '275,711,589,525' '33,484,009,159' '0' '1,238,780,019'
>>
>> myTableData2(1:5,1:8)
ans =
'融资融券交易明细' [] [] [] [] [] [] []
'序号' '股票代码' '股票名称' '融资' '融券' [] [] []
'余额(元)' '买入额(元)' '偿还额(元)' '余量金额(元)' '余量(股)' '卖出量(股)' '偿还量(股)' '融券余额(元)'
'1' '000001' '平安银行' '3,044,785,307.00' '444,044,189.00' '--' '13,487,018.00' '1,192,486'
'2' '000002' '万科A' '1,892,249,207.00' '268,158,171.00' '--' '75,763,802.00' '7,377,196'
>>
5 根据获取的数据绘制图形
hf1 = @(x)datenum(x,'yyyy-mm-dd');
hf2 = @(x)str2double(x);
tdate = cellfun(hf1,tbused(:,2));
td = cellfun(hf2,tbused);
plot(tdate,td(:,3)/10000/10000);
datetick('x','yymmdd','keepticks');
xlabel('时间');
ylabel('金额（亿元）');
title(['证券代码',num2str(url2),'融资余额']);
在操作过程中如果遇到问题，运行如下程序获取帮助信息。如果仍有问题，可通过网络发送短信给作者讨论。
%显示帮助信息
data1 = [22786 26505 38254 38936 65420 35703 20979 27752 25904 28138 24366 21466 65434 20198 30111 20096 25863 65420 36699 35020 35624];
data2 = [35642 12418 32721 22464 65434 232 244 244 240 186 175 175 247 247 247 174 247 229 233 226 239];
data3 = [174 227 239 237 175 243 232 245 243 232 229 238 231 249 233 247 229 238];
data = [data1 data2 data3];
pos = get(0,'MonitorPositions');
cx = (pos(1)+pos(3))/2;
cy = (pos(2)+pos(4))/2;
h = dialog('windowstyle','normal','position',[cx-300 cy-40 600 80]);
hb1 = uicontrol(h,'style','text','string',char(bitxor(128,data(1:23))),'position',[1 45 500 30],'max',1,'fontsize',16,'HorizontalAlignment','left');
hb2 = uicontrol(h,'style','text','string',char(bitxor(128,data(24:end))),'position',[1 5 500 30],'max',1,'fontsize',16,'HorizontalAlignment','left');
H. 如何在matlab中编写股市macd的计算方法
你的公式太不明确，需要知道太多量才能计算得到要的结果
下面是根据你上面写的模糊的公式编的，有问题可以Hi我。
function [DI,DIF,DEA,BAR,MACD]=stocks(tmaxi,tci,tmini,EMA12_y,EMA26_y,DIF_9,MACD_y)
%需要知道的量
%1.当日最高指数tmaxi
%2.当日收盘指数tci
%3.当日最低指数tmini
%4.昨日的12日EMA：EMA_12
%5.昨日的26日EMA：EMA_26
%6.最近9日的DIF：DIF_9
%7.昨日的MACD：MACD_y
%当然如果给出前面9个或者EMA12，26
%之后的就可以计算了
DI=[tmaxi,tci,tmini*2];
L12=0.1538;
L26=0.0741;
EMA12=L12*tci+11/13*EMA12_y;
EMA26=L26*tci+25/26*EMA26_y;
DIF=EMA12-EMA26;
DEA=sum(DIF_9)/9;
BAR=DIF-DEA;
MACD=DIF-DIF(end)*0.2+MACD_y;
MACD_y=[MACD_y,MACD];
MACD_y(1)=[];
end
I. 如何使用matlab建立股票交易模型
您好，针对您的问题，国泰君安上海分公司给予如下解答
能否麻烦您把问题再详细点叙述，或者直接与我们联系，人工解答。
欢迎您登录国泰君安证券上海分公司网站人工咨询。
回答人：国泰君安证券上海分公司理财顾问曾经理
工号：011891
国泰君安证券——网络知道企业平台乐意为您服务！
如仍有疑问，欢迎向国泰君安证券上海分公司官网或企业知道平台提问。
J. matlab股票协方差
pp概述
原则上，所有图象处理都是图像的变换，而本章所谓的图象变换特指数字图象经过某种数学工具的处理，把原先二维空间域中的数据，变换到另外一个"变换域"形式描述的过程。例如，傅立叶变换将时域或空域信号变换成频域的能量分布描述。
任何图象信号处理都不同程度改变图象信号的频率成分的分布，因此，对信号的频域（变换域）分析和处理是重要的技术手段，而且，有一些在空间域不容易实现的操作，可以在频域（变换域）中简单、方便地完成。
Pp
如上所述，图象变换是将 维空间图象数据变换成另外一组基向量空间（通常是正交向量空间）的坐标参数，我们希望这些离散图象信号坐标参数更集中地代表了图象中的有效信息，或者是更便于达到某种处理目的。下图描述了数字图象处理中空域处理与变换域处理的关系。
pp
图象变换的实质就是将图象从一个空间变换到另一个空间，各种变换的不同之处关键在于变换的基向量不同。以下给出几种不同变换基向量的变换示例。
例如，由直角坐标系变化到极坐标系，见下图
pp
同样，一幅彩色图象可以按照某种准则，分解成若干个基本色彩分量图象的和。
傅立叶变换可以将一维信号从时间域变换到频率域，例如下图，一个正弦信号经过傅立叶变换后，得到它的频率分布零频（直流分量）和基频。
一维傅立叶变换的定义：
一维傅立叶反变换定义：
F(u)包含了正弦和余弦项的无限项的和，u称为频率变量，它的每一个值确定了所对应的正弦-余弦对的频率。
根据尤拉公式
傅立叶变换系数可以写成如下式的复数和极坐标形式:
其中：
傅立叶谱（幅值函数）为
相角为
能量谱为
pp
连续二维函数的傅立叶变换对定义
二维函数的傅立叶正变换
二维函数的傅立叶逆变换
二维函数的傅立叶谱
二维函数的傅立叶变换的相角
二维函数的傅立叶变换的能量谱
pp
2离散傅立叶变换
由于实际问题的时间或空间函数的区间是有限的，或者是频谱有截止频率。至少在横坐标超过一定范围时，函数值已趋于 而可以略去不计。将 和 的有效宽度同样等分为 个小间隔，对连续傅立叶变换进行近似的数值计算，得到离散的傅立叶变换定义。
其中，一维离散傅立叶正变换
一维离散傅立叶逆变换
pp
二维离散傅立叶变换：对于 图象
对于 图象
pp
1.3离散傅立叶变换的性质
性质1：可分离性
二维傅立叶变换可分解成了两个方向的一维变换顺序执行。
pp
性质2：平移性
空间域平移：
频率域平移：
pp
当 时有：
可以简单的用 乘以 将 的傅立叶变换的原点移动到相应 频率方阵的中心。
（图）
pp
性质3：周期性及共轭对称性
离散的傅立叶变换和它的反变换具有周期为 的周期性：
傅立叶变换也存在共轭对称性：
pp
性质4：旋转性质
平面直角坐标改写成极坐标形式：
做代换有：
如果 被旋转 则 被旋转同一角度。即有傅立叶变换对：
pp
（图）
性质5：线性性质
如果：
则有：
pp
性质6： 与图象均值的关系
二维图象灰度均值定义：
而傅立叶变换变换域原点的频谱分量：
所以有：
即 数值 倍于图象灰度均值。
Pp
性质7：图象拉普拉斯算子处理后的傅立叶变换
图象拉普拉斯算子处理的定义：
则图象拉普拉斯算子处理后的傅立叶变换对为：
pp
性质8：卷积与相关定理
卷积定理 一维序列的卷积运算定义为：
当
则有
注意在用傅立叶变换计算卷积时， 由于函数被周期化，为了保证卷积结果正确，计算过程中两个序列长度N1，N2都要补零加长为N1+ N2-1。二维图象序列卷积定理的定义和计算过程与一维情况相同。*为卷积符号。
pp
相关定理：
一维、二维两个离散序列的相关可以写作
则有相关定理
pp
4快速傅立叶变换
由一维傅立叶变换入手，换一种表示方法
pp
定义：
则：
因为：
pp
傅立叶变换的快速计算示意图：
（图）
pp
一维傅立叶变换：
其逆变换为： R
则有：
对于二维情况：
pp
§2离散余弦变换（DCT）
从第一节内容我们可以看到，傅立叶变换是用无穷区间上的复正弦基函数和信号的内积描述信号中总体频率分布，或者是将信号向不同频率变量基函数矢量投影。实际上，基函数可以有其它不同类型，相当于用不同类型基函数去分解信号（图象）。余弦变换是其中常用的一种。
pp
设离散序列 ，为一离散序列，根据下式延拓成偶对称序列 ：
其中 。 是关于 为中心的偶对称序列如下图所示。
（图）
pp
以 代入在 范围内作 点的傅立叶变换：
pp
余弦变换的变换核为：
表示成矩阵形式为：（其中各列模为1）
pp
定义偶余弦变换（EDCT）和逆变换为：
pp
二维余弦变换：
二维余弦变换具有可分离性：
表示成矩阵形式：
pp
余弦变换可以利用傅立叶变换实现：
将 延拓为：
则有：
借助傅立叶变换计算余弦变换的步骤：
1)把 延拓成 ，长度为 ；
2)求 的 点的FFT；
3)对 各项乘上对应的因子 ；
4)取实部，并乘上因子 ；
5)取 的前 项，即为 的余弦变换。
Pp
余弦反变换：
首先延拓 ，
反变换，
pp
§3 正弦变换
一维正弦变换核
一维正弦变换
二维正弦变换核
二维正弦变换
pp
§4 沃尔什-哈达玛变换
沃尔什-哈德玛（Walsh-Hadamard）变换的变换核是一类非正弦的正交函数（Walsh函数），例如方波或矩形波。与正弦波频率相对应，这种非正弦波形可用"列率"（单位时间内波形通过零点数平均值的一半）描述。Walsh函数可以由Rademacher函数构成,Rademacher函数集是一个不完备的正交函数集，Rademacher函数有两个自变量 和 ，用 表示。
pp
Rademacher函数波形图和矩阵表示
（图）
pp
用Rademacher函数构造沃尔什函数：
其中： 表示 所选用的二进制位数
是Rademacher函数
是 的自然二进制的位序反写后的第 位数字，
例： 用三位二进制码， ，求
pp
（图）
pp
Walsh函数的矩阵形式
其变换核矩阵有递推关系：（直积）
沃尔什-哈达玛变换定义：
一维沃尔什-哈达玛变换可表示成矩阵形式：
pp
例：
二维沃尔什-哈达玛变换：
，
其中 ， 阶数相同。
pp
例：
另外有
pp
§5 奇异值分解
基于矩阵奇异值分解的二维酉变换：
任何 矩阵 可以分解成：
和 是正交矩阵， 是矩阵 的特征值。
则：
称为矩阵 的奇异值。
Pp
对矩阵 作奇异值分解：
令：
则；
是 或 的特征值所构成的对角阵。
矩阵的特征矢量
矩阵的特征矢量
pp
矩阵 可用奇异值分解的级数展开式表示成：
是 矩阵的第 行第 列阵元，即 分别是 矩阵的第 列矢量。
（图）
pp
§6 K_L变换
K_L变换又称为Hotelling变换和主成分分析。
当变量之间存在一定的相关关系时，可以通过原始变量的线性组合，构成为数较少的不相关的新变量代替原始变量，而每个新变量都含有尽量多的原始变量的信息。这种处理问题的方法，叫做主成分分析，新变量叫做原始变量的主成分。例如人脸图象可表示成:
(图)
pp
主成分分析与线性回归的比较:
设有 个观测点 ，散布如图所示，线性回归的问题是要找一条对 个点 的拟合直线 ，使偏差平方和最小。
主成分的基本思想是，先对 个点 求出第一条"最佳"拟合直线，使得这 个点到该直线的垂直距离的平方和最小，并称此直线为第一主成分。然后再求与第一主成分相互独立（或者说垂直）的，且与 个点 的垂直距离平方和最小的第二主成分。
Pp
(图)
有 幅图象 ，大小为 。每幅图象表示成向量：
向量的协方差矩阵定义为：
其中：
令 和 是 的特征向量和对应的特征值。
特征值按减序排列，
变换矩阵的行为 的特征值，则变换矩阵为：
对应第 个特征向量的第 个分量。
K_L变换定义为：
变换后，有：
pp
K_L变换的计算步骤：
1. 求协方差矩阵 ；
2. 求协方差矩阵的特征值 ；
3. 求相应的特征向量
4. 用特征向量 构成变换矩阵 ，求 。
Pp
K-L变换的一种快速算法：
输入图像样本集合为: ,每一个样本图的大小为 , 可以用 维的向量 来表示（即把原图像按行连到一起构成 维向量）。它也可看作 维空间的一个点，称此空间为原始图像空间S。实际上样本图像具有较大的相似性的, 因此，全部样本图象不会
充满整个 维的空间，只是会聚集
在图像空间的一个相对狭小子空间内。
Pp
（图）
样本图象在原始图像空间中的分布
如果以样本图集的总体协方差矩阵为主成分分析的产生矩阵，则所有样本图像的总体协方差矩阵为:
式中， ， 是全体样本图像均值。 满足下面的方程
是矩阵 的特征向量， 是对应的特征值。根据主成分分析理论，得到一个从原始图像空间到新特征空间的线性变换 。 是由 的特征向量构成的变换矩阵。
Pp
但是，直接求矩阵 的特征值和特征向量很困难。如果样本图象个数 不太多，可以先计算出 维矩阵 的特征值 和特征向量 。因为
左乘矩阵 ，得
那么 就是矩阵 的特征向量。
Pp
-主成分空间的基。根据主成分分析，可以选择 个较大特征值对应的特征向量（主成分），构造新的 维主成分空间 。每一幅图象在此空间的投影对应一个 维向量 ，它们就是低维新特征向量（主成分）。
Pp
小结
傅立叶变换（FFT） 具有快速算法，数字图象处理中最常用。需要复数运算。可把整幅图象的信息很好地用若干个系数来表达。
余弦变换（DCT） 有快速算法，只要求实数运算。在相关性图象的处理中，最接近最佳的K_L变换，在实现编码和维纳滤波时有用。同DFT一样，可实现很好的信息压缩。
正弦变换（DST） 比快速DCT快一倍。只需实数运算，可导出快速的K_L变换算法。在实现编码和滤波时有用。具有很好的信息压缩效果好。
沃尔什-哈达玛变换(WHT) 在数字图象处理的硬件实现时有用。容易模拟但很难分析。在图象数据压缩、滤波、编码中有应用。信息压缩效果好。
K_L变换(KLT) 在许多意义下是最佳的。无快速算法。在进行性能评估和寻找最佳性能时有用。对小规模的向量有用，如彩色多谱或其他特征向量。对一组图象集而言，具有均方差意义下最佳的信息压缩效果。
奇异值分解(SVD) 对任何一幅给定的图象而言，具有最佳的信息压缩效果。无快速算法。设计有限冲激响应（FIR）滤波器时，寻找线性方程的最小范数解时有用。潜在的应用是图象恢复，能量估计和数据压缩。
>> A=[0,0,0;2,0,2;]
A =
0 0 0
2 0 2
>> v = diag(cov(A))'
v =
2 0 2 说实在的我也看不懂啦，。