Matlab для c, реализующий цикл for

Я пытаюсь реализовать это

Fs = 2000;            % Sampling frequency                    
T = 1/Fs;             % Sampling period       
L = 4000;             % Length of signal
t = (0:L-1)*T;          %time samples
A=[16 18 28 33 38 41 59 43 40 58];% Amplitude of noise in dBA (refer to the measurements from consultant)
forigin=[25 31.5 40 50 63 80 100 125 160 200];% Frequency of the noise components (refer to the measurements from consultant)
S=zeros(1,length(t));% This will be the signal representing transformer noise, in your case, it will be the sound created from exciter 
%V=zeros(1,L);%not needed, you can remove this
for k=1:length(A)
S = S+10^(A(k)/20)*sin(2*pi*forigin(k)*t);%Creating the transformer noise from the amplitude and frequency components, in your case, it will be the sound created from exciter

end

Я получил это до сих пор, не похоже, что k увеличивается через массив forigin.

#define _USE_MATH_DEFINES  // to force M_PI to exist
#include <math.h>          //for pow, sin and M_PI
#define ASIZE 10   //for A 
#define TSIZE 40   //or 4000
int Fs = 2000;
int L = 40;
double A[ASIZE] = {16, 18, 28, 33, 38, 41, 59, 43, 40, 58};
double forigin[ASIZE] = {25.0, 31.5, 40.0, 50.0, 63.0, 80.0, 100.0, 125.0, 160.0, 200.0};
double S[TSIZE];
double t[TSIZE];
double T = 1/(double)Fs;
int i,k,  j;
//int k = 10;
 double  apow = pow(10, A[k]/20);
void setup() {

  Serial.begin(9600);
}

void loop()
{
   for ( i = 0; i < TSIZE ; i++ )  /* initialize elements of array t to 0 */  
       {
      t[i] = (double) i*T; /* set element at location i to i*T */
       }
      for (j = 0; j < TSIZE; j++ ) /* output each array element's value */
       {

//      Serial.print(t[j],15);
//     Serial.println("");
//     delay(500);
       } 

//initialise S to 0
for (i = 0; i < TSIZE; ++i) //S[i] = 0;

//implement the for k=1:length(A)

  for (k = 0; i < ASIZE; ++k){
   //devectorise the S calculation
   for (i = 0; i < TSIZE; ++i){
       S[i] = S[i] + pow(10, A[k]/20) * sin(2*M_PI*forigin[k]*t[i]);
       Serial.print(S[i]);
     Serial.println("");
     delay(500);

Может ли кто-нибудь помочь. Matlab выпустил

Фактические значения S в Matlab: 0 899.475071494380 1570.91017541897 1861.36356724185 1740.88440906200 1306.16803613629 739.441548916790 238.639601470524 -54.6370978848435 -102.514143930475

Как видите, я не получаю ничего подобного.

0
0
28
1

Ответы 1

У меня сегодня мало времени, поэтому я написал ответ для обычного компьютера, а не для встроенной системы (например, Arduino, я думаю, вы работаете с чем-то подобным).

Я думаю, что у вас неправильное приращение, а сам код довольно беспорядочный, поэтому я решил переписать его с нуля (ошибки в циклах for, циклы forr, которые плохо прокомментированы, плохие отступы ...):

/** @file: test.c */
#include <stdio.h>
#include <math.h>
#include <stddef.h>

#define T_LENGTH 4000
#define A_LENGTH 10

int main() {
  const double Fs = 2000.0;
  const double T = 1/Fs;

  const double A[A_LENGTH] = {
    16, 18, 28, 33, 38, 41, 59, 43, 40, 58
  };

  const double forigin[A_LENGTH] = {
    25, 31.5, 40, 50, 63, 80, 100, 125, 160, 200
  };

  double t[T_LENGTH];
  double S[T_LENGTH] = {0};

  for  (size_t i = 0; i < T_LENGTH; i++) 
    t[i] = (double)(i) * T;

  for (size_t i = 0; i < A_LENGTH; i++) { // it is i++ not ++i
    for (size_t j = 0; j < T_LENGTH; j++) {
      S[j] += pow(10, A[i] / 20) * sin(2 * M_PI * forigin[i] * t[j]);
    }
  }

  // Printing out the results
  for (size_t j = 0; j < T_LENGTH; j++)
    printf("%f\n", S[j]);

  return 0;
}

Вы можете скомпилировать этот код с помощью:

gcc test.c -lm -o test

и запустите его с помощью:

./test > output.txt

чтобы поместить результат в файл и сравнить его с вашим результатом в Matlab:

load output.txt -ascii
norm(S' - output)
% something in the order of 1e-7

Убедившись, что все в порядке, вы можете поместить его в Arduino IDE.

Другие вопросы по теме