matlab - IF 和 FOR 循环不打印正确的值
问题描述
我有一个可能很愚蠢的问题,但我无法克服它。
这一切的目的是让两个变量 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
解决方案
我不得不更改您的代码上的一些内容。
看来您在理解 for 循环在 Matlab 中的工作方式时遇到了一些麻烦。
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
推荐阅读
- matplotlib - 具有不同标签/类别的区域之间的整齐轮廓
- angular - Angular 嵌套表单组件
- jmeter - Jmeter command line set heap size
- homebrew - 自制错误“kafka:没有可用的瓶子!” 我的 m1 macbook 上的第一个错误
- android - 未使用数据/主机意图过滤器调用 Android 接收器
- java - 如何将android studio stacktrace日志保存到设备内部存储中的文本文件中
- qt - 转换后如何访问我的自定义 QGraphicsItem 的数据属性的内容
- reactjs - OTT平台如何限制用户(如3流/账号)【设备限制】
- r - RStudio 不会在 Anaconda-Navigator 中启动
- snowflake-cloud-data-platform - 使用映射表中的多个字段的雪花行级策略