【MATLAB库函数系列】resample(重采样函数)的C语言实现【姊妹篇1基于std::vector】

已于 2023-02-12 21:58:51 修改
阅读量5.3k 收藏 17
点赞数 5
分类专栏： VC算法与小工具实现 文章标签： resample matlab 重采样 c upfirdn
于 2020-09-20 17:14:12 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/wlwdecs_dn/article/details/108694980
版权

resample的原理

MATLAB中resample函数最简单的参数如下：
y = resample(x,p,q)
以p/q乘以原始采样速率对输入序列x重新采样。resample在x上应用抗锯齿FIR低通滤波器，补偿滤波器带来的延迟。
根据以下2篇文章
MATLAB库函数upfirdn(分数倍采样率变换)的C语言实现（采用了STL::vector）
MATLAB库函数firls(最小二乘线性相位FIR滤波器设计)的C语言实现
可以知道，resample的核心还是分数倍采样率变换，而分数倍采样率变换的核心是：插零=>低通滤波=>抽取

常量和结构体定义

常量和结构体定义的头文件ConstParam.h在下面文章中
MATLAB库函数hilbert(希尔伯特变换)的C语言实现(FFT采用FFTW库)

resample的C语言实现

bessel.h和upfirdn.h请参考之前的博客，下面将核心代码进行注释。
第一类修正贝塞尔函数的C语言实现


                

                
                
        

        

    




        

            

                

                

                了解本专栏
                
                
                  
                订阅专栏 解锁全文
                 
            

        

    

      

        

            

              
                
                  codersnote
                
              
            

            

                关注
              
            

        

        

          
    
            
  • 
              
                
                
                
                
                    5
                
              
              
    点赞
    
            
    • 
            
  • 
              
                
                
                
              
              
    踩
    
            
    • 
            
  • 
              
                
                
                
                
                    17
                
              
              
    
                
    
                  
                    收藏
                  
                
    
              
    
              
    
                
    
                  觉得还不错?
                  
                    一键收藏
                  
                 
                
    
              
    
            
    • 
                
  • 
                  
                    打赏
                    
                  
                  
    打赏
    
                
    • 
          
  • 
            
    
              
              
            
    
              
              
              
                  17
              
            
            
    评论
    
          
    • 
          
  • 
          
    • 
          
  • 
            
              
            
              
    
            
              
                
            
    
          
    • 
        

      

      

            

                专栏目录
          

                

                        订阅专栏
                










                



	

		

			

				
					
【MATLAB库函数系列】resample(重采样函数)的C语言实现【姊妹篇2纯C语言实现】

				
			

			

				

					业精于勤，荒于嬉
				

				

					09-28
					
					2460
					
				

			

		

		

			
				
常量和结构体定义的头文件ConstParam.h在下面文章中MATLAB库函数hilbert(希尔伯特变换)的C语言实现(FFT采用FFTW库)

			
		

	



                

              
                



	

		

			

				
					
分数倍转采样算法与多相滤波实现

				
			

			

				

					qq_33467472的博客
				

				

					11-20
					
					2205
					
				

			

		

		

			
				
多相滤波是利用内插、滤波和抽取运算在一定条件下的可交换性，保证输出结果相同的前提下减少运算次数，提高滤波效率的一种实现方式。

			
		

	



                


  
              

                


	

		

			

				
					
c语言实现线性重采样

				
			

			

				

					08-13
				

			

		

		

			
				
线性重采样demo代码，可实现8K/16K转32K/48K，16bit精度转32bit精度，可使用线性算法也可简单复制，自测有效，音频测试可采用。

			
		

	





	

		

			

				
					
MATLAB转C语言的步骤与解决方案

				
			

			

				

					weixin_44982036的博客
				

				

					12-11
					
					1万+
					
				

			

		

		

			
				
MATLAB转C语言的步骤
步骤
我们首先以求三阶矩阵的逆为例，首先在matlab中编写的代码如下：
  function[res] = inverse(src)
  res = inv(src)
  end

