二维弹性介质波动方程数值模拟动态图像

yeahamen

已于 2023-09-10 21:27:07 修改

阅读量468

点赞数 1

文章标签：算法开发语言 C 波动方程数值模拟

于 2022-01-24 12:00:38 首次发布

本文链接：https://blog.csdn.net/amen_yanwang/article/details/122664689

版权

/*Forward modeling*/
#include<stdio.h>
#include<math.h>
#include<graphics.h>
#include<conio.h>
#include<windows.h>
#include<stdlib.h>
void fun();
int main (void)
{
   	
	
	/*Background board*/
	initgraph(600,600);       
	setbkcolor(BLACK);         
	setlinecolor(WHITE);       
	cleardevice();             
	setorigin(2,598);

	fun();
	return 0;
}
void fun()
{
   
	/*Initializing Variables*/
	int x,y,t;
	double i;
	double j;
	double dx,dy,dt,dh,v;
	double b[179][179];
	double constant;
	double ricker = 100;
	double m,n;
	int m1,n1;
	double (*p)[179][179]=(double(*)[179][179])malloc(sizeof(double)*175*179*179);
	v = 0.5;
	dt = dx = dy = dh = 1;	
	constant = (v*dt/dh)*(v*dt/dh);

	/*Gird depiction*/
	for(j=2.0;j<=587.4;j+=3.3)
	{
   
		line(2.0,-j,586.1,-j);   
		line(j, -586.1,j,-2.0);       
	}


	/*Boundary initialization*/
	for(t=0;t<175;t++)
	{
   
		for(x=0;x<179;x++)
		{
   
			for(y=0;y<179;y++)
			{
   
				p[t][x][y] = 0;//the value of 'p