作业帮能解决的话，用maple也行。
来源：学生作业帮 编辑：拍题作业网作业帮 分类：综合作业 时间：2024/05/01 14:51:59
能解决的话，用maple也行。
如果matlab要很多个文件才能解决的话请分开写下：）
如果matlab要很多个文件才能解决的话请分开写下：）
clear
clc
s=[1,-1,2,4,2,1.5,2.5];
N=7;
a=1;b=7;
syms x p
fai_x=-x^2+3*x-3/2;
jj=0:4;
cp_x_p=sum(cos(2*p*j*pi/N))*fai_x;
for ii=0:6
cp_ip(ii+1)=limit(cp_x_p,ii); %cp_ip是i=0:6时CP(i,p)的7个值,为下面计算FCP(P)做准备
end
fcp_p=((b-a)/N)*sum(s.*cp_ip);
ccp_p=int(cp_x_p*cp_x_p,a,b); %int()为定积分
ap_p=fcp_p/ccp_p;
for p=0:3
Ap(p+1)=limit(ap_p,p); %计算p=0:3时AP(P)的4个值
end
AP=double(Ap)
将以上保存为脚本 运行后得
AP =
-0.0777 -0.0543 -0.0251 -0.0105
>>
不知道理解正确没有.
clc
s=[1,-1,2,4,2,1.5,2.5];
N=7;
a=1;b=7;
syms x p
fai_x=-x^2+3*x-3/2;
jj=0:4;
cp_x_p=sum(cos(2*p*j*pi/N))*fai_x;
for ii=0:6
cp_ip(ii+1)=limit(cp_x_p,ii); %cp_ip是i=0:6时CP(i,p)的7个值,为下面计算FCP(P)做准备
end
fcp_p=((b-a)/N)*sum(s.*cp_ip);
ccp_p=int(cp_x_p*cp_x_p,a,b); %int()为定积分
ap_p=fcp_p/ccp_p;
for p=0:3
Ap(p+1)=limit(ap_p,p); %计算p=0:3时AP(P)的4个值
end
AP=double(Ap)
将以上保存为脚本 运行后得
AP =
-0.0777 -0.0543 -0.0251 -0.0105
>>
不知道理解正确没有.