• 研究基于印度NIFTY 50指数49只股票的VWAP（成交量加权平均价）时间序列数据，采样时间窗口为近1000个交易日，确保时间序列的平稳性和多谱估计稳定性 [page::5][page::6]。

- 结合谱均值和二谱均值计算多种加权函数下的高阶谱特征，捕捉股票时间序列中的非线性和高阶相互作用，这些特征优于传统的一阶分析指标 [page::8]。

• 图示展示了不同股票的自相关函数（ACF）分布，验证无明显季节性存在，支持采用差分和多谱分析进行特征提取 [page::10]。

• 多组仿真测试验证了该聚类方法在不同趋势、分布及非线性条件下的高分类准确率，最高达100%，展示了谱和二谱均值特征的强区分能力。

• 使用Hopkins指数(0.7997, p<0.001)和多种指标（Elbow、Silhouette、Gap Statistic等）判断数据适合聚类，结合Gap统计量推荐k=5为最佳聚类数。

• 聚类分析结果展示不同k值下的聚类形态，确认k=5时划分合理且聚类结构明显 [page::18]

• 多种聚类算法（k-means、PAM、CLARA、FANNY）应用，对聚类结果进行对比，k-means算法形成的五类簇具有明确的行业及市场逻辑，分别涵盖大型财团与国有行、能源与银行混合组、IT与公共事业、防御性行业、以及多行业大盘股 [page::20][page::21]。

• 特征相关性分析显示二谱均值与价格差异特征（DiffEndStart, MeanDiff, MaxDiff）显著相关，表明非线性时序属性是划分的重要依据，股票聚类主要基于价格行为差异而非传统行业界定。

• 五个聚类的实地金融意义阐释：

- 簇1：大型工业集团与基础设施，体现基建与大宗商品周期关联，还有国企行SBI的金融链接。
- 簇2：能源、银行、汽车及消费混合，反映公共政策、利率及经济周期的综合影响。
- 簇3：IT龙头和电力公用事业，表现出稳定的营收与低波动性特征。
- 簇4：HCL科技与太阳药业，因国际业务和机构投资者偏好形成防御性组。
- 簇5：多行业大盘蓝筹组合，涵盖周期性与防御性，展现强流动性及机构持有特征。

• 本文重视量化特征及非线性结构挖掘，突破传统行业分类限制，提升对印度股市复杂动态的理解 [page::23][page::24][page::27][page::29]。

极其详尽与全面分析报告

Polyspectral Mean based Time Series Clustering of Indian Stock Market

作者：Dhrubajyoti Ghosh
机构：Duke University生物统计与生物信息学系
相关邮箱：dg302@duke.edu

---

1. 元数据与概览

报告标题：《Polyspectral Mean based Time Series Clustering of Indian Stock Market》
作者及机构：Dhrubajyoti Ghosh，Duke University
文献类型：学术研究论文 / 方法学应用研究
研究对象：印度股市中的49只股票
核心主题：基于高阶频谱特征（多谱平均）利用k-means算法对印度股票时间序列数据进行聚类分析
核心创新：采用传统一阶频谱无法捕获的高阶谱（尤其是双谱）特征，通过不同加权函数计算多谱平均，用以揭示时间序列中的非线性趋势与结构，实现股票的时间序列聚类
核心结论：

• 识别出5个独特的股票簇，分布上不仅涵盖行业聚类，还表现出集团、所有权结构（例如国有银行的非典型归属）等复杂的市场内在关系。

- 特别地，ADANI、BIRLA、TATA集团与国有银行SBI归入一同一簇；IT板块(TCS、WIPRO)归为一组；第三组包含私人银行、国有企业和RELIANCE；最后一个簇则包含分散所有权的上市公司。

• 该聚类揭示了印度股市中复杂金融关系，能为投资者和分析师提供深刻见解。

关键词：多谱（Polyspectra）、双谱（Bispectra）、聚类（Clustering）、Gap统计方法
[page::0]

---

2. 逐节深度解读

2.1 引言(Introduction)

关键论点：

• 印度股市作为国家金融体系关键部分，具有悠久历史和高度复杂的结构。

- 股市聚类是指股票基于相似特征或表现出价格共动形成的群组现象，这是全球市场普遍存在的。

• 聚类依据多样，包括行业板块（sectoral clustering）、市值大小（大盘股、中盘股、小盘股）、风险特征和地理因素。

- 行业间共同经济驱动因素使行业内部股票价格具有同步性，而市值聚类反映不同的风险-收益偏好。

• 识别聚类对投资组合管理、风险分散与投资决策至关重要。

作者推理依据：

研究通过概述已发表的文献引入聚类现象，提及印度具体市场特征与聚类存在的多维度，包括部门、市值、风险和地域。聚类的经济学解释结合行业经济驱动力，市场周期与股价动态，强调对免疫投资组合的重要性。
[page::1,2]

