IF 和 FOR 循环不打印正确的值

Question

我有一个可能很愚蠢的问题，但我无法克服它。

这一切的目的是让两个变量F_hor_alb和F_ver_alb根据if语句中报告的规则来定义。

如果该角度的余弦 (60wt) 在这种情况下大于等于 -0.37（阴影）的量，则这两个变量必须具有正弦趋势。否则它们的价值为零。

通过验证余弦和变量的输出，第一个是应该的，但变量的所有值都等于零，除了第一个，因此不遵循余弦本身假设的值。 IF 中有什么不工作的地方吗？

（Fhor 和 Fver 已经在另一段代码中定义）。

Example of output
Cos(60wt): 0,9980         0,9970 0,9960   .....
F_hor_alb: 0,9980*F_hor     0      0        0     0 0 0 0 0 0 0
F_ver_alb: 0,9980*F_ver     0      0        0     0 0 0 0 0


Re = 6378; % Radius of Earth (km)
z = 500; % Orbital Altitude (km)
mu = 3.986004415e5 % Gravitational constan of Earth (km3/s2)
OrbitalT = 2*pi*sqrt(((Re+z)^3)/mu); % Orbital Period (s)
orbits = 30; % Number of cycles for the simulation
simtime = (OrbitalT*orbits)/60; % Simulation time (min)
simstep = 1; % Sampling step (min)

w=2*pi/OrbitalT; 
F_hor_alb = zeros(ceil(simtime));
F_ver_alb = zeros(ceil(simtime));

for t = 1:simstep:ceil(simtime)
    if cos(w*t*60)>shadow
        F_hor_alb = F_hor*cos(w*t*60);
        F_ver_alb = F_ver*cos(w*t*60);
    else
        F_hor_alb = 0;
        F_ver_alb = 0;
    end

    F_hor_alb(t) = F_hor_alb;
    F_ver_alb(t) = F_ver_alb;
    coseno(t)=cos(w*t*60);
    t=t+1;
end

Answer 1

我不得不对您的代码进行一些更改。

您似乎在理解 Matlab 中 for 循环的工作原理时遇到了一些麻烦。

Re = 6378; % Radius of Earth (km)
z = 500; % Orbital Altitude (km)
mu = 3.986004415e5 % Gravitational constan of Earth (km3/s2)
OrbitalT = 2*pi*sqrt(((Re+z)^3)/mu); % Orbital Period (s)
orbits = 30; % Number of cycles for the simulation
simtime = (OrbitalT*orbits)/60; % Simulation time (min)
simstep = 1; % Sampling step (min)

w=2*pi/OrbitalT;

F_hor = 1;
F_ver = 2;
shadow = 0.5;

t = (1:simstep:ceil(simtime))';
F_hor_alb = nan(size(t));
F_ver_alb = nan(size(t));
coseno = nan(size(t));

for i = 1:length(t)
    if cos(w*t(i)*60)>shadow
        disp('here');
        F_hor_alb(i) = F_hor*cos(w*t(i)*60);
        F_ver_alb(i) = F_ver*cos(w*t(i)*60);
    else
        F_hor_alb(i) = 0;
        F_ver_alb(i) = 0;
    end

    coseno(i)=cos(w*t(i)*60);
end