此为新建的一个函数，函数的保存名称为 inverse.m ,需要注意的是这个.m文件必须是函数文件，而不能是脚本文件，而且文件的名称和函数的名称必须一致。
在matlab中输入一个三阶矩阵与输出的结果如下所示：
`
在本文的末尾会验证在VS2017C语言代码的运行结果。
进入主题，在已经写好了函数的

			
		

	



		

			
		



	

		

			

				
					
resample函数的C语言实现

				
			

			

				

					10-23
				

			

		

		

			
				
由于matlab程序将各信号处理函数封装，无法看到原始代码，这里是其中resample函数的C语言实现源码，可实现任意整数倍上采及下采

			
		

	





	

		

			

				
					
重采样：resample

				
			

			

				

					lbywyj的专栏
				

				

					05-24
					
					1414
					
				

			

		

		

			
				
resample(rule, how=None, axis=0, fill_method=None, closed=None, label=None, convention='start', kind=None, loffset=None, limit=None, base=0, on=None, level=None)






			
		

	





	

		

			

				
					
matlab中resample重采样函数

				
			

			

				

					weixin_46463398的博客
				

				

					12-18
					
					2万+
					
				

			

		

		

			
				
resample函数
对时间序列进行重采样
格式：
1、y=resample(x,p,q)
x–待重采样的时间序列；
p–重采样之后目标频率；
q–待重采样的时间序列频率
采用多相滤波器对时间序列进行重采样，得到的序列y的长度为原来的序列的长度的p/q倍，p和q都为正整数。此时，默认地采用使用FIR方法设计的抗混叠的低通滤波器
2、y=resample(x,p,q,n)
n–滤波器长度与n成正比；
采用chebyshevIIR型低通滤波器对时间序列进行重采样，滤波器的长度与n成比例，n缺省值为10.
3、y

			
		

	





	

		

			

				
					
通信算法之126：数字信号处理-采样库函数

				
			

			

				

					leegang12的专栏
				

				

					04-01
					
					561
					
				

			

		

		

			
				
MATLAB中resample函数使用一.

			
		

	





	

		

			

				
					
C语言库函数速查手册 C语言库函数速查手册

				
			

			

				

					06-16
				

			

		

		

			
				
C语言库函数速查手册C语言库函数速查手册C语言库函数速查手册C语言库函数速查手册C语言库函数速查手册C语言库函数速查手册C语言库函数速查手册C语言库函数速查手册C语言库函数速查手册C语言库函数速查手册C语言...

			
		

	





	

		

			

				
					
resample, 采样率转换C代码

				
			

			

				

					10-11
				

			

		

		

			
				
可对音频或数字信号进行重采样， 提供了抽取、插值、分数倍采样率转换的代码

代码在linux上实现，在windows上只需要建工程编译即可

更多代码，可访问：
http://code.google.com/p/falab

			
		

	





	

		

			

				
					
resample.cpp

				
			

			

				

					04-15
				

			

		

		

			
				
matlab 的resample函数的C语言版本，对一维信号进行重采样，与matlab结果是一致的。

			
		

	





	

		

			

				
					
C库函数手册_C函数库_函数库_C语言库函数手册_C++库函数手册_语言库_

				
			

			

				

					09-30
				

			

		

		

			
				
C库函数手册是程序员在使用C语言进行编程时的重要参考资料，它详尽地列出了C语言标准库中提供的各种函数，这些函数涵盖了输入输出、字符串处理、数学运算、内存管理等多个方面。C语言库函数手册不仅是初学者学习...

			
		

	





	

		

			

				
					
传递函数的C语言实现.zip_c 传递函数_substancef7x_传递函数_传递函数C语言实现_传递函数c语言

				
			

			

				

					07-14
				

			

		

		

			
				
在本文中，我们将深入探讨如何使用C语言来实现传递函数，这是控制理论中的一个重要概念，特别是在数字信号处理和系统分析中。传递函数是描述线性时不变系统动态特性的数学模型，它以频率域的形式表示输入与输出之间...

			
		

	





	

		

			

				
					
FFT快速算法的MATLAB示例：可以提供C语言的实现思路.zip

				
			

			

				

					08-20
				

			

		

		

			
				
对于想要将MATLAB代码转换为C语言实现的学生，需要注意的是，MATLAB的内置函数在C语言中可能没有直接对应的库函数。因此，需要找到相应的算法实现，比如FFTW库，这是一个高效的C语言编写的FFT库，可以作为C语言实现...

			
		

	





	

		

			

				
					
c语言函数库调用.rar_C函数_C函数库_c语言函数库_c语言库函数_库函数

				
			

			

				

					09-21
				

			

		

		

			
				
C语言函数库，也称为标准库，是C编程语言的核心组成部分，它提供了大量的预定义函数，使得程序员可以方便地执行各种任务，如输入输出、字符串处理、数学计算、内存管理等。C语言库函数是C语言编程的重要工具，极大地...

			
		

	





	

		

			

				
					
手撕重采样，考虑C的实现方式

					
最新发布

				
			

			

				

					weixin_40194697的博客
				

				

					01-24
					
					1008
					
				

			

		

		

			
				
3、具体来说，插值，频谱变窄，即信号频带压缩了，如果不做处理，信号会包含带宽以外的噪声，所以需要做低通只滤出变窄的信号频带，去掉噪声。考查原始信号时间点和重采样后信号的时间点，一般来说，重采样信号的一个时间点，是处在原始信号的两个信号点之间（若重合，则直接不用计算），根据此时间点到原信号两个信号点的距离，用原信号这两个信号点插值得到重采样信号此点的信号幅值。1、对于原采样率fs，需要重采样到fs1，一般fs和fs1都是整数哈，则先找fs和fs1的最小公倍数，设为m，设m/fs=M，m/fs1 = L。

			
		

	





	

		

			

				
					
C

				
			

			

				

					qq_44570750的博客
				

				

					03-02
					
					174
					
				

			

		

		

			
				
循环
#include<stdio.h>
int main ()
{
	int i,j;
	for(i=0;i<5;i++)
	  {
	    for(j=0;j<i+1;j++)
	       printf("*");
	    printf("\n");}
	return 0;
 } 


注意转义字符
//注意转义字符 
#include<stdio.h&g...

			
		

	





	

		

			

				
					
基于C语言实现PCM音频流或音频文件重采样（48K到16K）

					
热门推荐

				
			

			

				

					a18359388932的博客
				

				

					08-12
					
					23万+
					
				

			

		

		

			
				
基于C语言实现PCM音频流或音频文件重采样（48K到16K）
由于云厂商SDK需要的音频采样率是48K的，而SFU回调上来的流是48K的，所以我们还需要对PCM音频数据进行重采样处理。
​    转换的原理比较简单，从48KHz降到16KHz，降了3倍，也就是说在同一时间的单位区间内，48KHz采样了3个点，而16KHz只采样了1个点，即从48KHz的 PCM流中每读取3个数据，就要根据这3个数据去推算得到1个数据，而这个数据对应的就是16KHz PCM流中的一个数据。

我们有很多种方法可以由3个数据得到

			
		

	





	

		

			

				
					
用java实现C++中得库函数std::arg(std::complex)

				
			

			

				

					06-09
				

			

		

		

			
				
在Java中，可以使用 `java.util.function.Function` 接口来实现类似于 C++ 中的 `std::arg(std::complex)` 函数。

首先，我们可以定义一个 `Complex` 类表示复数，包含实部和虚部：

```java
class Complex {
    private double real;
    private double imaginary;

