主题:MATLAB繪圖問題
clear
clc
time=[0 10];
y=[0 1]';
[t1 y1]=ode45('finode3',time,y);
sum=0;
b=0;
for i=1:1:5
for k=0.0000000000001:1:i
a=1;
b=0;
b =b+((y1(max(find(t1<k)))+y1(max(find(t1<k+a))))*a)/2;
end
sum=sum+b
plot(t1,sum);
end
---------------------------------------------------------------------
副程式
function dy =finode3 (t,y)
dy(1)=y(2);
dy(2)=-y(1);
dy=dy';
end
-----------------------------------------------------------------------
我想繪製出時間和迴圈跑出SUM的那些值 0.3992,1.2677,1.8367 等等.......
但是我在workspace裡看到的sum只有一個值,而不是一個向量,不知道該怎麼會出圖出來
懇請各位大大可以教小弟一下,謝謝
clc
time=[0 10];
y=[0 1]';
[t1 y1]=ode45('finode3',time,y);
sum=0;
b=0;
for i=1:1:5
for k=0.0000000000001:1:i
a=1;
b=0;
b =b+((y1(max(find(t1<k)))+y1(max(find(t1<k+a))))*a)/2;
end
sum=sum+b
plot(t1,sum);
end
---------------------------------------------------------------------
副程式
function dy =finode3 (t,y)
dy(1)=y(2);
dy(2)=-y(1);
dy=dy';
end
-----------------------------------------------------------------------
我想繪製出時間和迴圈跑出SUM的那些值 0.3992,1.2677,1.8367 等等.......
但是我在workspace裡看到的sum只有一個值,而不是一個向量,不知道該怎麼會出圖出來
懇請各位大大可以教小弟一下,謝謝