Перевод кода Matlab в Си

@nkment · Регистрация: 11.12.2009

Студворк — интернет-сервис помощи студентам

Есть ПРАВИЛЬНЫЙ алгоритм в MATLABe:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
function PSHz = PSHzformule(Fn, Fv, P_1, a)
% ВЫЧИСЛЕНИЕ ИНТЕГРАЛЬНОГО УРОВНЯ ДАВЛЕНИЯ, Па 
Fn=1000*Fn;     % НИЖНЯЯ ЧАСТОТА СПЕКТРА,Гц
Fv=1000*Fv;     % ВЕРХНЯЯ ЧАСТОТА СПЕКТРА,Гц 
a=-a;           % ЗАКОН СПАДА СПЕКТРАЛЬНЫХ СОСТАВЛЯЮЩИХ, дБ/окт
n=Fv-Fn;
Sum=0;
Pn=P_1/10^(a*log2(1000/Fn)/20);
for i=1:n
   S=(Pn*10^(a*log2((Fn+i)/Fn)/20))^2;
   Sum=Sum+S;
end
PSHz=sqrt(Sum); % интегральная плотность сплошной части
return

написал тоже самое в Си

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
float PSHz (double Fn, double Fv, double P_1k, double a)
{
  double Pi;
  Fn*=1000;
  Fv*=1000;
  a=-a;
  double n=Fv-Fn;
  double Psum=0;
  double temp=a/20*Log2(1000/Fn);
  double Pn=P_1k/pow(10,temp);
  for (int i=1; i<n; i++)
  {
    temp=a/20*Log2((Fn+i)/Fn);
    Pi=Pn/pow(10,temp);
    Pi=pow(Pi,2);
    Psum=Psum+Pi;
  }
  Psum=sqrt(Psum);
  return Psum;
 
}

прогнал оба алгоритма для
PSHz (0.1, 10, 1, 6)

в результате получается:

S1=96,5046, S2=94.6271, S3=92.8048, S4=91.0347 => PSHz=98.7935
P1=100.4084, P2=102.3, P3=104., P4=106. => Psum=56446 => ПОЛНЫЙ БРЕД

помогите разобраться в чем проблема

@Somebody · 18.12.2009, 21:13

Code
1
2
3
for i=1:n
 
for (int i=1; i<n; i++)

Наверное, условие i<=n.

Code
1
2
3
4
5
S=(Pn*10^(a*log2((Fn+i)/Fn)/20))^2;
 
temp=a/20*Log2((Fn+i)/Fn);
Pi=Pn/pow(10,temp);
Pi=pow(Pi,2);

В первом случае умножение на десять в степени, во втором - деление.

Новые блоги и статьи Все статьи Все блоги /
Изучаю kubernetes lagorue 13.01.2026 А пригодятся-ли мне знания kubernetes в России?	Сукцессия микоризы: основная теория в виде двух уравнений. anaschu 11.01.2026 https:/ / rutube. ru/ video/ 7a537f578d808e67a3c6fd818a44a5c4/	WordPad для Windows 11 Jel 10.01.2026 WordPad для Windows 11 — это приложение, которое восстанавливает классический текстовый редактор WordPad в операционной системе Windows 11. После того как Microsoft исключила WordPad из. . .	Classic Notepad for Windows 11 Jel 10.01.2026 Old Classic Notepad for Windows 11 Приложение для Windows 11, позволяющее пользователям вернуть классическую версию текстового редактора «Блокнот» из Windows 10. Программа предоставляет более. . .
Почему дизайн решает? Neotwalker 09.01.2026 В современном мире, где конкуренция за внимание потребителя достигла пика, дизайн становится мощным инструментом для успеха бренда. Это не просто красивый внешний вид продукта или сайта — это. . .	Модель микоризы: классовый агентный подход 3 anaschu 06.01.2026 aa0a7f55b50dd51c5ec569d2d10c54f6/ O1rJuneU_ls https:/ / vkvideo. ru/ video-115721503_456239114	Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ФедосеевПавел 06.01.2026 Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .	Модель микоризы: классовый агентный подход 2 anaschu 06.01.2026 репозиторий https:/ / github. com/ shumilovas/ fungi ветка по-частям. коммит Create переделка под биомассу. txt вход sc, но sm считается внутри мицелия. кстати, обьем тоже должен там считаться. . . .