---

2.2 非线性时间序列聚类的必要性和方法(Introduction继续)

关键论点：

• 传统聚类多依赖于时间序列的“一阶特征”（如基于谱密度或ARMA模型系数），在处理非线性时间序列时信息缺失较多。

- 印度股市显示非线性行为，已有文献证明这种非线性存在且对回报预测有影响。

• 高阶统计量（多谱）特别是双谱能揭示非线性关系及高阶相互作用，提供比功率谱更多信息。

- 本研究首次尝试用多谱均值（spectral与bispectral mean）特征进行股票时间序列聚类。

• 这种基于时间序列的聚类对构建投资组合、风险评估、行业轮动策略以及行为金融学分析具有实际应用价值。

逻辑与证据：

• 作者从非线性时间序列的理论和实证研究谈起，指出铁路目前相关方法缺失，采用Ghosh等人2024年关于多谱均值的理论成果给出方法支撑。

- 明确给出双谱公式及其加权平均的表达式，指出高阶多谱处理非高斯非线性数据的优势。

• 作者强调目前无相关文献使用多谱均值做时间序列聚类，体现该方法的创新性。

[page::3,4]

---

2.3 数据与方法(Data and Methods)

数据描述：

• 选取NIFTY 50中除Bharti Infratel（巿场结构变动导致数据不可用）外的49只股票，时间跨度2000年1月1日到2021年4月30日。

- 选用最近1000天的VWAP做时间序列分析对象，权衡长期市场结构和短期噪音。

• VWAP定义明晰，强调它融合价格和成交量，为市场行为提供精准刻画，有利高阶谱特征提取。

特征提取与数学定义：

• 对股票差分序列进行多谱（k阶谱）均值计算，k=1时为谱均值，k=2时为双谱均值。

- 采用各类加权函数（如指示函数、三角窗、圆盘函数、余弦等）对频域信息加权，提取不同频段和相位依赖特征。

• 利用离散傅里叶变换，结合\(\Phi(\lambda)\)指标函数过滤频率组合，保证统计量合理性和无偏性。

- 差分操作确保序列平稳，满足多谱均值理论分布假设。

数据预处理与分析细节：

• 说明无显著季节性（通过ACF图像，缺少周期峰值），故无明显季节调整需求。

- 说明噪音对多谱估计影响有限，数据处理因为差分和平滑窗口长度（1000天）而有效减弱短期市场噪音干扰。

关键数据点：

• 表1揭示具体股票的统计特征，如均值收益、波动率、一阶自相关、赫斯特指数（揭示趋势或均值回归性质）、趋势强度、季节强度等，体现股票时间序列内生特征的多样性。

[page::5,6,7,8,9,10]

---

2.4 仿真实验(Simulation)

目的

评估k-means聚类算法基于多谱均值的特征组合能否成功区分不同类型的时间序列。

仿真设计

• Simulation I：两类平稳时间序列，分别生成ARMA(2,2)过程和GARCH(1,1)过程，加上高斯噪声，100长度；

- Simulation II：在Simulation I基础上增加不同时间趋势，使两类更明显区分；

• Simulation III：三类，分别为常数项加ARMA、二次函数加ARMA和二次函数倍乘正弦的GARCH过程。

结果

• 多次仿真显示k-means对不同类别有较高分类准确率，灵敏度与特异性普遍超过80%，部分达到100%；

- 仿真2分类准确率均为100%，说明趋势信息被有效捕获；

• 仿真3三类准确率依旧在高水平（F1分数0.89-0.93，平均AUC约为0.95）；

- 图表（Fig3,4,5,6）直观展示时间序列与特征重要性排名，验证谱均值和双谱均值特征在分类中的关键作用。

[page::11,12,13,14,15]

---

2.5 数据实证分析(Real Data Analysis)

聚类潜力验证

• Hopkins指数为0.7997，p<0.001，显示数据极具聚类性。

- 图7对比市场特征矩阵和随机数据的异质性矩阵，市场数据展示明显模式，随机数据无此特征。

聚类数目选择

• 采用多种指标（Elbow、Silhouette、Gap Statistic、Calinski-Harabasz、Davies-Bouldin、Dunn指数），

- 结果显示2-3簇过少，无明显Elbow，Gap Statistic提示可选择较大k，结合作者建议综合权衡，k=5为合理聚类数。

• 图8-9展示不同k的聚类形态，对比发现k>5簇过细，簇内样本数过少且意义不大，k=5在平衡内聚力和解释性上效果最佳。

聚类算法比较

• 采用k-means, PAM（基于中位数点的鲁棒聚类法）、CLARA（大数据子样本PAM）、FANNY（模糊聚类）四种算法。

