2021-05-06

最新推荐文章于 2024-04-18 10:00:40 发布
最新推荐文章于 2024-04-18 10:00:40 发布
阅读量906 收藏 7
点赞数 1
分类专栏： 计算流体力学 CFD 安德森 文章标签： 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ErnestChai/article/details/116458651
版权

二维超声速流动的数值解——普朗特-迈耶稀疏波程序分享

本文主要分享了安德森的《计算流体力学基础及其应用》一书中：“二维超声速流动的数值解——普朗特-迈耶稀疏波” 的程序求解。水平有限，请见谅！

文章目录

Matlab

x = 12.6905 m 处中心稀疏波超声速流的结果

在这里插入图片描述

Prandtl_Meyer.m

function [u_t, v_t, rho_t, p_t, T_t, Ma_t, Delta_xi] = Prandtl_Meyer(x, u, v, rho, p, Ma)
%  [u_t, v_t, rho_t, p_t, T_t, Ma_t] = Prandtl_Meyer(x, u, v, rho, p)
% 二维超声速流动的数值解——普朗特-迈耶稀疏波；
%  由x处的流场参数计算 x + Delta_x 处的参数；
%  输入是(i,j)处的流场参数，输出是(i+1,j)处的流场参数；
%  u、u_t —— x方向的速度；
%  v、v_t —— y方向的速度；
%  rho、rho_t —— 密度；
%  p、p_t —— 压力；
%  Ma、Ma_t —— 马赫数；
%  Delta_xi —— 从(i+1,j)到(i+2,j)的推进步长 Delta_xi；

%% 各种常数；
gama = 1.4; % 比热比；
CFL = 0.5;% CFL常数；
R = 8.314;% 理想气体常数 [J/(mol*K)]；
M = 0.029; %空气的摩尔质量 [kg/mol]；
Rg= R / M; % 空气的气体常数 [J/(kg*K)]；
Cy = 0.6; %人工粘性常数；

%% 计算域；
J = 41; % y方向的计算节点数；
theta = 5.352/180 * pi;% 转角大小；
E = 10;% 转角位置；

% 物理平面；
ys = @(x) 0*(x>=0 & x<=E) + (-(x - E)*tan(theta))*(x>=E);% 下边界的纵坐标；
% h —— 物理平面从下表面到上边界的距离；
h = @(x) 40 * (x>=0 & x<= 10) + (40 + (x - 10)*tan(theta))*(x>=10 & x<=45);
Delta_y = h(x)/(J - 1);
y = ys(x):Delta_y:40;

% 计算平面；
eta = @(x,y) (y - ys(x))/h(x);
eta_min = 0;
eta_max = 1;
Delta_eta = (eta_max - eta_min)/(J - 1);

%% 度量；
partial_eta_x = @(x,y) 0*(x>=0 & x<=E)...
              + ((1 - eta(x,y))*tan(theta)/h(x))*(x>=E);

%% 推进步长 Delta_xi；
% 计算马赫角 mu;
mu = asin(1 ./ Ma);
% 计算推进步长 Delta_xi；
tan_max = max(max(abs(tan(theta + mu))), max(abs(tan(theta - mu))));
Delta_xi = CFL * Delta_y / tan_max;

%% 计算解向量 F 和 向量 G;
F1 = rho.*u;
F2 = rho.*u.*u + p;
F3 = rho.*u.*v;
F4 = gama/(gama - 1)*p.*u + rho.*u.*(u.*u + v.*v)/2;

G1 = rho.*F3./F1;
G2 = F3;
G3 = rho.*(F3./F1).^2 + F2 - F1.^2./rho;
G4 = gama/(gama - 1)*(F2 - F1.^2./rho).*F3./F1...
   + rho/2.*F3./F1.*((F1./rho).^2 + (F3./F1).^2);

