matlab 离散点求导_【Spectral Method】8.谱求导矩阵求解波动方程

最新推荐文章于 2021-10-18 15:58:57 发布
最新推荐文章于 2021-10-18 15:58:57 发布
阅读量670 收藏
点赞数
文章标签： matlab 离散点求导
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_34010875/article/details/112065227
版权

  1. 1D情况

 想象这么一个情景：有一根两端固定的、柔韧度足够的弦，若它的平衡位置设置为x轴，u(x,t)表示弦上坐标为x的点在时刻t垂直于x方向的位移，则弦的自由震动方程为：

02deb559766537bf8fd08e01c60ccee5.png(可以取a=1)

   依旧拿之前的例子作为初始条件，它对于波动的表现较形象：

4fada1b6ba7349dfccfabcd9dce230c3.png

0a8671ce6a94ae03b9efceb00b0d311e.png

  《数理方程》中介绍过达朗贝尔公式，根据公式，这2个sech函数将分别分裂成两半，以速度a沿着x轴的正方向和负方向传播，即正行波和反行波，这些正行波和反行波的叠加就形成了弦的总位移。

基本思路就是：

最低0.47元/天 解锁文章
不爱提问的张学友
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
matlab对信号求导,符号函数的微分（求导）在Matlab中实现
weixin_28871989的博客
03-16 1646
MATLAB软件提供求函数导数的指令是diff，具体使用格式如下：(1)diff(f, x)表示对f(这里f是一个函数表达式)求关于符号变量x的一阶导数．若x缺省，则表示求f对预设独立变量的一阶导数．(2)diff(f, x, n) 表示对f求关于符号变量x的n阶导数．若x缺省，则表示求f对预设独立变量的n阶导数．例1 ax2+bx+c的一阶导二阶导解>>syms a b c x&...
离散数据求导
01-26
离散求导离散求导离散求导离散求导离散求导离散求导离散求导
第一话 统计计算之随机变量生成方法
scc的博客
10-26 5395
随机数生成方法1 引言2 随机数生成方法-逆变换方法(Inverse Transform Method) （本篇博客是自己根据导师的统计计算课程的课件自行整理可得，如有错误，请联系我，此笔记仅用于学习，如需转载，请注明来源，谢谢！） 1 引言 计算统计的基本工具就是要从指定的概率分布中生成随机变量，最基本最重要的是均匀分布的伪随机数生成器，其他概率分布的随机数生成方法大都依赖于均匀分布随机数生成器。 在R语言中，均匀分布的随机数的生成为： runif(n) #generate a vector of siz
Matlab求离散数据的一阶导
qq_43738047的博客
10-18 5601
Matlab求离散数据的一阶导 本测试使用Matlab-2019b–***-diff函数*** y=[离散数据，我的离散数据共1937个数，要加入的时候要调控你的x和dx的值]; x=0:0.01:19.36; % 设置X长度，0-19.36，步长为0.01 dy = diff(y)./diff(x); % 求取离散数据的导，会比y少一位 dx = 0:0.01:19.35; % 设置求导后dx的长度，会比x少一位 subplot(2,1,1); plot(x,y);
clr错误 8004005程序将立即终止_Matlab 的程序调试及其他功能
weixin_39630048的博客
11-21 452
（建议阅读原文）预备知识　Matlab 的函数程序调试 图 1：在行首设置 Breakpoint 　　 若要调试程序，可选择一行代码并单击该行前面的横线，这时会出现红色圆点 Breakpoint （图 1 ），程序运行到 Breakpoint 会暂停． 　　 此时要查看变量情况，可通过 Workspace 查看各个变量的情况，也可用光标悬停在某个变量上．还可以用 Command Window ...
快速切比雪夫微分：计算切比雪夫网格上的数值导数。-matlab开发
06-01
此脚本在 O(N*log(N)) 时间内对 Chebyshev-Gauss-Lobatto 网格上的函数进行数值微分。 等效于计算 D*x，其中 D 是切比雪夫搭配矩阵，但对于大向量 x 更快。
Higher-Order_Spectral_Analysis_toolbox.rar_matlab 分析_spectral_to
07-14
《深入探索高阶谱分析工具箱在MATLAB中的应用》 在现代信号处理领域，高阶谱分析（Higher-Order Spectral Analysis, HOSA）作为一种强大的技术，被广泛应用于非高斯信号的特征提取、噪声检测以及系统诊断等多个方面...
gpops.zip_matlab gpops_spectral method_zzzyyykkk_伪谱法_伪谱法 matlab
07-14
伪谱法安装程序，直接在MATLAB环境下运行安装文件，随后即可使用，效果不错的一种方法，希望有用。
Log_Spectral_Distance.rar_LSD_spectral distance_倒谱距离_对数谱距离lsd_对数
07-14
这里提到的"Log_Spectral_Distance"（LSD），即对数谱距离，是一种广泛使用的工具，用于衡量这种差异。本文将深入探讨对数谱距离的概念、计算方法以及其在实际应用中的重要性。 对数谱距离（LSD）是基于频谱分析的...
Classical_spectral_estimation.zip_matlab例程_matlab_
08-12
本示例主要探讨了两种经典的谱估计方法：自回归（AR）模型和自回归移动平均（ARMA）模型，并通过MATLAB编程实现进行比较。我们将详细解析这两种模型的原理、区别以及在MATLAB中的应用。 自回归（AR）模型是一种线性...
矩阵求导手册Matrix Differentiation
01-19
矩阵求导常见的公式，值得收藏，以备未来查询之用！！
MATLAB——求解矩形域内的波动方程
08-27
本代码主要利用MATLAB工具实现MATLAB——求解矩形域内的波动方程，简单明了，易于理解
gonglvpumidu.rar_matlab_power spectral_功率谱密度_功率谱密度midu_高斯功率谱
09-24
总结来说，`gonglvpumidu.rar_matlab_power spectral_功率谱密度_功率谱密度midu_高斯功率谱`这个资源可能是用于演示或教学如何在MATLAB中计算高斯信号的功率谱密度，尤其是使用Welch方法。通过对高斯信号的功率谱...
matlab db函数_MATLAB仿真阵列天线切比雪夫综合法（附代码）
weixin_39743369的博客
11-28 6601
来源：cnblogs在《自适应天线与相控阵》这门课中，了解到了关于理想低副瓣阵列设计的一些方法，其中切比雪夫等副瓣阵列设计方法是一种基础的方法，故将其设计流程写成maltab程序供以后学习使用。在此分享一下。此方法全称为道尔夫-切比雪夫综合法，简称为切比雪夫综合法，是一种工程实际中常用的可控制副瓣电平的阵列天线综合方法。切比雪夫阵列的特点是：(1)等副瓣电平；(2)在相同副瓣电平和相同...
matlab 切比雪夫矩源代码,matlab实现快速切比雪夫二阶微分
weixin_30140317的博客
03-18 423
Fast Chebyshev second order differentiationfchd2(V) computes the first derivative of the data in V located along the N+1 Chebyshev–Gauss–Lobatto points cos(pi*(0:N)/N).Example 1:Use fchd2 to find the ...
MATLAB中用差分法求解函数 or 离散数据的一阶导数和二阶导数
热门推荐
han____shuai的专栏
02-23 2万+
matlab在计算数值求导方面有很多函数。下面我来不用这些函数，介绍简单的函数求一阶导数二阶导数的差分方法以及其MATLAB实现。 工具/原料 matlab软件 一组数据，程序中已经准备好了 方法/步骤 求解一阶导数的公式：y'=[y(x0+h)-y(x0-h)]/(2h); 求解二阶导数的公式：y''=[y(x0+h
MATLAB对plot曲线进行求导
csdnpliu2021的博客
08-12 8591
t=0:0.01:5; v=exp(-t)+sin(t) figure plot(t,v); h=get(gca,'children') x = get(h,'XData'); y = get(h,'YData'); dy = diff(y) ./ diff(x); hold on plot(x(1:end-1),dy,'r') legend('原始图像','求导图像')
matlab 离散求偏导,matlab怎么求离散曲面的偏导数
weixin_39779032的博客
03-16 966
我想求以下坐标点构成的曲面的偏导数，不是先拟合出函数再求，而是对离散点用差分法求，不知道怎么操作，希望能得到解答x=linspace(4.97,7.47,12);y=linspace(1.34,3.84,12);z=[0.798244 0.901491 1.017275 1.143169 1.278748 1.428146 ...
讲一下tensorflow代码转为pytorch代码。def simulated_spectral9_camera_spectral_response_function(hyper_spectral_image,masked_response_function): masked_response_function = tf.cast(masked_response_function, dtype=hyper_spectral_image.dtype) ##9,31 hyper_spectral_image = hyper_spectral_image[:,1:-1, 1:-1,:] #裁切中心部分 batch_size,h,w,c = hyper_spectral_image.shape response3x3 = tf.reshape(masked_response_function,[3,3,31]) #padding到 h,w responsehxw = tf.tile(response3x3,[h//3,w//3,1]) response_img = hyper_spectral_image*responsehxw response_img = tf.reduce_sum(response_img,axis=-1)/ tf.reduce_sum(responsehxw,axis=-1) #,keepdims=True # 预期 shape: (batch, height, width, 1) return response_img
最新发布
06-01
def simulated_spectral9_camera_spectral_response_function(hyper_spectral_image, masked_response_function): masked_response_function = masked_response_function.to(hyper_spectral_image.dtype) hyper_...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值