- FANNY仅获得3簇，其他三种均获得5簇，K-means结果最直观，符合预期。

• 聚类结构表现稳定，算法表现支持所选聚类数和结果。

聚类结果解读

• Cluster 1：ADANI, BIRLA, TATA集团核心公司和SBI国有银行，共享基础设施建设、商品循环和政策依赖性。

- Cluster 2：跨能源、银行、汽车、消费品和基础设施，包含众多政府控股企业和私企，反映政策敏感性及经济周期影响。

• Cluster 3：TCS、WIPRO、POWERGRID，稳定收益、低波动性公司，机构投资者青睐的防御性股票。

- Cluster 4：HCL科技与Sun Pharma，全球市场布局相似及机构投资者偏好，非行业同类但表现相近。

• Cluster 5：多行业大盘股核心，高流动性和机构持股，结合周期性与防御性的投资特征。

统计特征驱动聚类的角色

• 图11显示双谱均值和VWAP差分特征（起终差、平均差、最大差）在聚类中占主导作用。

- 这些特征体现股票价格的非线性依赖与价格波动特征，不单纯依赖行业分类。

• 多维相关特征表明多谱均值与价格变动指标的中度正相关，说明非线性动态与趋势紧密相关。

[page::15,16,17,18,19,20,21,22,23]

---

2.6 讨论(Discussion)

• 聚类反映市场动态、机构投资以及宏观经济环境的共同作用，非单纯行业关系。

- 股价的高阶谱特征捕捉了传统财务指标难以体现的非线性、隐秘依赖关系。

• 各聚类板块背后隐含着资本流动、政策敏感和结构性经济联系。

- 聚类结果内涵复杂，尚需进一步结合高频数据、其他特征挖掘方式和不同聚类模型来验证与拓展。

• 研究所示聚类与传统行业分类部分重合，但更强调市场动态的交互作用。

[page::24,25,26,27,28]

---

2.7 结论(Conclusion)

• 利用k-means和PAM算法，结合多谱均值（谱与双谱）特征对印度49支股票进行聚类，揭示股票价格行为中的非线性关系。

- 共发现5个聚类，反映了市场驱动而非纯行业驱动的群组结构。

• 各聚类分别对应大工业集团、政府与私企混合板块、IT与公用事业、防御性创新科技和多元化大盘龙头。

- 结果证实高阶频谱分析能更深入地揭露市场隐含结构及投资者行为，补足传统分析不足。

• 建议今后研究进一步探索多样特征提取、不同聚类技术和高频数据的应用，以完善股票市场结构的理解和风险管理。

[page::29]

---

3. 图表深度解读

3.1 表格

表1：VWAP缩放序列统计量（时间序列主要特征描述）

• 统计指标：均值收益率（Mean Return），波动率（Volatility），一阶自相关（ACF Lag-1），赫斯特指数（Hurst Exp），趋势强度（Trend Strength），季节强度（Seasonal Strength）。

- 如ADANIPORTS呈现较高的均值收益和趋势强度，且Hurst指数>0.5，显示趋势显著；波动率中等，季节成分弱。

• 不同股票的数据体现其行业属性与股价动态差异，为谱特征提取提供统计基础。

- 该表体现每只股票的结构及其稳定性，为后续聚类提供量化基础指标。
[page::7]

3.2 图1：谱均值与双谱均值权重函数

• 左列列出了谱均值权重函数（如恒定1、指标函数、三角窗等）；右列列出了双谱均值中常用权重（如圆盘指标函数、径向函数、余弦三角组合等）。

- 这些权重逻辑是提取对应频率范围和交互成分，重点捕获股价的非线性复杂关系。

• 权重组合有助于增强多谱均值特征的稳定性和解释性。

[page::9]

3.3 图2：ADANIPORTS和ASIANPAINT的ACF图

• 展示两支样本股票差分后的自相关函数，横轴为滞后期数，纵轴为自相关值。

- 图中无明显周期峰值，说明短期内无显著季节性循环。

• 支撑作者关于不需要显式季节调整的论断。

[page::10]

3.4 图3-6：仿真时间序列及特征重要性

• 图3-5左侧显示多类仿真时间序列变化趋势，右侧为对应特征在聚类中重要性排序。

- 在仿真1和2中特征如谱均值、双谱均值和最大差异依赖度高。

• 仿真3的准确率优异，显示该方法有效区分多类复杂趋势和过程区别。

- 图6呈现三组仿真数据在二维降维空间的聚类效果，形态明显区分。
[page::12,13,14,15]

3.5 图7：真实数据与随机数据异质性矩阵对比

• 左图为股票市场特征矩阵距离（色值表示差异度），复杂结构明显，显示数据内在群落特征；

- 右图为随机数据，结构无规无序，边界明显差异。

• 突出股票数据存在内在可聚类结构，为聚类分析提供依据。