%% 预估步计算 (i , j) 处的partial_F；
% 向前差分； 
partial_F1 = zeros(1,J);
partial_F2 = zeros(1,J);
partial_F3 = zeros(1,J);
partial_F4 = zeros(1,J);

for j = 2:J-1
    partial_F1(j) = partial_eta_x(x,y(j))*(F1(j) - F1(j+1))/Delta_eta...
                  + 1/h(x)*(G1(j) - G1(j+1))/Delta_eta;
              
    partial_F2(j) = partial_eta_x(x,y(j))*(F2(j) - F2(j+1))/Delta_eta...
                  + 1/h(x)*(G2(j) - G2(j+1))/Delta_eta;
              
    partial_F3(j) = partial_eta_x(x,y(j))*(F3(j) - F3(j+1))/Delta_eta...
                  + 1/h(x)*(G3(j) - G3(j+1))/Delta_eta;
              
    partial_F4(j) = partial_eta_x(x,y(j))*(F4(j) - F4(j+1))/Delta_eta...
                  + 1/h(x)*(G4(j) - G4(j+1))/Delta_eta;          
end

%% 计算预估步的人工粘性 SF；
SF1 = zeros(1,J);
SF2 = zeros(1,J);
SF3 = zeros(1,J);
SF4 = zeros(1,J);

for j = 2 : J-1
    SF1(j) = Cy*abs(p(j+1) - 2*p(j) + p(j-1))/(p(j+1) + 2*p(j) + p(j-1))...
              * (F1(j+1) - 2*F1(j) + F1(j-1));
          
    SF2(j) = Cy*abs(p(j+1) - 2*p(j) + p(j-1))/(p(j+1) + 2*p(j) + p(j-1))...
              * (F2(j+1) - 2*F2(j) + F2(j-1));
          
    SF3(j) = Cy*abs(p(j+1) - 2*p(j) + p(j-1))/(p(j+1) + 2*p(j) + p(j-1))...
              * (F3(j+1) - 2*F3(j) + F3(j-1));
          
    SF4(j) = Cy*abs(p(j+1) - 2*p(j) + p(j-1))/(p(j+1) + 2*p(j) + p(j-1))...
              * (F4(j+1) - 2*F4(j) + F4(j-1));
end

%% (i+1 , j) 处的预估值；
bar_F1 = F1 + partial_F1 * Delta_xi + SF1;
bar_F2 = F2 + partial_F2 * Delta_xi + SF2;
bar_F3 = F3 + partial_F3 * Delta_xi + SF3;
bar_F4 = F4 + partial_F4 * Delta_xi + SF4;

%% 求原变量 rho的 (i+1 , j) 处预估值 bar_rho； 
% A、B和C为中间变量； 
bar_A = bar_F3.^2 ./ (2 * bar_F1) - bar_F4;
bar_B = gama/(gama - 1) * bar_F1 .* bar_F2;
bar_C = - (gama + 1)/(2*(gama-1)) * bar_F1 .^3;

bar_rho = (- bar_B  + sqrt(bar_B.^2 - 4*bar_A.*bar_C)) ./ (2*bar_A);

%% 计算bar_p，用于计算校正步的人工粘性 bar_SF；
bar_u = bar_F1 ./ bar_rho;
bar_p = bar_F2 - bar_F1 .* bar_u;

%% 求Ｇ的 (i+1 , j) 处的预估值 bar_G;
bar_G1 = bar_rho .* bar_F3 ./ bar_F1;
bar_G2 = bar_F3;
bar_G3 = bar_rho .* (bar_F3 ./ bar_F1).^2 ...
            + bar_F2 - bar_F1.^2 ./ bar_rho;
bar_G4 = gama/(gama - 1)*(bar_F2 - bar_F1.^2 ./ bar_rho) .* bar_F3 ./ bar_F1...
            + bar_rho/2 .* bar_F3 ./ bar_F1 .* ((bar_F1./bar_rho).^2 + (bar_F3./bar_F1).^2);

