matlab和Python线性规划

先说matlab，数据用csv文件，

打开数据

用csvread函数

注意：csvread函数只试用与用逗号分隔的纯数字文件

第一种：M = CSVREAD(‘FILENAME’) ，直接读取csv文件的数据，并返回给M

第二种：M = CSVREAD(‘FILENAME’,R,C) ，读取csv文件中从第R-1行，第C-1列的数据开始的数据，这对带有头文件说明的csv文件(如示波器等采集的文件)的读取是很重要的。

第三种：M = CSVREAD(‘FILENAME’,R,C,RNG)，其中 RNG = [R1 C1 R2 C2]，读取左上角为索引为(R1,C1) ,右下角索引为(R2,C2)的矩阵中的数据。

注意：matlab认为CSV第1行第1列的单元格坐标为（0,0）

例1.1 读取整个文件
csvread(‘csvlist.csv’)
例1.2 读取第2行以下，第0列以右区域的数据
m = csvread(‘csvlist.dat’, 2, 0)
例1.3 读取第2行以下，第0列以右，第3行以上，第3列以左区域的数据
m = csvread(‘csvlist.dat’, 2, 0, [2,0,3,3])

使用textscan函数

在使用textscan函数前必须用fopen函数打开CSV文件。textscan函数读取的结果会存在cell数组中。
help textscan
可以知道该函数可以通过控制参数，从文本文件或字符串文件中读取数据
C = textscan(FID,’FORMAT’) reads data from an open text file identified by FID into cell array C. Use FOPEN to open the file and obtain FID. The FORMAT is a string of conversion specifiers enclosed in single quotation marks（由单括号括起来的一串字符说明）. The number of specifiers determines the number of cells in the cell array C. （说明字符的数量决定了输出cell的数量）For more information, see “Format Options.”

C = textscan ( FID, ‘FORMAT’ ,N) reads data from the file, using the FORMAT N times, where N is a positive integer. To read additional data from
the file after N cycles, call textscan again using the original FID.

C = textscan (FID, ‘FORMAT’, ‘PARAM’, VALUE) accepts one or more comma-separated parameter name/value pairs. For a list of parameters and values, see “Parameter Options.”

C = textscan(FID,’FORMAT’,N,’PARAM’,VALUE) reads data from the file, using the FORMAT N times, and using settings specified by pairs of PARAM/VALUE arguments.

C = textscan(STR,…) reads data from string STR. You can use the FORMAT, N, and PARAM/VALUE arguments described above with this syntax. However, for strings, repeated calls to textscan restart the scan from the beginning each time.（重复调用时每次从文件首部开始） (To restart a scan from the last position, request a POSITION output. See also Example 3.)

[C, POSITION] = textscan(…) returns the file or string position at the end of the scan as the second output argument. For a file, this is the value that FTELL(FID) would return after calling textscan（将读取位置作为第二个参数输出）. For a string, POSITION indicates how many characters textscan read.（对于字符串而言，位置参数显示了以读取的字符数）

Supported values for SPECIFIER:

    Numeric Input Type   Specifier   Output Class    ------------------   ---------   ------------    Integer, signed        %d          int32                           %d8         int8                           %d16        int16                           %d32        int32                           %d64        int64    Integer, unsigned      %u          uint32                           %u8         uint8                           %u16        uint16                           %u32        uint32                           %u64        uint64    Floating-point number  %f          double                           %f32        single                           %f64        double                           %n          double 

其他格式可以看帮助文档
Examples:

Example 1: Read each column of a text file.    Suppose the text file 'mydata.dat' contains the following:        Sally Level1 12.34 45 1.23e10 inf Nan Yes 5.1+3i        Joe   Level2 23.54 60 9e19 -inf  0.001 No 2.2-.5i        Bill  Level3 34.90 12 2e5   10  100   No 3.1+.1i    Read the file:        fid = fopen('mydata.dat');        C = textscan(fid, '%s%s%f32%d8%u%f%f%s%f');        fclose(fid);    textscan returns a 1-by-9 cell array C with the following cells:        C{1} = {'Sally','Joe','Bill'}            %class cell        C{2} = {'Level1'; 'Level2'; 'Level3'}    %class cell        C{3} = [12.34;23.54;34.9]                %class single        C{4} = [45;60;12]                        %class int8        C{5} = [4294967295; 4294967295; 200000]  %class uint32        C{6} = [Inf;-Inf;10]                     %class double        C{7} = [NaN;0.001;100]                   %class double         C{8} = {'Yes','No','No'}                 %class cell        C{9} = [5.1+3.0i; 2.2-0.5i; 3.1+0.1i]    %class double

感谢：http://blog.sina.com.cn/s/blog_a0246c110101v1k3.html

拟合

cftool

可以进行拟合，很方便的看到拟合效果
1. 右键还能一键保存在底部的table of fits 中选中行右键，save … to Workspace
2. 也可以生成M文件，函数输出fitresult，和gofl两个量。其中fitresult为函数拟合的变量结果，用fitresult.参数名就可以查看参数。
3. 比如你想模拟 y = p（1）*x + p（2）；
通过lsqcurvefit函数得到了 p（1）和p（2）；
你把方程右边变成字符串类型: str = [num2str(p(1)),’*x+’,num2str(p(2))];
然后global y； y = inline（str）；
然后再另一个控件的callback中，global y； 假设自变量x=1，输入y（1）就会得到因变量y的值

命令行调用

fit:用曲线或曲面拟合数据
fitobject = fit(x,y,fitType)按fittype说明的参数对x,y进行拟合
fitobject = fit([x,y],z,fitType)三元
fitobject = fit(x,y,fitType,fitOptions)按fitoptions设置的参数用fittype对xy进行拟合
fitobject = fit(x,y,fitType,Name,Value)
[fitobject,gof] = fit(x,y,fitType)
[fitobject,gof,output] = fit(x,y,fitType)

感谢：http://cn.mathworks.com/help/curvefit/fit.html?s_tid=gn_loc_drop