矢量场

clear all; close all; clc;

h = 0.25;

tlist = 0 : 0.3: 2*pi;

Clist = -5: 0.1 : 5 ;

h1 = figure;

hold on

for s1 = 1 : length(tlist)

    for s2 = 1 : length(Clist)

        t = tlist(s1);

        C = Clist(s2);

        if abs(C + sin(t)) < 2

            quiver(t, C+ sin(t),  h, h * cos(t), "b")

        end

    end

end

plot(tlist, sin(tlist),"r")