%% 校正步计算 (i+1 , j) 处的 bar_partial_F;
% 向后差分；
bar_partial_F1 = zeros(1,J);
bar_partial_F2 = zeros(1,J);
bar_partial_F3 = zeros(1,J);
bar_partial_F4 = zeros(1,J);

for j = 2 : J-1
    bar_partial_F1(j) = partial_eta_x(x,y(j))*(bar_F1(j-1) - bar_F1(j))/Delta_eta...
                              + 1/h(x)*(bar_G1(j-1) - bar_G1(j))/Delta_eta; 
                          
    bar_partial_F2(j) = partial_eta_x(x,y(j))*(bar_F2(j-1) - bar_F2(j))/Delta_eta...
                              + 1/h(x)*(bar_G2(j-1) - bar_G2(j))/Delta_eta; 
                          
    bar_partial_F3(j) = partial_eta_x(x,y(j))*(bar_F3(j-1) - bar_F3(j))/Delta_eta...
                              + 1/h(x)*(bar_G3(j-1) - bar_G3(j))/Delta_eta; 
                          
    bar_partial_F4(j) = partial_eta_x(x,y(j))*(bar_F4(j-1) - bar_F4(j))/Delta_eta...
                              + 1/h(x)*(bar_G4(j-1) - bar_G4(j))/Delta_eta;                          
end

%% 校正步的bar_partial_F改为向前差分；
bar_partial_F1(1) = partial_eta_x(x,y(1))*(bar_F1(1) - bar_F1(2))/Delta_eta...
                              + 1/h(x)*(bar_G1(1) - bar_G1(2))/Delta_eta; 
                          
bar_partial_F2(1) = partial_eta_x(x,y(1))*(bar_F2(1) - bar_F2(2))/Delta_eta...
                              + 1/h(x)*(bar_G2(1) - bar_G2(2))/Delta_eta; 
                          
bar_partial_F3(1) = partial_eta_x(x,y(1))*(bar_F3(1) - bar_F3(2))/Delta_eta...
                              + 1/h(x)*(bar_G3(1) - bar_G3(2))/Delta_eta;   
                          
bar_partial_F4(1) = partial_eta_x(x,y(1))*(bar_F4(1) - bar_F4(2))/Delta_eta...
                              + 1/h(x)*(bar_G4(1) - bar_G4(2))/Delta_eta; 

%% 求 (i , j) 和 (i+1 , j) 之间的平均导数； 
bar_partial_F1_av = (partial_F1 + bar_partial_F1) / 2;
bar_partial_F2_av = (partial_F2 + bar_partial_F2) / 2;
bar_partial_F3_av = (partial_F3 + bar_partial_F3) / 2;
bar_partial_F4_av = (partial_F4 + bar_partial_F4) / 2;

%% 计算校正步的人工粘性 bar_SF；
bar_SF1 = zeros(1,J);
bar_SF2 = zeros(1,J);
bar_SF3 = zeros(1,J);
bar_SF4 = zeros(1,J);

for j = 2 : J-1
    bar_SF1(j) = Cy*abs(bar_p(j+1) - 2*bar_p(j) + bar_p(j-1))/(bar_p(j+1) + 2*bar_p(j) + bar_p(j-1))...
              * (bar_F1(j+1) - 2*bar_F1(j) + bar_F1(j-1));
          
    bar_SF2(j) = Cy*abs(bar_p(j+1) - 2*bar_p(j) + bar_p(j-1))/(bar_p(j+1) + 2*bar_p(j) + bar_p(j-1))...
              * (bar_F2(j+1) - 2*bar_F2(j) + bar_F2(j-1));
          
    bar_SF3(j) = Cy*abs(bar_p(j+1) - 2*bar_p(j) + bar_p(j-1))/(bar_p(j+1) + 2*bar_p(j) + bar_p(j-1))...
              * (bar_F3(j+1) - 2*bar_F3(j) + bar_F3(j-1));
          
    bar_SF4(j) = Cy*abs(bar_p(j+1) - 2*bar_p(j) + bar_p(j-1))/(bar_p(j+1) + 2*bar_p(j) + bar_p(j-1))...
              * (bar_F4(j+1) - 2*bar_F4(j) + bar_F4(j-1));    
