注册 登录
编程论坛 Matlab

请大家帮我看看这个程序该怎样修改

Tianbian1127 发布于 2008-12-14 14:04, 870 次点击
我按附件中的公式中的内容编了一个程序,对于两次的求和用matlab不知道应该怎样写,程序可以运行,但是得到的图却不对,最终得到的图应该是峰变窄,面积不变,大家帮我看看应该怎么修改,在下感激不尽!程序如下:
clear
clc
x(1)=40;x(2)=60;
m=2;
x=1:0.1:100;
f=50*exp(-(x-40).^2/100)+70*exp(-(x-60).^2/100);

N1=zeros(size(x));
N2=zeros(size(x));
N3=zeros(size(x));
N4=zeros(size(x));
t=zeros(size(x));

N=length(x);
for k=1:1:m
    t=x-x(k)
    for j=1:N
        if t(j)>=0 & t(j)<1;
            N1(k,j)=1/2*(t(j)).^2
        elseif t(j)>=1 & t(j)<2
            N2(k,j)=3/4-(t(j)-3/2).^2
        elseif t(j)>=2 & t(j)<=3
            N3(k,j)=1/2*(3-t(j)).^2
        else
            N4(k,j)=0
        end
    end
    N=N1+N2+N3+N4;
    F=N.*f;
end
plot(F)
3 回复
#2
Tianbian11272008-12-14 14:16
在线等高人啊,急啊
#3
neoxhl2008-12-14 14:53
2次微分应该有两个定义域的限制?
问题应该在
k=1:1:m
#4
Tianbian11272008-12-14 15:06
谢谢,还有就是公式里两个求和的运算我不知道怎么写进去,你能指点我一下吗?
1