声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1306|回复: 5

[综合讨论] 请教如何绘制该函数图

[复制链接]
发表于 2008-4-10 13:55 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
请教高手如何绘制图中函数曲线图,编了如下计算,运行总出问题,请教高手指点一下。其中式中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)
QLV.JPG
回复
分享到:

使用道具 举报

发表于 2008-4-10 15:15 | 显示全部楼层

我来帮LZ贴错误吧

??? 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 编辑 ]
 楼主| 发表于 2008-4-10 16:00 | 显示全部楼层
呵呵,还未解决,期待高手进一步指点
发表于 2008-4-10 16:22 | 显示全部楼层
发表于 2008-4-10 16:27 | 显示全部楼层
nod啊,大不了后面直接添个0

[ 本帖最后由 eight 于 2008-4-10 17:28 编辑 ]
发表于 2008-4-10 16:46 | 显示全部楼层

回复 3楼 的帖子

求导的话,表达式不复杂,还是自己算的好
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-9 10:22 , Processed in 0.077301 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表