end

%% (i+1 , j) 处的校正值；
F1_t = F1 + bar_partial_F1_av * Delta_xi + bar_SF1;
F2_t = F2 + bar_partial_F2_av * Delta_xi + bar_SF2;
F3_t = F3 + bar_partial_F3_av * Delta_xi + bar_SF3;
F4_t = F4 + bar_partial_F4_av * Delta_xi + bar_SF4;

%% 边界条件：计算(i+1,41)上边界处的流场参数；
% 注：此处采用沿垂直方向线性外插的方式是不合适的，而应该从内部网格点沿特征线外插；
F1_t(J) = 2*F1_t(J-1) - F1_t(J-2);
F2_t(J) = 2*F2_t(J-1) - F2_t(J-2);
F3_t(J) = 2*F3_t(J-1) - F3_t(J-2);
F4_t(J) = 2*F4_t(J-1) - F4_t(J-2);

%% (i + 1, j)处的原变量计算； 
% 中间变量 A、B 和 C 的计算；
A_t = F3_t.^2 ./ (2 * F1_t) - F4_t;
B_t = gama/(gama - 1) * F1_t .* F2_t;
C_t = -(gama + 1) / (2*(gama - 1)) * F1_t.^3;

% 原变量的计算；
rho_t = (-B_t + sqrt(B_t.^2 - 4*A_t.*C_t)) ./ (2*A_t);
u_t = F1_t ./ rho_t;
v_t = F3_t ./ F1_t;
p_t = F2_t - F1_t.*u_t;
T_t = p_t ./ (rho_t * Rg);
Ma_t = sqrt((u_t.^2 + v_t.^2) ./ (gama * Rg * T_t));

%% 边界条件：计算(i+1,1)下边界处的流场参数；
psi = atan(abs(v_t(1))/u_t(1));% 速度向量的夹角；
% 分为扩张前(x<=E)的壁面和扩张后(x>E)的壁面的计算；
phi = psi*(x<=E) + (theta - psi)*(x>E);

f_cal = sqrt((gama+1)/(gama-1))*atan(sqrt((gama-1)/(gama+1)*(Ma_t(1)^2-1)))...
        - atan(sqrt(Ma_t(1)^2-1));
f_act = f_cal + phi; % 普朗特-迈耶关系式；
% 由f_act求解Ma_act;
Fun = @(Ma)sqrt((gama+1)/(gama-1))*atan(sqrt((gama-1)/(gama+1)*(Ma^2-1)))...
       - atan(sqrt(Ma^2-1)) - f_act;
Ma_act = fsolve(Fun, 2.2);

% 计算真实的压力、密度与温度；
% Ma_t(1) 即 Ma_cal;
p_t(1) = p_t(1) * ((1 + ((gama-1)/2) * Ma_t(1)^2) / (1 + ((gama-1)/2) * Ma_act^2))^(gama/(gama-1));
T_t(1) = T_t(1) * ((1 + ((gama-1)/2) * Ma_t(1)^2) / (1 + ((gama-1)/2) * Ma_act^2));
rho_t(1) = p_t(1)/(Rg*T_t(1));

% 计算 y 方向的速度 v;
% u_t(1) = u_cal;
v_t(1) = 0 * (x<=E) + (- u_t(1) * tan(theta))*(x>E);
Ma_t(1) = Ma_act; % Ma-act才是真正的马赫数；

