matlab循环每次循环都把结果输出来,怎么把循环的每一步结果都保存下来啊

本帖最后由 芒点DK 于 2016-12-9 11:04 编辑

各位前辈请看我的代码

% clear workspaces

clear

clc

%define variables

rmin = 0;

rmax = 8;

N = 9;

dt = 1;

t = 0;

tmax = 10;

b = 1.18;

y2 = [100 90 80 70 60 50 60 60 80 0]';

dr = 1;

r = rmin : dr : rmax;

% set initial conditions

v = [1586163 2529096 3481993 4301050 5275240 6084689 7439647 8497063 9465949 0];

unp2 = v;

% loop through time

nsteps = tmax/dt;

for n = 1 : nsteps;

% calculate the FOU scheme of population

for i = 2 : N + 1

unp2(i) = v(i) - dt/dr*(v(i) - v(i-1));

end

% update t and u and v

t = t + dt;

q = 1 : 1 : 9;

v = 0.9*unp2;

% calculate boundary condations

F = b*0.001*v*y2;

v(1) = 0.9*F;

v(N) = 0;

% plot solution

hold on;

a2 = sum(v(7:9))/sum(v);

if n == 1

hline = plot(t, a2, 'k-');

a1 = a2;

else

hline.Marker = 'none';

plot([t-1 t], [a1 a2], 'k-');

hline = plot(t, a2, 'k-');

a1 = a2;

end

axis([0 10 0 1])

shg

pause(0.01);

前辈请看，我把我的代码简化之后是这样，有两个问题想解决

一个是我想在F的位置分别乘以1.1 1.2 1.3 1.4

然后分别画四条曲线，怎么把它弄在一张图里呢？

第二个问题是我想把每一步的结果都保存下来，要怎么弄呢？

求前辈指点迷津！