请教如何绘制该函数图

hjydhu

请教高手如何绘制图中函数曲线图，编了如下计算，运行总出问题，请教高手指点一下。其中式中numbda在我计算过程中写为s。先多谢了
clear all
clc
syms x;
i=1;
for t=0:0.2:14
    s=0.5*t;
    C=[0.1946 0.0213 0.0125];
    S=[s s.^2 s.^3];
T=@(t) C.*S, % instantaneous elastic response
G=@(t) 1-0.562*(1-exp(-t/1.6483)), %relaxation function
g=@(x) T(t-x).*diff(1-0.562*(1-exp(-t./1.6483))), %integration function
F(i,1)=T+quad(g,0,t); % force response
i=i+1;
end
t=0:0.2:14;
plot(t,F,'c')
figure(2),plot(t,G)
figure(3),plot(t,T)

sigma665

??? Error using ==> times
Matrix dimensions must agree.

Error in ==> @(x)T(t-x).*diff(1-0.562*(1-exp(-t./1.6483)))


Error in ==> quad at 77
y = f(x, varargin{:});

diff以后，数据长度就变小了

[ 本帖最后由 sigma665 于 2008-4-10 15:18 编辑 ]

hjydhu

呵呵，还未解决，期待高手进一步指点

eight

原帖由 hjydhu 于 2008-4-10 16:00 发表

                               
登录/注册后可看大图

呵呵，还未解决，期待高手进一步指点

请根据2楼的提示自行解决

wwbeyondww

nod啊，大不了后面直接添个０

[ 本帖最后由 eight 于 2008-4-10 17:28 编辑 ]

sigma665

求导的话，表达式不复杂，还是自己算的好