matlab基础学习3

PLOT函数

x=[1.5 3 8 4]; plot(x)
绘制的是由（1，1.5）（2，3) (3,8) (4,4)这几个点连成的直线

横坐标是x向量几个元素一一对应的位置序号1，2，3，4

纵坐标是x向量

此时x向量为横坐标，y向量为纵坐标

TXT数据读入

读取纯数字文件

读取纯字符串的文本
1.记事本的编码要改为ansi

2.
（这里注意:‘r’是只读，不能用‘w’，会清除原txt文件中的内容，编码一定要改为ansi）

再继续下去就会得到第四行，第五行…的内容
附原txt文件



第一行是将tline这个字符串以’,'为分割点，分成若干个字符串
第二行是将分割后形成的三个字符串分别赋给date，lat和long

text文本数据写

多项式计算

多项式;

在matlab中：
p(x)表示为向量形式：[a_n,a_n-1,a_n-2,…,a₁,a₀]
关于如何打出上下标
多项式乘法
conv（p₁,p₂)其中p1，p2为多项式的系数向量
返回的结果也是一个系数向量
多项式除法
[q,r]=deconv(p₁,p₂)
q返回p₁除以p₂的商式
r返回p₁除以p₂的余式
q，r仍是多项式系数向量
----------------------------分隔--------------------------
由deconv是conv的逆运算，得
p₁=conv（q，p₂）+r
多项式求导
1.polyder（p）求多项式p的导函数
2.polyder（p,q)求多项式p和q乘积的导函数
3.[m,n]=polyder(p,q)求p/q的导函数，导函数的分子存入m，分母存入n
例：
求f(x)=3x³+x²-6和h(x)=x+2的乘积的导函数和商的导函数

%%
p=[3 1 0 -6];          %f(x)的系数向量
q=[1 2];                  %h(x)的系数向量
d1=polyder(p,q);      %二者乘积的导函数的系数向量
d1
[d2,d3]=polyder(p,q);     %f(x)/h(x）的导函数的分子，分母系数向量
d2                  %分子
d3                  %分母

多项式的求值
1.polyval(p,x)
p为多项式系数向量，x可以是标量，向量，矩阵
结果为多项式在x矩阵每个元素处所对应的函数值构成的矩阵
2.polyvalm(p,x)
多的字母m代表矩阵，所以这个函数要求x为方阵，以方阵为自变量，执行的运算为矩阵运算，不是点运算

注意一个是点运算，一个不是

(当x为标量时，二者计算结果相同）
注：size(A)返回一个行向量，第一个元素是矩阵的行数，第二个元素是矩阵的列数
多项式的求根
1.roots( p ) ,p为多项式的系数向量
2.若已知多项式的全部根，可用poly函数建立起该多项式。
p=poly（x），x为多项式的全部根

实例：
计算f(x)=-38.89x²+126.11x-3.42在[0,2]区间内的极值

%%
x=0:0.01:2;
y=-38.89*x.^2+126.11*x-3.42;      %函数方程
p=[-38.89,126.11,-3.42];              %函数的系数向量
daohanshu=polyder(p);                           %对函数求导         
gen=roots(daohanshu);                     %导函数的根
gen
jizhi=polyval(p,gen);          %根这点处的函数值，即极值
jizhi
plot(x,polyval(p,x),1.6214,98.8154,'rp');  %在极值点处标*号
text(gen,jizhi,'最大值')

平均数 mean
中位数 median
连乘 prod
累加和 cumsum
累乘积 cumprod
标准差 std
相关系数 corrcoef