%% 从(i+1,j)到(i+2,j)的推进步长 Delta_xi；
% 计算马赫角 mu;
mu = asin(1 ./ Ma_t);
% 计算推进步长 Delta_xi；
tan_max = max(max(abs(tan(theta + mu))), max(abs(tan(theta - mu))));
Delta_xi = CFL * Delta_y / tan_max;
end

test_Prandtl_Meyer.m

% Author: CXT
% Date: 2021/5/6
% Theme: 二维超声速流动的数值解——普朗特-迈耶稀疏波；
clc;
clear;
close;

%% 初始条件(x = 0)；
J = 41;
u = 678 * ones(1,J); % [m/s];
v = zeros(1,J);% [m/s];
rho = 1.23 * ones(1,J); % [kg/m3];
p = 0.101 * 10^6 * ones(1,J); % [N/m2];
T = 286 * ones(1,J); % [K]
Ma = 2 * ones(1,J); % [1];

%% test;
x = 0;
for i = 1:100
    [u, v, rho, p, T, Ma, Delta_xi] = Prandtl_Meyer(x, u, v, rho, p, Ma);  
    x = x + Delta_xi;
    if (x>=12.928)
        break;
    end
end
%% 绘图；
% 计算 x 所在的 y 的分布；
x = x - Delta_xi;
theta = 5.352/180 * pi;% 转角大小；
E = 10;% 转角位置；
ys = @(x) 0*(x>=0 & x<=E) + (-(x - E)*tan(theta))*(x>=E);% 下边界的纵坐标；
% h —— 物理平面从下表面到上边界的距离；
h = @(x) 40 * (x>=0 & x<= 10) + (40 + (x - 10)*tan(theta))*(x>=10 & x<=45);
Delta_y = h(x)/(J - 1);
y = ys(x):Delta_y:40;

figure;
plot(u,y,'*-');
xlabel('u/(m/s)');
ylabel('y/m');
title("x = 12.6905 m 处中心稀疏波超声速流的结果");

额外说明

在计算到 x = 45.0246 m 处，程序就爆掉了。错误提示信息如下：

索引超出数组元素的数目(0)。

出错 Prandtl_Meyer (第 70 行)
    partial_F1(j) = partial_eta_x(x,y(j))*(F1(j) - F1(j+1))/Delta_eta...

出错 test_Prandtl_Meyer (第 20 行)
    [u, v, rho, p, T, Ma, Delta_xi] = Prandtl_Meyer(x, u, v, rho, p, Ma);

笔者未找出原因，也未对该问题进行深究！

