当教室突然安静·二——衰减与“静默波”

#include 
using namespace std;
int length=10;//人阵规模，不得超出10000*10000
int width=10;
//下为f二三项
double A(double a,double b)
{
	double y=log1p(a)/log(2)*0.5*(1+tanh(5*(b+1.5)));
	return y;
}
//下为随机函数
double g(double t)
{
	double y=1+tanh(7*(sin((sqrt(3))*0.7*t)+sin((sqrt(5))*0.7*t)+0.332*sin((sqrt(16))*0.7*t)+sin((sqrt(14))*0.7*t)+1.02*sin((sqrt(2.5803)*0.7*t))));
	return y;
}
//这个g2是用来错开各个f的随机初始值，为了省事直接把g里面的部分拿出来了
double g2(double t)
{
	double y=(sin((sqrt(3))*0.7*t)+sin((sqrt(5))*0.7*t)+0.332*sin((sqrt(16))*0.7*t)+sin((sqrt(14))*0.7*t)+1.02*sin((sqrt(2.5803)*0.7*t)));
	return y;
}
//递推步进函数
double Nextf(int step,vector ramvec,double deltat)
{
	double t=0.05*step;
	double fl=(ramvec[(step+19)%20]+ramvec[(step+18)%20]+ramvec[(step+17)%20]+ramvec[(step+16)%20]+ramvec[(step+15)%20]+ramvec[(step+14)%20])/6;
	double fv=(((ramvec[(step+19)%20]+ramvec[(step+18)%20]+ramvec[(step+17)%20]+ramvec[(step+16)%20]+ramvec[(step+15)%20])/5)-((ramvec[(step+9)%20]+ramvec[(step+8)%20]+ramvec[(step+7)%20]+ramvec[(step+6)%20]+ramvec[(step+5)%20])/5))*2;
	double ft=g(t-deltat)*A(fl,fv)+0.001*g(t-deltat);
	return ft;
}

double fdata[10000][10000];//传出响度定义
//卷积核(迫真)就是模型里的S(……),决定声音如何传播，衰减
double core(int i,int j)//接收坐标
{
	double fsum=0;
	double f1sum=0;
	double daor=0;
	for(int m=0;m fsave[width][length];//历史接受响度初始化
	vector example(20,1);
	for(int i=0;i tdata={};//突变时间记录
//	int count=0;//突变次数
//	int test=0;
//	int last=0;
	long long int step;//这里step指第几步而非步长，步长固定在0.05
	//———————————————————————————————————递推循环————————————————醒目的分界线————————————————————————————■■■■■■■■■■———— 
	for(step=0;step<=(20*1000);step=step+1)//迭代循环,1000s用15s
	{
		double t=0.05*step;
		//这里放迭代主体
		//这里放卷积
		for(int i=0;i=0.1)
//		{
//			test=0;
//		}
//		else
//		{
//			test=1;
//		}
//		
//		if(test-last==1)
//		{
//			count=count+1;
//			tdata.push_back(t);
//		}
//		last=test;
		
		for (int i = 0; i < width; ++i) {//把每个step的矩阵拍扁后输出，这里代码排版莫名其妙炸了
    	    for (int j = 0; j < length; ++j) {
   	         cout<

import os
import math
import numpy as np
import matplotlib.pyplot as plt
current_dir = os.getcwd()

folder_path = "/storage/emulated/0/Documents/Pydroid3"
file_name = "long.txt"
 
full_path = os.path.join(current_dir, folder_path, file_name)

with open(full_path, 'r') as file:
	content = file.read().splitlines() 
	list = [[float(x) for x in line.split()]for line in content] 
fig=plt.figure(figsize=(10,10))
fig,ax = plt.subplots()
for step in range(0,1200):#把被拍扁的矩阵立起来
	for i in range(0,1000):
		a=math.tanh(list[step][i])
		rect=plt.Rectangle((i%200,(i//200)),1,1,alpha=a)
		ax.add_patch(rect)
		plt.axis([0, 200, 0, 5])
	
	plt.savefig('/storage/emulated/0/Documents/Pydroid3/longaa/'+str(step)+'.png')
	ax.cla()