    public Complex(double real, double imaginary) {
        this.real = real;
        this.imaginary = imaginary;
    }

    public double getReal() {
        return this.real;
    }

    public double getImaginary() {
        return this.imaginary;
    }
}
```

然后，我们可以定义一个 `arg` 函数，该函数接受一个 `Complex` 对象作为参数，并返回该复数的幅角（弧度制）：

```java
import java.util.function.Function;

public class Main {
    public static void main(String[] args) {
        // 创建一个复数对象
        Complex z = new Complex(1.0, 1.0);

        // 调用 arg 函数获取幅角
        double arg = arg(z);

        // 输出结果
        System.out.println(arg);
    }

    public static double arg(Complex z) {
        // 使用 Function 接口计算幅角
        Function<Complex, Double> argFunc = (c) -> Math.atan2(c.getImaginary(), c.getReal());

        // 调用 apply 方法获取结果
        return argFunc.apply(z);
    }
}
```

在 `arg` 函数中，我们定义了一个 `Function` 对象 `argFunc`，该对象接受一个 `Complex` 对象作为参数，并使用 `Math.atan2()` 方法计算幅角。然后，我们调用 `apply()` 方法，将复数对象传递给 `argFunc`，得到幅角的结果。

注意，`Math.atan2()` 方法返回的是弧度制的幅角，与 C++ 中的 `std::arg()` 函数一致。如果需要得到角度制的幅角，可以将结果乘以 `180.0/Math.PI`。

			
		

	



              



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

  

      
    
          
  • 
              
    
                  
                  
              
    
              
    非常没帮助
    
          
    • 
          
  • 
              
    
                  
                  
              
    
              
    没帮助
    
          
    • 
          
  • 
              
    
                  
                  
              
    
              
    一般
    
          
    • 
          
  • 
              
    
                  
                  
              
    
              
    有帮助
    
          
    • 
          
  • 
              
    
                  
                  
              
    
              
    非常有帮助
    
          
    • 
      

      

      

          
          提交
      

      

  




          




        



    





    



      

      

        
      

      





	

		
评论 17

	

	
	




  

    

      添加红包
      
    

    

      

        
        

          
          
        

        
请填写红包祝福语或标题

      

      

        
        

          
          个
        

        
红包个数最小为10个

      

      

        
        

          
          元
        

        
红包金额最低5元

      

      

        
        

          当前余额3.43元
          前往充值 >
        

      

      

        

          需支付：10.00元
        

        
        
      

    

  





  

    
    
  

  

    
    
  








  

    

      

        

          

            
            
成就一亿技术人!

          

          

        

        

          

          

            领取后你会自动成为博主和红包主的粉丝
            规则
          

        

      

      

        

          

            
              
            
            

              
hope_wisdom
 发出的红包
            

          

        

        

          

          

          

        

      

    

    

  





	
打赏作者

	

		

		
			
		
		

		

			
codersnote

			
对学生党 赞赏是鼓励也是鞭策！

		

	

	

			

	

	

    ¥1
    ¥2
    ¥4
    ¥6
    ¥10
    ¥20
	

	

		

			

				

					扫码支付：¥1
				

				

					

					
					获取中
					

				

				

					
					
					扫码支付
				

			

		

		

			
您的余额不足，请更换扫码支付或充值

			
打赏作者

		

	



      
      

      
实付元

      
使用余额支付

      

        

          

            
            点击重新获取
          

        

        
扫码支付

      

      

        
          
            
          
          
        
      

      

        
        钱包余额
          0
          

            
            

              

                
抵扣说明：

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

              

            

          

      

      余额充值