main.m

来自「用matlab求两个向量的相关函数」· M 代码 · 共 20 行

clear;clc;
load f00d50.txt;
load f00d100.txt;
n1=500;n2=20;
s1=f00d50(1:n1,1);
s2=f00d100(1:n2,1);
%s2=s1;
n3=n1+n2;
s3=zeros(n3,1);

for k=1:n3
    tao(k)=k-n1;
    for i=1:n2
        j=i+tao(k);
        if i>=1&i<=n1&j>=1&j<=n2
        s3(k)=s3(k)+1/n3*s2(i)*s1(j);
        end 
    end
end
plot(tao, s3, 'DisplayName', 's3 vs tao', 'XDataSource', 'tao', 'YDataSource', 's3');