|
这样写好像可以,呵呵
function chapter7()
clc
clear all
close all
global tc w
tc=0.1;w=6.28;
final_time=1;
a=[0,final_time];
y0=0;
[t,y]= ode45(@eq4,a,y0);
[F,a]=eq4(t,y);
% plot(t,y) %
plot (t,y,t,F,':')
%% plot (t,y,t,F,':') % 这个地方就出错了,为什么找不到求解方程的F值那?可是
%% 运行时去掉F表达式后面的分号,在command窗口有F的求解结果,怎么就画不出来哪?
%% 希望明白人给个明白的解释,谢谢
xlabel('t'),ylabel('y(t)')
axis tight
end
function [F,ydot]=eq4(t,y)
global tc w
F=t.*exp(-t/tc).*cos(w*t);
ydot=5*(F-y);
end |
|