24小时热门版块排行榜    

查看: 1312  |  回复: 1

在捂汗

铁虫 (小有名气)

[求助] matlab中怎么给四条曲线作出图例啊?求大神指导 已有1人参与

clear all; clc
disp('插值求取平滑曲线后的数值' )
t = [5 10 15 20 30 40 60 80 120] ;
A = [29.84 40.01 46.71 53.87 63.32 67.01 67.33 67.54 67.87];

knots = 3; K = 3 ;          % 三次B样条
sp = spap2(knots,K,t,A);
pp = fnder(sp);         % 计算B样条函数的导函数
dAdt = fnval(pp,t);    % 计算t处的导函数值
rAm = dAdt;

% 绘制浓度拟合曲线
ti = linspace(t(1),t(end),200);
Ai = fnval(sp,ti);
plot(t,A,'ro',ti,Ai,'b:')
hold on
t = [5 10 15 20 30 40 60 80 120] ;
B = [36.32 48.21 55.81 62.21 69.81 72.39 72.72 72.91 73.92];
knots = 3; K = 3 ;          % 三次B样条
sp = spap2(knots,K,t,B);
pp = fnder(sp);         % 计算B样条函数的导函数
dBdt = fnval(pp,t);    % 计算t处的导函数值
rBm = dBdt;
ti = linspace(t(1),t(end),200);
Bi = fnval(sp,ti);
plot(t,B,'rx',ti,Bi,'b:')
hold on
t = [5 10 15 20 30 40 60 80 120] ;
C = [41.91 55.62 64.13 70.46 77.64 79.80 80.04 80.57 83.92];

knots = 3; K = 3 ;          % 三次B样条
sp = spap2(knots,K,t,C);
pp = fnder(sp);         % 计算B样条函数的导函数
dCdt = fnval(pp,t);    % 计算t处的导函数值
rCm = dCdt;

% 绘制浓度拟合曲线
ti = linspace(t(1),t(end),200);
Ci = fnval(sp,ti);
plot(t,C,'rs',ti,Ci,'b:')
hold on
t = [5 10 15 20 30 40 60 80 120] ;
D = [42.16 62.13 72.91 79.87 84.21 85.61 85.84 86.21 86.22];

knots = 3; K = 3 ;          % 三次B样条
sp = spap2(knots,K,t,D);
pp = fnder(sp);         % 计算B样条函数的导函数
dDdt = fnval(pp,t);    % 计算t处的导函数值
rDm = dDdt;

% 绘制浓度拟合曲线
ti = linspace(t(1),t(end),200);
Di = fnval(sp,ti);
plot(t,D,'rv',ti,Di,'b:')
Xlabel('t')
ylabel('X_A')
legend(o,'100%',x,'90%',s,'80%',v,'70%')
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (职业作家)

Matlab专家

感谢参与,应助指数 +1
本帖仅楼主可见
2楼2015-03-09 20:48:58
已阅   申请程序强帖   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 在捂汗 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 295复试调剂 +7 简木ChuFront 2026-03-19 7/350 2026-03-20 14:11 by 星空星月
[考研] 本人考085602 化学工程 专硕 +18 不知道叫什么! 2026-03-15 20/1000 2026-03-20 13:52 by danney002
[考研] 一志愿武汉理工材料工程专硕调剂 +7 Doleres 2026-03-19 7/350 2026-03-20 10:39 by EBSD
[考研] 08工学调剂 +4 用户573181 2026-03-20 4/200 2026-03-20 10:32 by 朗月清风1
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 0817调剂 +3 没有答案_ 2026-03-14 3/150 2026-03-19 09:51 by Xu de nuo
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
信息提示
请填处理意见