#matlab
#矩阵乘法
#信号采样
根据你的公式来看,应该是要用对应元素相乘,不是矩阵乘法,可将乘号改为.*
步长为5a
clc,clear;
%% 读取数据
excelFile='c\奥特曼\被检序列.xlsx'; % 文件路径
myData=xlsread(excelFile); % 读取数据
timeSeries=myData(:,1); % 时间序列数据
dataSeries=myData(:,2); %
dataCount=length(dataSeries); %%设置步长与检验值
step=5; % 步长
v=step+step-2; % 计算自由度
ttest=2.87; % 查表得t检验值,修改
len1=step;
len2=step;
x=timeSeries(step:dataCount-step);
for i=step:dataCount-step
n1=dataSeries(i-step+1:i);
n2=dataSeries(i+1:i+step);
mean1=mean(n1);
mean2=mean(n2);
c=(len1+len2)/(len1*len2);
var1=1/len1*sum((n1-mean1).^2);
var2=1/len2*sum((n2-mean2).^2);
delta1=len1*var1+len2*var2;
delta=delta1/(len1+len2-2);
t(i-step+1)=(mean1-mean2)/sqrt(delta*c);
end;
%% 制图
figure(1);
plot(x,t,'r-','linewidth',1.5);
xlabel('年份','FontName','TimesNewRoman','FontSize',12);
ylabel('统计量z','FontName','TimesNewRoman','Fontsize',12);
axis([min(x),max(x),-4,4]);
hold on
plot(x,0*ones(i-step+1,1),'-.','linewidth',1);
plot(x,ttest*ones(i-step+1,1),'-','linewidth',1);%更改数字
plot(x,-ttest*ones(i-step+1,1),':','linewidth',1);% 更改数字
hl=legend('t统计量','0.05显著水平','上限','下限');
set(hl,'Box','off');
应用案例
已知长时间水文序列带入检验
检验结果如图所示
在MATLAB中,矩阵乘法的维度要满足:左边矩阵的列数等于右边矩阵的行数 具体可以一起讨论