自定义波源-Gaussian脉冲.doc

《自定义波源-Gaussian脉冲.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、高斯函数：clearall;closeall;x=-4:0.1:4;y=(2*pi)^(-1/2).*exp(-x.^2/2);yy=diff(y);plot(y);f=gaussmf(x,[25]);ff=diff(f);plot(f);假设波源为f(x)=gaussmf(x,[sigc]),f’(x)=cos(x),其中f’(x)为f(x)的导数，为了能非常详细的说明，下面我从如何算出数据开始一步一步的来1、编写matlab代码(计算波源函数及其导数在离散点的值)如下：clearall;closeall;x=0:0.1:10

2、;y=gaussmf(x,[25]);fid=fopen('y.txt','w');fprintf(fid,'%g',y);fclose(fid);fx=diff(gaussmf(x,[25]));fid=fopen('fx.txt','w');fprintf(fid,'%g',fx);fclose(fid);对于其他的函数如高斯函数，只用换掉其中的波源函数及其导数即可。运行上面的函数，在你的函数保存文件夹中将会生成y.txt和yy.txt两个文件。其中y.txt为波源函数在各个离散点的值，yy.txt为波源函数导数的值

3、。现在回到你的matlab的workspace中输入a=[y’yy’]如下图就可以把波源函数在各个离散点的值和波源函数导数的值合并成两列的矩阵，运用size（a）可以获得矩阵的大小，即离散个数如下图表明有1885个离散值现在就开始来制作我们需要的userdefined文件有了上面几步在matlab的workspace左上角有a矩阵的值，图下图所示双击其中的a得如下图其中第一列为波源函数值，第二列为波源函数导数值点击第一列和第二列，然后右键copy把上面copy的数据paste到一个你建立好的**.txt文件中，本示例为文件打开为

4、然后根据xfdtd自定义波源，要在上面文件的第一行加上离散个数大小加上后为下图其中1885由上面的size(a)而得至此，自定义波源文件建立完成，不是很麻烦吧，哈哈下面导入到xfdtd中看看其显示是否正确上面就是导入文件后在xfdtd软件中的显示，完全正确，ok！搞定了