pasuka 发表于 2014-10-14 09:14
国防项目的话，通过单位或组织去联系中国工程物理研究院、海军装备研究院、中国船舶科技中心等部门
汽车 ...

非常感谢.

hanshikai 发表于 2014-10-14 08:55
看来是内行啊!
我的问题和编程真的没有半毛钱关系?请指示

国防项目的话，通过单位或组织去联系中国工程物理研究院、海军装备研究院、中国船舶科技中心等部门
汽车、建筑的不了解

pasuka 发表于 2014-10-14 08:44
1、冲击谱相关标准本来就是先有实验数据再有拟合经验公式，若是船舶的冲击谱，国外不知道炸多少条船做实验 ...

看来是内行啊!
我的问题和编程真的没有半毛钱关系?请指示

1、冲击谱相关标准本来就是先有实验数据再有拟合经验公式，若是船舶的冲击谱，国外不知道炸多少条船做实验，汽车行业就更不用提了；
2、lz的问题和编程看不出有任何联系，理论问题还是请教原文作者或者翻阅标准和相关文献

aliouying 发表于 2014-10-13 10:20
我对你的理论不是很清楚，是否需要添加误差，然后再进行修正？

看了下您的程序，貌似只是对理论的数据进行 ...

感谢您的关注!
1. 误差已经在module NORM_中给出即实常量Tol=3;
2. 本程序所要做的就是对幅值A进行修正,A(:)这个数组中各个元素是不发生联系的;
3. 时域波形是通过一系列特定频率的基波合成的,或者说把时域波形离散成一些列基波的和,在每个时间点上都要做这样的合成.程序中WavMat这个数组存储了所有的时域信息,他的行代表时间,列代表基波序列,按照列来代入有关计算公式(即subroutine Recs_Filt,一次仅代入一列)计算,得到是时间上的结果,挑选最大值就是该基频对应的响应,依次代入WavMat各列,得到各频率对应的响应.附图即是频率和响应关系图(横纵坐标均取常用对数);
4. 上面得到的响应不可避免的会带有误差,SelectCase(nUpDat)中的Case(2)就是修正数据的过程,在检查误差之前有一次强行修正.

再次感谢,不知道还有什么地方没讲清楚.

我对你的理论不是很清楚，是否需要添加误差，然后再进行修正？

看了下您的程序，貌似只是对理论的数据进行修正。