西西西西柚
关注
专栏目录
二维超声速流——普朗特-迈耶稀疏波的流场CFD解（附完整代码）
weixin_45205765的博客
08-12 3940
入门CFD，主要参考书目《计算流体力学基础及其应用》（John D.Anderson著，吴颂平等译） 实现了第 8.3节普朗特-迈耶稀疏波流场的数值解 的代码，采用的是MacCormack 方法，对守恒型方程求解。 注意： 代码最后运行的结果与作者在书中所提供的并不完全相同： 1、当采用作者在书中提供的思路，即在流场下边界壁面处均采用向前差分计算时，结果是在壁面处速度大于稀疏波后速度，稀疏波后速度约710m/s，近壁面速度可达720m/s以上，而作者在书中指出近壁面速度约为706-70...
基于matlab求Prandtl-Meyer膨胀波流场数值解
qq_59747472的博客
01-26 122
含有水蒸气的气流,经过快速膨胀,温度降低,水蒸气凝结放热,影响流场结构.这一点在Prandtl-Meyer膨胀波中尤为明显.为此,通过流场显示的方法,观察了不同来流湿度下,后台阶流所形成膨胀波中,凝结对流场结构的影响.发现随着来流湿度的增加,凝结区域中会产生一道由弱到强的凝结激波.而湿度过大时,又会在上游产生一道正激波,反而抑制了下游的凝结.同时通过数值模拟,与实验研究进行比较验证,结果基本一致.数值模拟还发现,一定来流条件下,凝 结流场出现自激振荡现象.
探秘 Supersonic：一款超音速的 Python 性能优化库
最新发布
gitblog_00079的博客
04-18 546
探秘 Supersonic：一款超音速的 Python 性能优化库 supersonicA lightweight and full-featured cross-platform desktop client for self-hosted music servers项目地址:https://gitcode.com/gh_mirrors/sup/supersonic 项目简介 是一个开源的 P...
基于MATLAB在流体力学中的数值分析.pdf
07-10
基于MATLAB在流体力学中的数值分析.pdf
普朗特迈耶稀疏波_埃里克·迈耶的访谈
culuo8053的博客
08-02 739
普朗特迈耶稀疏波I've always wanted to interview Eric Meyer. His early CSS books are a big reason this blog exists today and the reason why I'm a web developer. Eric gave me some time to hit the history of C...
查询所有Redis集群配置信息-data-2021-05-06 11_08_34.zip
05-06
查询所有Redis集群配置信息-data-2021-05-06 11_08_34.zip
查询所有Redis集群配置信息-data-2021-05-06 09_56_40.zip
05-06
查询所有Redis集群配置信息-data-2021-05-06 09_56_40.zip
Video_2021-06-28_103102.wmv
06-28
Video_2021-06-28_103102.wmv
二维超声流 普朗特迈耶稀疏波
12-28
二维超声流 普朗特迈耶稀疏波
亚声速—超声速等熵喷管流动的CFD解法
08-28
安德森《计算流体力学及其应用》第七章拟一维喷管流动的数值解 7.3节 亚声速—超声速等熵喷管流动的CFD解法 。
流过平面的超声速流动数值模拟
05-15
二维完全纳维斯托克斯方程的显显式有限差分解法，对流过平板的超声速流动进行数值模拟，Fortran语言
Matlab 编程 《计算流体力学基础及应用（约翰D安德森）》 亚声速-超声速等熵喷管流动CFD解法 拉瓦尔喷管 守恒形式方程解法
马儿不吃草
12-17 5921
《计算流体力学基础及应用（约翰D安德森）》 亚声速-超声速等熵喷管流动CFD解法，拉瓦尔喷管 守恒形式方程解法 Matlab 编程 问题之 亚声速一超声速等熵喷管流动的CFD解法初始化参数迭代过程绘图结果 最近学习经典的计算流体力学入门书籍《计算流体力学基础及应用（约翰D安德森）》，跟着作者给出的过程写了计算程序，程序写的比较初级。 问题之 亚声速一超声速等熵喷管流动的CFD解法 问题的提法见《...
声音的本质 The Nature of Sound
qq_38149276的博客
03-29 734
引言声由振动产生，所以我们将从物理学开始讲起。压力波声音通常是由物体的循环运动发出的，比如：扬声器的膜片不断地跳动，声带之间的间隙时而缩小时而扩大，或者吉他弦来回振动。正是这些运动向周围介质(固体、液体或气体)的传播及其在介质中的传播过程构成了声音。在某些情况下，运动的就是介质本身，比如有人在瓶口上方吹气时瓶颈里的空气就是如此。非运动源包括突然释放的热能(如爆炸或火花)和快速振荡的热源。当扬声器的膜片运动时，膜片产生的声波具有和电信号相同的变化规律。每当膜片向外移动时，它同步挤压前面的空气分子，使它们彼此靠
音频基础学习一——声音的本质、术语与特性
山河君的分享博客
08-01 4036
百度词条中介绍声音的概念声音（sound)是由物体振动产生的声波。是通过介质（空气或固体、液体）传播并能被人或动物听觉器官所感知的波动现象。最初发出振动（震动）的物体叫声源。声音以波的形式振动（震动）传播。声音是声波通过任何介质传播形成的运动。声音是一种波。可以被人耳识别的声（频率在20Hz~20000Hz之间），我们称之为声音。本篇文章以最根本的物理角度来看待声音，其中的内容都是笔者自己学习整理的。如果有错误或者需要补充的地方希望大家多多指正。提示这里对文章进行总结例如https。............
流过平板的超声速流动_高超声速空气动力学 | 基本特征
weixin_39824834的博客
11-26 1675
内容结构指引学习路线 | 高超声速流动特征 | 高超声速流基本激波关系式 | 高超声速流基本膨胀波关系式主要名词检索高超声速 | 对流层 | 平流层 | 同温层 | 电离层 | 位势高度 | 气动力 | 动压 | 声速 | 马赫数 | 压力系数 | 升力系数 | 阻力系数 | 法向力系数 | 轴向力系数 | 力矩系数 | 特征长度 | 小密度比 | 薄激波层 | 粘性效应 | 高熵层 | 涡干扰 ...
Deep Learning论文笔记之（二）Sparse Filtering稀疏滤波
热门推荐
zouxy09的专栏
08-15 4万+
Deep Learning论文笔记之（二）Sparse Filtering稀疏滤波zouxy09@qq.comhttp://blog.csdn.net/zouxy09          自己平时看了一些论文，但老感觉看完过后就会慢慢的淡忘，某一天重新拾起来的时候又好像没有看过一样。所以想习惯地把一些感觉有用的论文中的知识点总结整理一下，一方面在整理过程中，自己的理解也会更深，另一方面也方便未来自己
2021-05-08
ErnestChai的博客
05-08 1332
不可压库埃特流的数值解 —— 隐式克兰克-尼科尔森（Crank-Nicolson）方法 本文主要分享了安德森的《计算流体力学基础及其应用》一书中：“不可压库埃特 ( Couette ) 流的数值解” 的程序求解。水平有限，请见谅！ 文章目录不可压库埃特流的数值解 —— 隐式克兰克-尼科尔森（Crank-Nicolson）方法最终结果绘图Matlab 程序testCouette.mCouette.m 最终结果绘图 Matlab 程序 testCouette.m %% % Author: 柴星涛 % Dat
java 集合 2023-01,2023-02,2023-03,2023-04 2023-05,2023-06, 2022-01,2022-02,2022-03,2022-04,2022-05,2022-06，2021-01,2021-02,2021-03,2021-04,2021-05,2021-06 这样的数据变成 2023-01,2022-01,2021-01,2023-02,2022-02,2021-02
06-09
您可以使用Java中的集合类（如ArrayList）对日期进行排序，然后遍历并提取年份和月份，将其组合成新的日期字符串。以下是一个示例代码： ```java import java.util.ArrayList; import java.util.Collections; import java.util.List; public class DateSorter { public static void main(String[] args) { List<String> dates = new ArrayList<>(); dates.add("2023-01"); dates.add("2023-02"); dates.add("2023-03"); dates.add("2023-04"); dates.add("2023-05"); dates.add("2023-06"); dates.add("2022-01"); dates.add("2022-02"); dates.add("2022-03"); dates.add("2022-04"); dates.add("2022-05"); dates.add("2022-06"); dates.add("2021-01"); dates.add("2021-02"); dates.add("2021-03"); dates.add("2021-04"); dates.add("2021-05"); dates.add("2021-06"); // 对日期进行排序 Collections.sort(dates); // 遍历并提取年份和月份，组合成新的日期字符串 List<String> newDates = new ArrayList<>(); String prevMonth = ""; for (String date : dates) { String[] parts = date.split("-"); String year = parts[0]; String month = parts[1]; if (!month.equals(prevMonth)) { newDates.add(date); } prevMonth = month; } // 输出结果 for (String date : newDates) { System.out.println(date); } } } ``` 输出结果为： ``` 2021-01 2022-01 2023-01 2021-02 2022-02 2023-02 ```
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值