[page::16]

3.6 图8：多种聚类数判定指标图

• 包含Elbow、Silhouette平均值、Gap统计、Calinski-Harabasz、Davies-Bouldin、Dunn指数。

- Elbow法不明显，Silhouette和Gap呈上升趋势，Dunn指数趋于2，Calinski-Harabasz认为5较合适。

• Gap方法依据Tibshirani规则选出k=5聚类数合理。

[page::17]

3.7 图9：k=2至9时聚类二维表示

• 图示K值变化时簇的可分性和稳定性，k=5时簇分布较合理，无过度碎片化。

- k=2或3簇过少，k>=6簇过细，聚类意义减弱。
[page::18]

3.8 图10：k=5时四种聚类算法结果对比

• k-means、PAM、CLARA显示5个清晰簇，FANNY显示3个簇。

- K-means聚类分配符合财团、行业及股权结构基础。

• 不同算法对簇的稳健性验证，显示核心聚类结构稳定。

[page::21]

3.9 图11：关键特征间散点、密度、相关性分析

• 对角线为各特征跨簇分布密度图，说明各特征对聚类的区分度；

- 散点图显示特征间的分布及簇内关系；

• 相关性矩阵显示双谱均值与VWAP差分类特征正相关，突出非线性动态与趋势的联动。

- 说明聚类基于更复杂的价格动态，而非仅传统经济分类。
[page::23]

---

4. 估值分析与风险因素

本报告无传统财务估值（如DCF、PE、EV/EBITDA）的讨论，聚焦于时间序列特征空间的无监督学习聚类框架，不涉及具体估值模型或价格预测，风险分析亦未详述，未来可结合此聚类探索风险敞口及波动性分布。

---

5. 风险因素评估

报告未专门列出风险因素条目，但可从方法角度挖掘潜在风险：

• 时间序列窗口长度选择（当前1,000天），未系统验证对结果的敏感性，可能影响稳定性。

- 只使用VWAP差分价，未考虑成交量、宏观经济数据联合建模信息可能丢失。

• 聚类模型（k-means、PAM等）对初始参数敏感，结果存在一定波动，未详细讨论算法选择风险。

- 复杂金融市场行为及结构变动可能导致聚类结构的时变与漂移。

• 数据处理未显式滤除高频噪声，虽双谱抗噪性强，但可能仍受极端事件影响。

未来应结合更多稳定性和敏感性分析。

---

6. 批判性视角与细微差别

• 作者强调多谱均值在发现非线性时序依赖方面的优势，但对比基准（如其他高阶谱方法或机器学习方法）不足；缺乏与经典因子模型对比验证。

- 聚类结果解释虽贴合行业逻辑，但存在部分“意外”分类（如政府银行归于工业集团簇），未深入探讨异类样本的形成机制。

• 选用聚类算法多，文本对结果差异讨论较弱，FANNY聚类少簇但未详解原因。

- 功能权重函数选择多样，但权重对聚类结果的具体影响未展开敏感性分析。

• 股票市场的高频交易、非对称信息可能导致非平稳，基于平稳假设的谱方法或受限。

- 数据时间窗未变更实验，影响聚类结构稳定性的研究缺失。

• 总体研究创新性强，方法解释清晰，实验支持充分，应用说明丰富。

---

7. 结论性综合

本研究基于多谱均值提取股票价格时间序列中非线性结构特征，结合k-means及其他聚类算法，首次实证展示印度NIFTY 50股票的五大内在聚类结构，揭示市场动态下的复杂关联。这些聚类不仅对应传统行业，更体现股权结构、机构持股和宏观政策脉动带来的价格行为共性。时间序列的双谱与差分特征在区分股票分布中发挥关键作用，反映投资者情绪与非线性市场机制深刻影响股票走势。本研究为基于频域的复杂金融时间序列分析提供新视角，助力投资组合优化、风险评估与市场结构理解。图表直观支持理论推断，尤其统计描述表与特征重要性图反映数据本质。未来研究应进一步结合高频数据，拓展模型及特征空间，提升聚类稳定性与应用深度。

总体来看，作者利用创新高阶谱方法，有效挖掘出印度股票市场的时间序列聚类特征，兼顾理论严谨与实证有效，具有较高学术及实务价值。

[page::0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]

---

附：主要图表展示

图2：ADANIPORTS与ASIANPAINT的差分序列自相关图


图3：仿真1时间序列样本与特征重要性


图4：仿真2时间序列样本与特征重要性


图5：仿真3三组时间序列曲线


图6：仿真3聚类二维表示


图7：股票市场与随机数据异质性矩阵对比


图8：多种指标确定聚类簇数


图9：k=2至9的聚类分布比较


图10：k=5时四种聚类算法结果


图11：关键特征间散点、密度及相关性

报告