Question

matlab矩阵积分

matlab.
有一个多变量的时间序列数据，F是其中几个列变量关于z的函数，对F(z)求积分。（假设有100行参数数据，则得到100行积分结果）
下面是自己写的前面部分，积分算不出来，请教高手~！

par= importdata('kou.txt');
no=par(:,1);
ws2=par(:,2);
ws10=par(:,3);
wd=par(:,4);
p=par(:,5);
t2=par(:,6);
t10=par(:,7);
vapor2=par(:,8);
vapor10=par(:,9);
ulr=par(:,10);
rnet=par(:,11);
soilt5=par(:,12);
vwc5=par(:,13);
sd=par(:,14);
% 各列变量倒入

ls=2.838*10^6;m=18.01;
langmdat=0.00063.*(t2+273)+0.0673;r=8313;
rous=0.611.*2.17./(t2+273).*exp(17.3.*t2./(t2+237.3));
d=2.06.*10^(-5).*((t2+273)./273).^1.75;
fi2=1-vapor2./(0.611.*exp(17.3.*t2./(t2+237.3)));
% define parmeters

ft=ls.*(ls.*m./(r.*(t2+273))-1)./(langmdat.*(t2+273))+1./(d.*rous);
% define ft

hf=0.08436.*0.02264.*ws10.^1.27;
uf=0.02264.*ws10;
% hf is the lower boundary for suspension

定义函数nu(z,ws10)和mt(z,ft,ws10,fi2)，
引入变量z，在z的区间0-6之间求函数F=nu(z,ws10)*mt(z,ft,ws10,fi2)的积分如何求？

Answer 1

假设你的nu=z^3*ws10 mt=sin(z)*ft.*w10.*fi2
那么F=z^3*ws10*sin(z).*ft.*w10.*fi2
积分这样写：
syms z ;这一行是定义符号z
int(z^3.*ws10*sin(z).*ft.*w10.*fi2,0,6)
;由于你的ft，ws10，fi2都已经导入了 也就是变量，可以直接用
;对了 补充一句 数组可以直接运算,但是要注意点乘和*的区别