0

我在 matlab 文件中看到了这个。频率.m 文件

   h = dividenowarn(fft([b zeros(1,s*nfft-nb)]),...
                   fft([a zeros(1,s*nfft-na)])).';

例子： 

x = fft([1.5,0,1,0,0,0,1,3]')

x =

   6.5000          
   3.6213 + 2.1213i
  -0.5000 + 3.0000i
  -0.6213 + 2.1213i
   0.5000          
  -0.6213 - 2.1213i
  -0.5000 - 3.0000i
   3.6213 - 2.1213i

现在

y = fft([1,1,2,3,1,0,9,3]')

y =

20.0000          
 0.7071 + 6.2929i
-9.0000 + 5.0000i
-0.7071 - 7.7071i
 6.0000          
-0.7071 + 7.7071i
-9.0000 - 5.0000i
 0.7071 - 6.2929i

fft 并不重要，我需要如何执行此操作..

z = (x./y)

z =

   0.3250          
   0.3968 - 0.5309i
   0.1840 - 0.2311i
  -0.2656 - 0.1050i
   0.0833          
  -0.2656 + 0.1050i
   0.1840 + 0.2311i
   0.3968 + 0.5309i

我需要一个算法（没有 matlab 代码），我需要一些 Java 的东西，或者一步一步的微积分......

» a

a =

   1.0000 + 2.0000i   3.0000 + 4.0000i   5.0000 + 6.0000i        0          

» b

b =

   5.0000 + 2.0000i   1.0000 - 2.0000i        0                  0          

» c = a./b
Warning: Divide by zero.

c =

   0.3103 + 0.2759i  -1.0000 + 2.0000i      Inf +    Infi      NaN -    NaNi

»
4

1 回答 1

1

./运算符执行逐元素除法。您可以从除号之前的元素判断它是一个元素运算符.。这意味着结果将是一个向量，其中包含使用规则获得的元素x[i] / y[i]。

如果您想在 Java 中执行此操作，您将需要自己编写除法代码来实现自己的Complex数字类，或者您可以使用Apache commons 数学类。 Complex

假设您使用 apache commons，Java 中的元素划分将如下所示：

List<Complex> elementWiseDivision(List<Complex> dividend, List<Complex> divisor) 
{
  if (dividend.size() != divisor.size()) 
  { 
    throw new IllegalArgumentException("Must have same size"); 
  }
  List<Complex> result = new ArrayList<>();
  // using iterators to get O(n) with both LinkedList and ArrayList inputs
  for (Iterator<Complex> xit = dividend.iterator(), yit = divisor.iterator(); xit.hasNext();) 
  {
    result.add(xit.next().divide(yit.next()));
  }
  return result;
}
于 2015-01-08T04:48:25.367 回答