clear all;
close all;
syms T;
syms k;
syms k1;
syms k2;
syms i;
A=3;
H=40;
syms t;
y1=t*(2*pi)^(-0.5)*exp(-(t-10)^2/2);
y2=(2*pi)^(-0.5)*exp(-(t-10)^2/2);
fy1=int(y1,k*T,k*T+T);
fy2=int(y2,k*T,k*T+T);
a=ceil((k*T-T)/0.5)*A/T;
b=ceil(k*T/0.5)*A;
c=ceil(k*T/0.5)*A/T;
for T=0.05:0.05:5
k2=ceil(H*0.5/(A*T));
if(T>=0.5)
k1=1;
i=round(T/0.05); % 此行增加round函数，保证i为整数
E(i)=symsum((a*(fy1-k*T)+b*fy2+c*(fy1-k*T)),k,k1,k2)+0+40*int((2*pi)^(-0.5)*exp(-(t-1)^2/2),k*T,k*T+T);
disp(E(i));
else
k1=ceil(0.5/T);
i=round(T/0.05); % 此行增加round函数，保证i为整数
E(i)=symsum((a*(fy1-k*T)+b*fy2+c*(fy1-k*T)),k,k1,k2)+0+40*int((2*pi)^(-0.5)*exp(-(t-1)^2/2),k*T,k*T+T);
disp(E(i));
end;
end
T = round(0.05:0.05:5); % 绘图x轴修改
plot( T,double(E) ) % 绘图S转换为double量，符号量不能绘图
xlim([0.05 5])
你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，目前超出我们的服务范围，暂时无法为您解答。
首次提问人员可免费体验一次有问必答服务。目前首次提问的问题服务范围为：编程语言、Java开发、python、数据库、前端开发 领域专业技术问题，为您提供问题的解决思路和指导。不提供源码代写、项目文档代写、论文代写、安装包资源发送或安装、软件使用指导等服务。
我们后续会持续优化，扩大我们的服务范围，为您带来更好地服务。