山野莽夫

  • 归档
    • 随笔
    • 建站资源
    • 分享
    • 代码
  • 地球物理学
    • 专业课
    • 概念解释
  • 计算机
  • 互联网
  • 教程
  • 规划
  • 实验室
    • 珍藏的软件
    • 贴吧云签到
    • A1账号自助申请
山野莽夫
小学生的挣扎的点点滴滴
  1. 首页
  2. 代码
  3. c
  4. 正文

地震记录合成代码

2016年4月27日 2324点热度 0人点赞 0条评论

 

工程物探上有一个实验要求合成地震记录,还是比较简单的,就是反射系数序列和雷克子波卷积就可。卷积公式这儿就不写出来了。反射系数序列可以自己设置,200个点同时赋值为0,地下有几层,就把几个点设置为其他的值。对于雷克子波,有公式就可以做出来,对于工程物探,频率尽量大一些,这样地震记录的图要好看一些。

QQ截图20160427000839

下边把代码奉献给大家。

 

#include<stdio.h>
#include<math.h>
#include<stdlib.h>
//#define dz1 "file1.txt"
//#define dz2 "file2.txt"

main ()
{ FILE *fp1,*fp2;
double dt=0.002,t[71];int fm1=70,i,j;double zibo[71],xishu[200],dzjl[300];
// fp1=fopen(dz1,"w+");
// fp2=fopen(dz2,"w+");
// fclose(fp1);
//fclose(fp2);

fp1=fopen("file1.txt","r+");
fp2=fopen("file2.txt","r+");
for(i=0;i<=35;i++)

{
t[i+35]=i*dt;

zibo[i+35]=exp(-pow(3.141592653*fm1*t[i+35],2))*(1-2*pow(3.141592653*fm1*t[i+35],2));

}
//for(i=0;i<71;i++)
// printf("%f\n",zibo[i]);
for(i=0;i<=34;i++)
{
t[i]=-t[70-i];
zibo[i]=zibo[70-i];

}
for(i=0;i<=70;i++)
{printf("%f\t%f\n",t[i],zibo[i]);
fprintf(fp1," %f\t%f\n",t[i],zibo[i]);

}

for(i=0;i<=199;i++)
{
xishu[i]=0;
}
xishu[10]=0.2;xishu[65]=0.4;xishu[110]=-0.5;xishu[155]=0.7;xishu[190]=-0.36;
for(i=0;i<200;i++)

{ dzjl[i]=0.0;

for(j=0;j<70;j++)

{ if((i-j)>0)

dzjl[i]=dzjl[i]+zibo[j]*xishu[i-j] ; }

}
for(i=0;i<200;i++)
{
// printf("%d\t%f\n",i,dzjl[i]);
fprintf(fp2,"%d\t%f\n",i,dzjl[i]);
}
fclose(fp1);
fclose(fp2);

}

标签: 卷积 地震记录 雷克子波
最后更新:2016年6月22日

小菜菜

菜鸟

打赏 点赞
< 上一篇
下一篇 >

文章评论

razz evil exclaim smile redface biggrin eek confused idea lol mad twisted rolleyes wink cool arrow neutral cry mrgreen drooling persevering
取消回复

此站点使用Akismet来减少垃圾评论。了解我们如何处理您的评论数据。

标签聚合
c语言 模板 onedrive 虚拟机 wordpress 宝塔面板 ppt 地震学程序
最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
数据结构(C语言版)第2版参考答案 写一写本人买斐讯k3的经历 八位二进制数为什么表示范围是:-128~~+127? 基于Gassmann方程的流体替换技术研究 Winrar 压缩软件 去弹窗广告注册纯净版含注册key文件 宝塔面板配置onedrive索引工具OLAINDEX

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2