{"id":1749,"date":"2023-08-10T14:14:04","date_gmt":"2023-08-10T06:14:04","guid":{"rendered":"https:\/\/www.huangrongzhen.ink\/?p=1749"},"modified":"2023-08-12T17:32:20","modified_gmt":"2023-08-12T09:32:20","slug":"c-%e8%af%ad%e8%a8%80%e7%94%9f%e6%88%90%e9%9a%8f%e6%9c%ba%e6%95%b0","status":"publish","type":"post","link":"https:\/\/www.huangrongzhen.ink\/?p=1749","title":{"rendered":"C \u8bed\u8a00\u751f\u6210\u968f\u673a\u6570"},"content":{"rendered":"<div class=\"wp-block-post-excerpt\"><p class=\"wp-block-post-excerpt__excerpt\">\u7b97\u6cd5\u7b14\u8bb0 <\/p><\/div>\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E5%8C%85%E5%90%AB%E7%9A%84%E5%A4%B4%E6%96%87%E4%BB%B6\"><\/span>\u5305\u542b\u7684\u5934\u6587\u4ef6<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">#include \"math.h\"\n#include \"malloc.h\"<\/code><\/pre>\n\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E5%9D%87%E5%8C%80%E5%88%86%E5%B8%83%E7%9A%84%E9%9A%8F%E6%9C%BA%E6%95%B0\"><\/span>\u5747\u5300\u5206\u5e03\u7684\u968f\u673a\u6570<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Uniform\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u5747\u5300\u5206\u5e03\u7684\u968f\u673a\u6570\n* \u8f93\u5165\u53c2\u6570\uff1a a\uff1a\u533a\u95f4\u4e0b\u9650\n*            b\uff1a\u533a\u95f4\u4e0a\u9650\n*            seed\uff1a\u968f\u673a\u6570\u79cd\u5b50\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a \u968f\u673a\u6570\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\ndouble Uniform(double a, double b, long int* seed)\n{\n  double t;\n  *seed = 2045 * (*seed) + 1;\n  *seed = *seed - (*seed \/ 1048576) * 1048576;\n  t = (*seed) \/ 1048576.0;\n  t = a + (b - a) * t;\n  return t;\n}<\/code><\/pre>\n\n\n\n<p>\u5728 -50~50 \u533a\u95f4\u5185\uff0c\u968f\u673a\u6570\u79cd\u5b50\u4e3a 13579 \u7684\u6761\u4ef6\u4e0b\uff0c\u8fde\u7eed\u83b7\u53d6\u4e00\u5343\u4e07\u4e2a\u968f\u673a\u6570\uff0c\u6d4b\u91cf\u7ed3\u679c\u5982\u4e0b\u3002\u7b2c\u4e00\u6761\u66f2\u7ebf\u4e3a\u6982\u7387\u5bc6\u5ea6\u51fd\u6570\u6240\u5bf9\u5e94\u7684\u66f2\u7ebf\u3002\u7b2c\u4e8c\u6761\u66f2\u7ebf\u8868\u793a\u5728 -50~50 \u8fd9\u4e2a\u533a\u95f4\u5185\uff0c\u53d6 1000 \u7b49\u4efd\uff0c\u7edf\u8ba1\u843d\u5230\u8fd9\u4e2a\u533a\u95f4\u5185\u7684\u6d4b\u91cf\u503c\u7684\u6570\u91cf\u3002\u4ece\u7b2c\u4e8c\u6761\u66f2\u7ebf\u53ef\u4ee5\u770b\u51fa\uff0c\u843d\u5230\u6bcf\u4e2a\u533a\u95f4\u5185\u7684\u6d4b\u91cf\u503c\u7684\u6570\u91cf\u5f88\u5747\u5300\uff0c\u5927\u81f4\u5728 9960~10040 \u8fd9\u4e2a\u8303\u56f4\u4e4b\u5185\uff0c\u5e73\u5747\u503c\u4e3a 10000\u3002\u5e73\u5747\u503c\u662f 10000\uff0c\u4e58\u4e0a\u533a\u95f4\u6570\u91cf 1000\uff0c\u6b63\u597d\u662f\u7b49\u4e8e\u6d4b\u91cf\u6b21\u6570\u4e00\u5343\u4e07\u3002<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-large\"><img decoding=\"async\" loading=\"lazy\" width=\"1024\" height=\"440\" src=\"http:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u5747\u5300\u5206\u5e03\u7684\u968f\u673a\u6570\u6d4b\u91cf\u7ed3\u679c-20230812B-1024x440.png\" alt=\"\" class=\"wp-image-1781\" srcset=\"https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u5747\u5300\u5206\u5e03\u7684\u968f\u673a\u6570\u6d4b\u91cf\u7ed3\u679c-20230812B-1024x440.png 1024w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u5747\u5300\u5206\u5e03\u7684\u968f\u673a\u6570\u6d4b\u91cf\u7ed3\u679c-20230812B-300x129.png 300w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u5747\u5300\u5206\u5e03\u7684\u968f\u673a\u6570\u6d4b\u91cf\u7ed3\u679c-20230812B-768x330.png 768w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u5747\u5300\u5206\u5e03\u7684\u968f\u673a\u6570\u6d4b\u91cf\u7ed3\u679c-20230812B-548x236.png 548w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u5747\u5300\u5206\u5e03\u7684\u968f\u673a\u6570\u6d4b\u91cf\u7ed3\u679c-20230812B.png 1195w\" sizes=\"(max-width: 1024px) 100vw, 1024px\" \/><\/figure><\/div>\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E6%AD%A3%E6%80%81%E5%88%86%E5%B8%83%E7%9A%84%E9%9A%8F%E6%9C%BA%E6%95%B0\"><\/span>\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Gauss\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570\n* \u8f93\u5165\u53c2\u6570\uff1a mean\uff1a\u6b63\u6001\u5206\u5e03\u7684\u5747\u503c\u03bc\n*            sigma\uff1a\u6b63\u6001\u5206\u5e03\u7684\u5747\u65b9\u5dee\u03c3\n*            seed\uff1a\u968f\u673a\u6570\u79cd\u5b50\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a \u968f\u673a\u6570\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\ndouble Gauss(double mean, double sigma, long int* seed)\n{\n  int i;\n  double x, y;\n  for (x = 0, i = 0; i &lt; 12; i++)\n  {\n    x += Uniform(0.0, 1.0, seed);\n  }\n  x = x - 6.0;\n  y = mean + x * sigma;\n  return y;\n}\n<\/code><\/pre>\n\n\n\n<p>\u5728 -50~50 \u533a\u95f4\u5185\uff0c\u968f\u673a\u6570\u79cd\u5b50\u4e3a 13579 \u7684\u6761\u4ef6\u4e0b\uff0c\u8fde\u7eed\u83b7\u53d6\u4e00\u5343\u4e07\u4e2a\u968f\u673a\u6570\uff0c\u6d4b\u91cf\u7ed3\u679c\u5982\u4e0b\u3002<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-large\"><img decoding=\"async\" loading=\"lazy\" width=\"1024\" height=\"440\" src=\"http:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png\" alt=\"\" class=\"wp-image-1783\" srcset=\"https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png 1024w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-300x129.png 300w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-768x330.png 768w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-548x236.png 548w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570-20230812.png 1195w\" sizes=\"(max-width: 1024px) 100vw, 1024px\" \/><\/figure><\/div>\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E6%8C%87%E6%95%B0%E5%88%86%E5%B8%83%E7%9A%84%E9%9A%8F%E6%9C%BA%E6%95%B0\"><\/span>\u6307\u6570\u5206\u5e03\u7684\u968f\u673a\u6570<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Exponent\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u6307\u6570\u5206\u5e03\u7684\u968f\u673a\u6570\n* \u8f93\u5165\u53c2\u6570\uff1a beta\uff1a\u6307\u6570\u5206\u5e03\u7684\u5747\u503c\u03b2\n*            seed\uff1a\u968f\u673a\u6570\u79cd\u5b50\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a \u968f\u673a\u6570\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\ndouble Exponent(double beta, long int* seed)\n{\n  double u, x;\n  u = Uniform(0.0, 1.0, seed);\n  x = -beta * log(u);\n  return x;\n}<\/code><\/pre>\n\n\n\n<p>\u5728 -50~50 \u533a\u95f4\u5185\uff0c\u968f\u673a\u6570\u79cd\u5b50\u4e3a 13579 \u7684\u6761\u4ef6\u4e0b\uff0c\u8fde\u7eed\u83b7\u53d6\u4e00\u5343\u4e07\u4e2a\u968f\u673a\u6570\uff0c\u6d4b\u91cf\u7ed3\u679c\u5982\u4e0b\u3002<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-large\"><img decoding=\"async\" loading=\"lazy\" width=\"1024\" height=\"440\" src=\"http:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u6307\u6570\u5206\u5e03\u7684\u968f\u673a\u6570-20230812B-1024x440.png\" alt=\"\" class=\"wp-image-1791\" srcset=\"https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u6307\u6570\u5206\u5e03\u7684\u968f\u673a\u6570-20230812B-1024x440.png 1024w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u6307\u6570\u5206\u5e03\u7684\u968f\u673a\u6570-20230812B-300x129.png 300w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u6307\u6570\u5206\u5e03\u7684\u968f\u673a\u6570-20230812B-768x330.png 768w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u6307\u6570\u5206\u5e03\u7684\u968f\u673a\u6570-20230812B-548x236.png 548w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u6307\u6570\u5206\u5e03\u7684\u968f\u673a\u6570-20230812B.png 1195w\" sizes=\"(max-width: 1024px) 100vw, 1024px\" \/><\/figure><\/div>\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E6%8B%89%E6%99%AE%E6%8B%89%E6%96%AF%E5%88%86%E5%B8%83%E7%9A%84%E9%9A%8F%E6%9C%BA%E6%95%B0\"><\/span>\u62c9\u666e\u62c9\u65af\u5206\u5e03\u7684\u968f\u673a\u6570<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Laplace\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u62c9\u666e\u62c9\u65af\u5206\u5e03\u7684\u968f\u673a\u6570\n* \u8f93\u5165\u53c2\u6570\uff1a beta\uff1a\u62c9\u666e\u62c9\u65af\u5206\u5e03\u7684\u53c2\u6570\u03b2\n*            seed\uff1a\u968f\u673a\u6570\u79cd\u5b50\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a \u968f\u673a\u6570\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\ndouble Laplace(double beta, long int* seed)\n{\n  double u1, u2, x;\n  u1 = Uniform(0.0, 1.0, seed);\n  u2 = Uniform(0.0, 1.0, seed);\n  if (u1 &lt;= 0.5)\n  {\n    x = -beta * log(1.0 - u2);\n  }\n  else\n  {\n    x = beta * log(u2);\n  }\n  return x;\n}<\/code><\/pre>\n\n\n\n<p>\u5728 -50~50 \u533a\u95f4\u5185\uff0c\u968f\u673a\u6570\u79cd\u5b50\u4e3a 13579 \u7684\u6761\u4ef6\u4e0b\uff0c\u8fde\u7eed\u83b7\u53d6\u4e00\u5343\u4e07\u4e2a\u968f\u673a\u6570\uff0c\u6d4b\u91cf\u7ed3\u679c\u5982\u4e0b\u3002<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-large\"><img decoding=\"async\" loading=\"lazy\" width=\"1024\" height=\"440\" src=\"http:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u62c9\u666e\u62c9\u65af\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png\" alt=\"\" class=\"wp-image-1790\" srcset=\"https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u62c9\u666e\u62c9\u65af\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png 1024w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u62c9\u666e\u62c9\u65af\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-300x129.png 300w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u62c9\u666e\u62c9\u65af\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-768x330.png 768w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u62c9\u666e\u62c9\u65af\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-548x236.png 548w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u62c9\u666e\u62c9\u65af\u5206\u5e03\u7684\u968f\u673a\u6570-20230812.png 1195w\" sizes=\"(max-width: 1024px) 100vw, 1024px\" \/><\/figure><\/div>\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E7%91%9E%E5%88%A9%E5%88%86%E5%B8%83%E7%9A%84%E9%9A%8F%E6%9C%BA%E6%95%B0\"><\/span>\u745e\u5229\u5206\u5e03\u7684\u968f\u673a\u6570<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Rayleigh\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u745e\u5229\u5206\u5e03\u7684\u968f\u673a\u6570\n* \u8f93\u5165\u53c2\u6570\uff1a beta\uff1a\u745e\u5229\u5206\u5e03\u7684\u53c2\u6570\u03b2\n*            seed\uff1a\u968f\u673a\u6570\u79cd\u5b50\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a \u968f\u673a\u6570\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\ndouble Rayleigh(double sigma, long int* seed)\n{\n  double u, x;\n  u = Uniform(0.0, 1.0, seed);\n  x = -2.0 * log(u);\n  x = sigma * sqrt(x);\n  return x;\n}<\/code><\/pre>\n\n\n\n<p>\u5728 -50~50 \u533a\u95f4\u5185\uff0c\u968f\u673a\u6570\u79cd\u5b50\u4e3a 13579 \u7684\u6761\u4ef6\u4e0b\uff0c\u8fde\u7eed\u83b7\u53d6\u4e00\u5343\u4e07\u4e2a\u968f\u673a\u6570\uff0c\u6d4b\u91cf\u7ed3\u679c\u5982\u4e0b\u3002<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-large\"><img decoding=\"async\" loading=\"lazy\" width=\"1024\" height=\"440\" src=\"http:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u745e\u5229\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png\" alt=\"\" class=\"wp-image-1793\" srcset=\"https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u745e\u5229\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png 1024w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u745e\u5229\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-300x129.png 300w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u745e\u5229\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-768x330.png 768w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u745e\u5229\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-548x236.png 548w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u745e\u5229\u5206\u5e03\u7684\u968f\u673a\u6570-20230812.png 1195w\" sizes=\"(max-width: 1024px) 100vw, 1024px\" \/><\/figure><\/div>\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E5%AF%B9%E6%95%B0%E6%AD%A3%E6%80%81%E5%88%86%E5%B8%83%E7%9A%84%E9%9A%8F%E6%9C%BA%E6%95%B0\"><\/span>\u5bf9\u6570\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Lognorm\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u5bf9\u6570\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570\n* \u8f93\u5165\u53c2\u6570\uff1a mean\uff1a\u5bf9\u6570\u6b63\u6001\u5206\u5e03\u7684\u53c2\u6570\u03bc\n*            sigma\uff1a\u5bf9\u6570\u6b63\u6001\u5206\u5e03\u7684\u53c2\u6570\u03c3\n*            seed\uff1a\u968f\u673a\u6570\u79cd\u5b50\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a \u968f\u673a\u6570\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\ndouble Lognorm(double mean, double sigma, long int* seed)\n{\n  double x, y;\n  y = Gauss(mean, sigma, seed);\n  x = exp(y);\n  return x;\n}<\/code><\/pre>\n\n\n\n<p>\u5728 0~10 \u533a\u95f4\u5185\uff0c\u968f\u673a\u6570\u79cd\u5b50\u4e3a 13579 \u7684\u6761\u4ef6\u4e0b\uff0c\u8fde\u7eed\u83b7\u53d6\u4e00\u5343\u4e07\u4e2a\u968f\u673a\u6570\uff0c\u6d4b\u91cf\u7ed3\u679c\u5982\u4e0b\u3002<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-large\"><img decoding=\"async\" loading=\"lazy\" width=\"1024\" height=\"440\" src=\"http:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u5bf9\u6570\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png\" alt=\"\" class=\"wp-image-1796\" srcset=\"https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u5bf9\u6570\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png 1024w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u5bf9\u6570\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-300x129.png 300w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u5bf9\u6570\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-768x330.png 768w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u5bf9\u6570\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-548x236.png 548w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u5bf9\u6570\u6b63\u6001\u5206\u5e03\u7684\u968f\u673a\u6570-20230812.png 1195w\" sizes=\"(max-width: 1024px) 100vw, 1024px\" \/><\/figure><\/div>\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E6%9F%AF%E8%A5%BF%E5%88%86%E5%B8%83%E7%9A%84%E9%9A%8F%E6%9C%BA%E6%95%B0\"><\/span>\u67ef\u897f\u5206\u5e03\u7684\u968f\u673a\u6570<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Cauchy\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u67ef\u897f\u5206\u5e03\u7684\u968f\u673a\u6570\n* \u8f93\u5165\u53c2\u6570\uff1a alpha\uff1a\u67ef\u897f\u5206\u5e03\u7684\u53c2\u6570\u03b1\n*            beta\uff1a\u67ef\u897f\u5206\u5e03\u7684\u53c2\u6570\u03b2\n*            seed\uff1a\u968f\u673a\u6570\u79cd\u5b50\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a \u968f\u673a\u6570\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\ndouble Cauchy(double alpha, double beta, long int* seed)\n{\n  double u, x;\n  u = Uniform(0.0, 1.0, seed);\n  x = alpha - beta \/ tan(3.1415926 * u);\n  return x;\n}<\/code><\/pre>\n\n\n\n<p>\u5728 -50~50 \u533a\u95f4\u5185\uff0c\u968f\u673a\u6570\u79cd\u5b50\u4e3a 13579 \u7684\u6761\u4ef6\u4e0b\uff0c\u8fde\u7eed\u83b7\u53d6\u4e00\u5343\u4e07\u4e2a\u968f\u673a\u6570\uff0c\u6d4b\u91cf\u7ed3\u679c\u5982\u4e0b\u3002<\/p>\n\n\n\n<figure class=\"wp-block-image size-large\"><img decoding=\"async\" loading=\"lazy\" width=\"1024\" height=\"440\" src=\"http:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u67ef\u897f\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png\" alt=\"\" class=\"wp-image-1798\" srcset=\"https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u67ef\u897f\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png 1024w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u67ef\u897f\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-300x129.png 300w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u67ef\u897f\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-768x330.png 768w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u67ef\u897f\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-548x236.png 548w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u67ef\u897f\u5206\u5e03\u7684\u968f\u673a\u6570-20230812.png 1195w\" sizes=\"(max-width: 1024px) 100vw, 1024px\" \/><\/figure>\n\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E9%9F%A6%E4%BC%AF%E5%88%86%E5%B8%83%E7%9A%84%E9%9A%8F%E6%9C%BA%E6%95%B0\"><\/span>\u97e6\u4f2f\u5206\u5e03\u7684\u968f\u673a\u6570<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Weibull\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u97e6\u4f2f\u5206\u5e03\u7684\u968f\u673a\u6570\n* \u8f93\u5165\u53c2\u6570\uff1a alpha\uff1a\u97e6\u4f2f\u5206\u5e03\u7684\u53c2\u6570\u03b1\n*            beta\uff1a\u97e6\u4f2f\u5206\u5e03\u7684\u53c2\u6570\u03b2\n*            seed\uff1a\u968f\u673a\u6570\u79cd\u5b50\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a \u968f\u673a\u6570\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\ndouble Weibull(double alpha, double beta, long int* seed)\n{\n  double u, x;\n  u = Uniform(0.0, 1.0, seed);\n  u = -log(u);\n  x = beta * pow(u, 1.0 \/ alpha);\n  return x;\n}<\/code><\/pre>\n\n\n\n<p>\u5728 0~5 \u533a\u95f4\u5185\uff0c\u968f\u673a\u6570\u79cd\u5b50\u4e3a 13579 \u7684\u6761\u4ef6\u4e0b\uff0c\u8fde\u7eed\u83b7\u53d6\u4e00\u5343\u4e07\u4e2a\u968f\u673a\u6570\uff0c\u6d4b\u91cf\u7ed3\u679c\u5982\u4e0b\u3002<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-large\"><img decoding=\"async\" loading=\"lazy\" width=\"1024\" height=\"440\" src=\"http:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u97e6\u4f2f\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png\" alt=\"\" class=\"wp-image-1799\" srcset=\"https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u97e6\u4f2f\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png 1024w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u97e6\u4f2f\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-300x129.png 300w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u97e6\u4f2f\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-768x330.png 768w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u97e6\u4f2f\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-548x236.png 548w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u97e6\u4f2f\u5206\u5e03\u7684\u968f\u673a\u6570-20230812.png 1195w\" sizes=\"(max-width: 1024px) 100vw, 1024px\" \/><\/figure><\/div>\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E5%9F%83%E5%B0%94%E6%9C%97%E5%88%86%E5%B8%83%E7%9A%84%E9%9A%8F%E6%9C%BA%E6%95%B0\"><\/span>\u57c3\u5c14\u6717\u5206\u5e03\u7684\u968f\u673a\u6570<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Erlang\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u57c3\u5c14\u6717\u5206\u5e03\u7684\u968f\u673a\u6570\n* \u8f93\u5165\u53c2\u6570\uff1a m\uff1a\u57c3\u5c14\u6717\u5206\u5e03\u7684\u53c2\u6570m\n*            beta\uff1a\u57c3\u5c14\u6717\u5206\u5e03\u7684\u53c2\u6570\u03b2\n*            seed\uff1a\u968f\u673a\u6570\u79cd\u5b50\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a \u968f\u673a\u6570\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\ndouble Erlang(double m, double beta, long int* seed)\n{\n  int i;\n  double u, x;\n  for (u = 1.0, i = 0; i &lt; m; i++)\n  {\n    u *= Uniform(0.0, 1.0, seed);\n  }\n  x = -beta * log(u);\n  return x;\n}<\/code><\/pre>\n\n\n\n<p>\u5728 0~20 \u533a\u95f4\u5185\uff0c\u968f\u673a\u6570\u79cd\u5b50\u4e3a 13579 \u7684\u6761\u4ef6\u4e0b\uff0c\u8fde\u7eed\u83b7\u53d6\u4e00\u5343\u4e07\u4e2a\u968f\u673a\u6570\uff0c\u6d4b\u91cf\u7ed3\u679c\u5982\u4e0b\u3002<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-large\"><img decoding=\"async\" loading=\"lazy\" width=\"1024\" height=\"440\" src=\"http:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u57c3\u5c14\u6717\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png\" alt=\"\" class=\"wp-image-1802\" srcset=\"https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u57c3\u5c14\u6717\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-1024x440.png 1024w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u57c3\u5c14\u6717\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-300x129.png 300w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u57c3\u5c14\u6717\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-768x330.png 768w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u57c3\u5c14\u6717\u5206\u5e03\u7684\u968f\u673a\u6570-20230812-548x236.png 548w, https:\/\/www.huangrongzhen.ink\/wp-content\/uploads\/2023\/08\/\u57c3\u5c14\u6717\u5206\u5e03\u7684\u968f\u673a\u6570-20230812.png 1195w\" sizes=\"(max-width: 1024px) 100vw, 1024px\" \/><\/figure><\/div>\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E8%B4%9D%E5%8A%AA%E9%87%8C%E5%88%86%E5%B8%83%E7%9A%84%E9%9A%8F%E6%9C%BA%E6%95%B0\"><\/span>\u8d1d\u52aa\u91cc\u5206\u5e03\u7684\u968f\u673a\u6570<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Bn\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u8d1d\u52aa\u91cc\u5206\u5e03\u7684\u968f\u673a\u6570\n* \u8f93\u5165\u53c2\u6570\uff1a m\uff1a\u8d1d\u52aa\u91cc\u5206\u5e03\u7684\u53c2\u6570p\n*            seed\uff1a\u968f\u673a\u6570\u79cd\u5b50\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a \u968f\u673a\u6570\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\nint Bn(double p, long int* seed)\n{\n  int x;\n  double u;\n  u = Uniform(0.0, 1.0, seed);\n  x = (u &lt;= p) ? 1:0;\n  return x;\n}<\/code><\/pre>\n\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E8%B4%9D%E5%8A%AA%E9%87%8C-%E9%AB%98%E6%96%AF%E5%88%86%E5%B8%83%E7%9A%84%E9%9A%8F%E6%9C%BA%E6%95%B0\"><\/span>\u8d1d\u52aa\u91cc-\u9ad8\u65af\u5206\u5e03\u7684\u968f\u673a\u6570<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Bg\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u8d1d\u52aa\u91cc-\u9ad8\u65af\u5206\u5e03\u7684\u968f\u673a\u6570\n* \u8f93\u5165\u53c2\u6570\uff1a m\uff1a\u8d1d\u52aa\u91cc\u5206\u5e03\u7684\u53c2\u6570p\n*            mean\uff1a\u9ad8\u65af\u5206\u5e03\u7684\u5747\u503c\u03bc\n*            sigma\uff1a\u9ad8\u65af\u5206\u5e03\u7684\u5747\u65b9\u5dee\u03c3\n*            seed\uff1a\u968f\u673a\u6570\u79cd\u5b50\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a \u968f\u673a\u6570\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\ndouble Bg(double p, double mean, double sigma, long int* seed)\n{\n  double u, x;\n  u = Uniform(0.0, 1.0, seed);\n  if (u &lt;= p)\n  {\n    x = Gauss(mean, sigma, seed);\n  }\n  else\n  {\n    x = 0.0;\n  }\n  return x;\n}\n<\/code><\/pre>\n\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E4%BA%8C%E9%A1%B9%E5%BC%8F%E5%88%86%E5%B8%83%E7%9A%84%E9%9A%8F%E6%9C%BA%E6%95%B0\"><\/span>\u4e8c\u9879\u5f0f\u5206\u5e03\u7684\u968f\u673a\u6570<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Bin\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u4e8c\u9879\u5f0f\u5206\u5e03\u7684\u968f\u673a\u6570\n* \u8f93\u5165\u53c2\u6570\uff1a n\uff1a\u4e8c\u9879\u5f0f\u5206\u5e03\u7684\u53c2\u6570n\n*            p\uff1a\u4e8c\u9879\u5f0f\u5206\u5e03\u7684\u53c2\u6570\u03bc\n*            seed\uff1a\u968f\u673a\u6570\u79cd\u5b50\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a \u968f\u673a\u6570\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\nint Bin(int n, double p, long int* seed)\n{\n  int i;\n  double x;\n  for (x = 0.0, i = 0; i &lt; n; i++)\n  {\n    x += Bn(p, seed);\n  }\n  return x;\n}<\/code><\/pre>\n\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E6%B3%8A%E6%9D%BE%E5%88%86%E5%B8%83%E7%9A%84%E9%9A%8F%E6%9C%BA%E6%95%B0\"><\/span>\u6cca\u677e\u5206\u5e03\u7684\u968f\u673a\u6570<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Poisson\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u6cca\u677e\u5206\u5e03\u7684\u968f\u673a\u6570\n* \u8f93\u5165\u53c2\u6570\uff1a n\uff1a\u6cca\u677e\u5206\u5e03\u7684\u5747\u503c\u03bb\n*            seed\uff1a\u968f\u673a\u6570\u79cd\u5b50\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a \u968f\u673a\u6570\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\nint Poisson(double lambda, long int* seed)\n{\n  int i, x;\n  double a, b, u;\n  a = exp(-lambda);\n  i = 0;\n  b = 1.0;\n  do\n  {\n    u = Uniform(0.0, 1.0, seed);\n    b *= u;\n    i++;\n  }while(b &gt;= a);\n  x = i - 1;\n  return x;\n}<\/code><\/pre>\n\n\n\n<p>\u5728 -50~50 \u533a\u95f4\u5185\uff0c\u968f\u673a\u6570\u79cd\u5b50\u4e3a 13579 \u7684\u6761\u4ef6\u4e0b\uff0c\u8fde\u7eed\u83b7\u53d6\u4e00\u5343\u4e07\u4e2a\u968f\u673a\u6570\uff0c\u6d4b\u91cf\u7ed3\u679c\u5982\u4e0b\u3002<\/p>\n\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"ARMApq_%E6%A8%A1%E5%9E%8B%E6%95%B0%E6%8D%AE\"><\/span>ARMA(p,q) \u6a21\u578b\u6570\u636e<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a ARMA\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210 ARMA(p,q) \u6a21\u578b\u6570\u636e\n* \u8f93\u5165\u53c2\u6570\uff1a a\uff1a\u53cc\u7cbe\u5ea6\u5b9e\u578b\u4e00\u7ef4\u6570\u7ec4\uff0c\u957f\u5ea6\u4e3a\uff08p+1\uff09\u3002ARMA(p,q) \u6a21\u578b\u7684\u81ea\u56de\u5f52\u7cfb\u6570\u3002\n*            b\uff1a\u53cc\u7cbe\u5ea6\u5b9e\u578b\u4e00\u7ef4\u6570\u7ec4\uff0c\u957f\u5ea6\u4e3a\uff08p+1\uff09\u3002ARMA(p,q) \u6a21\u578b\u7684\u6ed1\u52a8\u5e73\u5747\u7cfb\u6570\u3002\n*            p\uff1a\u6574\u5f62\u53d8\u91cf\u3002ARMA(p,q) \u6a21\u578b\u7684\u81ea\u56de\u5f52\u9636\u6570\u3002\n*            q\uff1a\u6574\u5f62\u53d8\u91cf\u3002ARMA(p,q) \u6a21\u578b\u7684\u6ed1\u52a8\u5e73\u5747\u9636\u6570\u3002\n*            mean\uff1a\u53cc\u7cbe\u5ea6\u5b9e\u578b\u53d8\u91cf\u3002\u4ea7\u751f\u767d\u566a\u58f0\u6240\u7528\u7684\u6b63\u6001\u5206\u5e03\u7684\u5747\u503c\u03bc\u3002\n*            sigma\uff1a\u53cc\u7cbe\u5ea6\u5b9e\u578b\u53d8\u91cf\u3002\u4ea7\u751f\u767d\u566a\u58f0\u6240\u7528\u7684\u6b63\u6001\u5206\u5e03\u7684\u5747\u65b9\u5dee\u03c3\n*            seed\uff1a\u957f\u6574\u578b\u6307\u9488\u53d8\u91cf\u3002*seed \u4e3a\u968f\u673a\u6570\u79cd\u5b50\u3002\n*            x\uff1a\u53cc\u7cbe\u5ea6\u5b9e\u578b\u4e00\u7ef4\u6570\u7ec4\uff0c\u957f\u5ea6\u4e3a n\u3002\u5b58\u653e ARMA(p,q) \u6a21\u578b\u7684\u6570\u636e\u3002\n*            n\uff1a\u6574\u578b\u53d8\u91cf\u3002ARMA(p,q) \u6a21\u578b\u6570\u636e\u7684\u957f\u5ea6\u3002\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a void\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670810\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\nvoid ARMA(double* a, double* b, int p, int q, double mean, double sigma, long int* seed, double* x, int n)\n{\n  int i, k, m;\n  double s;\n  double* w;\n  w = malloc(n * sizeof(double));\n  for (k = 0; k &lt; n; k++)\n  {\n    w[k] = Gauss(mean, sigma, seed);\n  }\n  x[0] = b[0] * w[0];\n  for (k = 1; k &lt;= p; k++)\n  {\n    s = 0.0;\n    for (i = 1; i &lt;= k; i++)\n    {\n      s += a[i] * x[k - i];\n    }\n    s = b[0] * w[k] - s;\n    if (q == 0)\n    {\n      x[k] = s;\n      continue;\n    }\n    m = (k &gt; q)? q:k;\n    for (i = 1; i &lt;= m; i++)\n    {\n      s += b[i]*w[k - i];\n    }\n    x[k] = s;\n  }\n  for (k = (p + 1); k &lt;= n; k++)\n  {\n    s = 0.0;\n    for (i = 1; i &lt;= p; i++)\n    {\n      s += a[i] * x[k - i];\n    }\n    s = b[0] * w[k] - s;\n    if (q == 0)\n    {\n      x[k] = s;\n      continue;\n    }\n    for (i = 1; i &lt;= q; i++)\n    {\n      s += b[i] * w[k - i];\n    }\n    x[k] = s;\n  }\n  free(w);\n}<\/code><\/pre>\n\n\n\n<h2 class=\"wp-block-heading\"><span class=\"ez-toc-section\" id=\"%E5%90%AB%E6%9C%89%E9%AB%98%E6%96%AF%E7%99%BD%E5%99%AA%E5%A3%B0%E7%9A%84%E6%AD%A3%E5%BC%A6%E7%BB%84%E5%90%88%E4%BF%A1%E5%8F%B7\"><\/span>\u542b\u6709\u9ad8\u65af\u767d\u566a\u58f0\u7684\u6b63\u5f26\u7ec4\u5408\u4fe1\u53f7<span class=\"ez-toc-section-end\"><\/span><\/h2>\n\n\n\n<p>\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code lang=\"c\" class=\"language-c\">\/*********************************************************************************************************\n* \u51fd\u6570\u540d\u79f0\uff1a Sinwn\n* \u51fd\u6570\u529f\u80fd\uff1a \u751f\u6210\u542b\u6709\u9ad8\u65af\u767d\u566a\u58f0\u7684\u6b63\u5f26\u7ec4\u5408\u4fe1\u53f7\n* \u8f93\u5165\u53c2\u6570\uff1a a\uff1a\u53cc\u7cbe\u5ea6\u5b9e\u578b\u4e00\u7ef4\u6570\u7ec4\uff0c\u957f\u5ea6\u4e3a m\u3002\u5404\u6b63\u5f26\u4fe1\u53f7\u7684\u632f\u5e45\u3002\n*            f\uff1a\u53cc\u7cbe\u5ea6\u5b9e\u578b\u4e00\u7ef4\u6570\u7ec4\uff0c\u957f\u5ea6\u4e3a m\u3002\u5404\u6b63\u5f26\u4fe1\u53f7\u7684\u9891\u7387\u3002\n*            ph\uff1a\u53cc\u7cbe\u5ea6\u5b9e\u578b\u4e00\u7ef4\u6570\u7ec4\uff0c\u957f\u5ea6\u4e3a m\u3002\u5404\u6b63\u5f26\u4fe1\u53f7\u7684\u76f8\u4f4d\u3002\n*            m\uff1a\u6574\u578b\u53d8\u91cf\u3002\u6b63\u5f26\u4fe1\u53f7\u7684\u4e2a\u6570\u3002\n*            fs\uff1a\u53cc\u7cbe\u5ea6\u5b9e\u578b\u53d8\u91cf\u3002\u91c7\u6837\u9891\u7387\uff0c\u5355\u4f4d\u662f Hz\u3002\n*            snr\uff1a\u53cc\u7cbe\u5ea6\u5b9e\u578b\u53d8\u91cf\u3002\u4fe1\u566a\u6bd4\uff0c\u5355\u4f4d\u662f dB\u3002\n*            seed\uff1a\u957f\u6574\u578b\u53d8\u91cf\u3002\u968f\u673a\u6570\u7684\u79cd\u5b50\u3002\n*            x\uff1a\u53cc\u7cbe\u5ea6\u5b9e\u578b\u4e00\u7ef4\u6570\u7ec4\uff0c\u957f\u5ea6\u4e3a n\u3002\u5b58\u653e\u6240\u4ea7\u751f\u7684\u6570\u636e\n*            n\uff1a\u6574\u578b\u53d8\u91cf\u3002\u6570\u636e\u957f\u5ea6\n* \u8f93\u51fa\u53c2\u6570\uff1a void\n* \u8fd4 \u56de \u503c\uff1a void\n* \u521b\u5efa\u65e5\u671f\uff1a 2023\u5e7408\u670811\u65e5\n* \u6ce8    \u610f\uff1a\n*********************************************************************************************************\/\nvoid Sinwn(double* a, double* f, double* ph, int m, double fs, double snr, long int* seed, double* x, int n)\n{\n  int i, k;\n  double z, pi, nsr;\n  pi = 4.0 * atan(1.0);\n  z = snr \/ 10.0;\n  z = pow(10.0, z);\n  z = 1.0 \/ (2.0 * z);\n  nsr = sqrt(z);\n  for (i = 0; i &lt; m; i++)\n  {\n    f[i] = 2.0 * pi * f[i] \/ fs;\n    ph[i] = ph[i] * pi \/ 180.0;\n  }\n  for (k = 0; k &lt; n; k++)\n  {\n    x[k] = 0.0;\n    for (i = 0; i &lt; m; i++)\n    {\n      x[k] = x[k] + a[i] * sin(k * f[i] + ph[i]);\n    }\n    x[k] = x[k] + nsr * Gauss(0.0, 1.0, seed);\n  }\n}<\/code><\/pre>\n","protected":false},"excerpt":{"rendered":"<p>\u7b97\u6cd5\u7b14\u8bb0<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[6],"tags":[],"_links":{"self":[{"href":"https:\/\/www.huangrongzhen.ink\/index.php?rest_route=\/wp\/v2\/posts\/1749"}],"collection":[{"href":"https:\/\/www.huangrongzhen.ink\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.huangrongzhen.ink\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.huangrongzhen.ink\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.huangrongzhen.ink\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=1749"}],"version-history":[{"count":20,"href":"https:\/\/www.huangrongzhen.ink\/index.php?rest_route=\/wp\/v2\/posts\/1749\/revisions"}],"predecessor-version":[{"id":1804,"href":"https:\/\/www.huangrongzhen.ink\/index.php?rest_route=\/wp\/v2\/posts\/1749\/revisions\/1804"}],"wp:attachment":[{"href":"https:\/\/www.huangrongzhen.ink\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=1749"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.huangrongzhen.ink\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=1749"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.huangrongzhen.ink\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=1749"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}