Каменная лампа

3 Страниц
←
1
2
3

Вы не можете создать новую тему
Вы не можете ответить в тему

Каменная лампа или о том, как из ОУ сделать лампу

#61 Маклауд

Злой модер

Отправлено 16 мая 2010 - 19:49

Rst7, я щас писать в линию не могу, звуковуха косячит. Ещё какой-то писк пробивается. Наверное от моника.
Тут Морячок оговорился, что периодически записывается в комп)
Этот семпл покруче будет!

"...Металлисты - это самый развитой и передовой класс, и никто не может отрицать, что это и есть передовой отряд всего пролетариата." (В.И. Ленин, "Полное Собрание Сочинений", том 24)

Наверх of the page up there ^

#62 Rst7

Отправлено 16 мая 2010 - 19:53

Цитата

Этот семпл покруче будет!

Который take3?

"Практика выше (теоретического) познания, ибо она имеет не только достоинство всеобщности, но и непосредственной действительности." - В.И. Ленин

Наверх of the page up there ^

#63 Маклауд

Злой модер

Отправлено 16 мая 2010 - 19:56

Rst7, я пост написал раньше, чем take 3 увидел)

Наверх of the page up there ^

#64 Rst7

Отправлено 16 мая 2010 - 21:48

Маклауд, ну так что, take3 назначаем годным вариантом?

Пока суть да дело, попропускал другие семплы из гитаррига через девайс (это, Аллах акбар, намного быстрее, чем симулить). Игрался только регулятором гейна.

Пара "чистых":

clean1.mp3 (980,72К)
Количество загрузок:815

clean2.mp3 (921,53К)
Количество загрузок:836
Овердрайв несильный:

overdrive.mp3 (493,98К)
Количество загрузок:1000
Посильнее:

overdrive2.mp3 (851,13К)
Количество загрузок:723

Наверх of the page up there ^

#65 Dr_Zlo

Отправлено 17 мая 2010 - 01:14

Маклауд, Если писать семплы на гитарке с синглами, в комп - это пистец, фон аццкий от моника. Я пишу либо вырубая моник либо хамбом либо двумя противофазно включенными синглами. (Гитарка - фендер, бридж хамб + сингл + нек сингл)

И не забывай, юный падаван, страх, гнев, и гитары FENDER ведут тебя к темной стороне силы...

Наверх of the page up there ^

#66 Rst7

Отправлено 18 мая 2010 - 11:00

Еще немножко поковырялся и, видимо, пока хватит:

take5.mp3 (1,1МБ)
Количество загрузок:743

Подсказали место, где можно разжиться подходящими семплами. Кое-какие покормил девайсу. Вот результаты -

samples.zip (2,92МБ)
Количество загрузок:435

Наверх of the page up there ^

#67 Маклауд

Злой модер

Отправлено 18 мая 2010 - 17:51

Rst7, хорошие семплы получились. Последний скорее на фузз смахивает, но звучит достойно.

Наверх of the page up there ^

#68 Rst7

Отправлено 18 мая 2010 - 18:48

Цитата

Последний скорее на фузз смахивает, но звучит достойно.

Эээ, их там уже дофига, последний - это какой конкретно?

Наверх of the page up there ^

#69 Маклауд

Злой модер

Отправлено 18 мая 2010 - 18:50

8.мрз

Наверх of the page up there ^

#70 Rst7

Отправлено 19 мая 2010 - 20:29

Цитата

хорошие семплы получились

Да, кое-кого из моих знакомых они тоже порадовали. Ну что, хотите конкретики по технике?

Сообщение отредактировал Rst7: 19 мая 2010 - 20:33

Наверх of the page up there ^

#71 SnigE

Отправлено 19 мая 2010 - 21:21

звук однозначно хороший, мне семпл 6 понравился. Хотелось бы конкретике по технике)

=*\*=*\*=*\*=*\*=*\*=*\*=*\*=*\*=*\*=*\*=*\ *=*\*=\*=*\=*\=*\

Наверх of the page up there ^

#72 Rst7

Отправлено 20 мая 2010 - 11:39

Цитата

Хотелось бы конкретике по технике)

Конкретика веселая :clapping:

Спойлер

//---------------------------------------------------------------------------
#include <stdio.h>
#include <stdlib.h>
#include <windows.h>
#include <winbase.h>
#include <string.h>
#include <math.h>

#pragma hdrstop

//---------------------------------------------------------------------------

typedef struct
{
  long Y;
  unsigned long TAU;
}HPF;

typedef struct
{
  long Y;
  unsigned long TAU;
}LPF;

typedef struct
{
  unsigned long GAIN;
  unsigned long ODDEVEN;
  LPF zero;
}SQ_DIST;

typedef struct
{
  signed long a0;
  signed long a1;
  signed long a2;
  signed long b1;
  signed long b2;
  signed long mem0;
  signed long mem1;
  signed long mem2;
  signed long mem3;
}BIQUAD;

signed long DebugSample;

static signed long LowPassFilter(signed long s, LPF *f)
{
  signed long long acc;
  acc=f->Y-s;
  acc*=f->TAU;
  if (acc>0)
  {
    acc+=0x80000000;
    acc>>=32;
  }
  else
  {
    acc=-acc;
    acc+=0x80000000;
    acc>>=32;
    acc=-acc;
  }
  acc=f->Y-acc;
  if (acc>0x7FFFFFFF) acc=0x7FFFFFFF;
  if (acc<-0x7FFFFFFF) acc=-0x7FFFFFFF;
  f->Y=acc;
  return acc;
}

static unsigned long int_sqrt(unsigned long v)
{
  unsigned int mask=0x00008000;
  unsigned int result=0;
  do
  {
    result|=mask;
    if (result*result>v) result&=~mask;
  }
  while(mask>>=1);
  return result;
}

static signed long HighPassFilter(signed long s, HPF *f)
{
  signed long long acc;
  acc=f->Y-s;
  acc*=f->TAU;
  acc>>=32;
  acc=f->Y-acc;
  if (acc>0x7FFFFFFF) acc=0x7FFFFFFF;
  if (acc<-0x7FFFFFFF) acc=-0x7FFFFFFF;
  f->Y=acc;
  return s-f->Y;
}

static signed long SquareDist(signed long s, SQ_DIST *d)
{
  signed long long acc;
  signed long z;
  if (s>d->zero.Y)
  {
    //Пиковый детектор
    z=LowPassFilter(s*3,&d->zero);
  }
  else
  {
    z=LowPassFilter(s,&d->zero);
  }
  DebugSample=z;
  acc=(signed long long)z*d->ODDEVEN;
  acc+=1ULL<<15;
  acc>>=16;
  acc+=s;
  if (acc> 0x7FFFFFFF) acc= 0x7FFFFFFF;
  if (acc<-0x7FFFFFFF) acc=-0x7FFFFFFF;
  acc=acc*d->GAIN;
  if (acc> 0x7FFFFFFF) acc= 0x7FFFFFFF;
  if (acc<-0x7FFFFFFF) acc=-0x7FFFFFFF;
  if (acc>;)/>)
  {
    acc=0x7FFFFFFFULL-acc;
    acc=acc*acc;
    acc+=1ULL<<30;
    acc>>=31;
    /*acc=acc*acc;
    acc+=1ULL<<30;
    acc>>=31;*/
    acc=0x7FFFFFFFULL-acc;
  }
  else
  {
    acc=0x7FFFFFFFULL+acc;
    acc=acc*acc;
    acc+=1ULL<<30;
    acc>>=31;
    /*acc=acc*acc;
    acc+=1ULL<<30;
    acc>>=31;*/
    acc=acc-0x7FFFFFFFULL;
  }
  if (acc>0)
  {
    acc+=1ULL<<8;
    acc>>=9;
  }
  else
  {
    acc=-acc;
    acc+=1ULL<<8;
    acc>>=9;
    acc=-acc;
  }
  if (acc>0x7FFFFFFF) acc=0x7FFFFFFF;
  if (acc<-0x7FFFFFFF) acc=-0x7FFFFFFF;
  return acc;
}

static void SetBiQuad(double Fs, double Fc, double BW, double G, BIQUAD *f)
{
    double          k,
                    om,
                    x;
    double          fi;
    double a0,a1,a2,b1,b2;
    k = pow(10, G / 40);	/* relative gain */
    BW = BW / (Fc - BW / 2);	/* bandwidth in octaves */
    om = 2 * M_PI * Fc / Fs;	/* normalized frequency in radians */
    fi = sinh(log(2) / 2 * BW * om / sin(om)) * sin(om);	/* stuff */
    x = 1 + fi / k;		/* b0 */
    a0 = (1 + fi * k) / x;
    a1 = -2 * cos(om) / x;
    a2 = (1 - fi * k) / x;
    b1 = -a1;
    b2 = -(1 - fi / k) / x;
    f->a0=65536.0*a0;
    f->a1=65536.0*a1;
    f->a2=65536.0*a2;
    f->b1=65536.0*b1;
    f->b2=65536.0*b2;
}

static void SetChebyshev(double Fs, double Fc, double ripple, int lowpass,
	       BIQUAD *f)
{
    double          x,
                    y,
                    z,
                    c,
                    v,
                    t,
                    r,
                    om,
                    m,
                    x0,
                    x1,
                    x2,
                    y1p,
                    y2,
                    k,
                    d,
                    tt,
                    tt2;
    double a0,a1,a2,b1,b2;
    // x=-cos(M_PI/4+M_PI/2);
    // y=sin(M_PI/4+M_PI/2);
    // c=-0.99915455413031497832540334286332;
    // v=0.041111761828599317357934264608497;
    c = -cos(M_PI / 4);
    v = sin(M_PI / 4);
    if (ripple > 0) {
	t = 100.0 / (100.0 - ripple);
	x = sqrt(t * t - 1);
	t = 1 / x;
	r = t + sqrt(t / x);
	y = 0.5 * log(r + 1);
	z = 0.5 * log(r - 1);
	t = exp(z);
	z = (t + 1 / t) / 2;
	t = exp(y);
	c *= (t - 1 / t) / 2 / z;
	v *= (t + 1 / t) / 2 / z;
    }
    tt = 2 * tan(0.5);
    tt2 = tt * tt;
    om = 2 * M_PI * Fc / Fs;
    m = c * c + v * v;
    d = 4 - 4 * c * tt + m * tt2;
    x0 = tt2 / d;
    x1 = x0 * 2;
    x2 = x0;
    y1p = (8 - 2 * m * tt2) / d;
    y2 = (-4 - 4 * c * tt - m * tt2) / d;
    if (lowpass)
	k = sin(0.5 - om / 2) / sin(0.5 + om / 2);
    else
	k = -cos(om / 2 + 0.5) / cos(om / 2 - 0.5);
    d = 1 + k * (y1p - y2 * k);
    a0 = (x0 - k * (x1 - x2 * k)) / d;
    // coeff[1]=(k*(-2*(x0+x2)+x1*k)+x1)/d;
    a1 = 2 * a0;
    a2 = a0;
    // coeff[2]=(k*(x0*k-x1)+x2)/d;
    b1 = (k * (2 + y1p * k - 2 * y2) + y1p) / d;
    b2 = (-k * (k + y1p) + y2) / d;
    if (!lowpass) {
	a1 = -a1;
	b1 = -b1;
    }
    f->a0=65536.0*a0;
    f->a1=65536.0*a1;
    f->a2=65536.0*a2;
    f->b1=65536.0*b1;
    f->b2=65536.0*b2;
}

static signed long BiQuad(signed long s, BIQUAD *q)
{
  signed long long acc;
  acc=(signed long long)s*q->a0;
  acc+=(signed long long)q->mem0*q->a1;
  acc+=(signed long long)q->mem1*q->a2;
  acc+=(signed long long)q->mem2*q->b1;
  acc+=(signed long long)q->mem3*q->b2;
  q->mem1=q->mem0;
  q->mem0=s;
  q->mem3=q->mem2;
  acc>>=16;
  if (acc>0x7FFFFF) acc=0x7FFFFF;
  if (acc<-0x7FFFFF) acc=-0x7FFFFF;
  return q->mem2=acc;
}

signed long stereo_delay[440];
int stereo_delay_p;

#pragma argsused
int main(int argc, char* argv[])
{
  static const char hdr1[]={
    0x52,0x49,0x46,0x46,0x00,0x00,0x00,0x00,0x57,0x41,0x56,0x45,0x66,0x6D,0x74,0x20,
    0x10,0x00,0x00,0x00,0x01,0x00,0x02,0x00,0x44,0xAC,0x00,0x00,0x10,0xB1,0x02,0x00,
    0x04,0x00,0x10,0x00,0x64,0x61,0x74,0x61,0x00,0x00,0x00,0x00
  };


  HANDLE fin;
  HANDLE fout;
  signed short *inbuf;
  signed short *outbuf;
  unsigned long total_bytes;
  unsigned long br;
  unsigned long pos=0;
  unsigned long rpos=0;
  HPF dc_offset_remove={0,0x100000000ULL*50/7022};
  HPF dc_offset_remove1={0,0x100000000ULL*10/7022};
  HPF dc_offset_remove2={0,0x100000000ULL*10/7022};

  BIQUAD preeq_b={0,0,0,0,0,0,0,0};

  SQ_DIST dist1={256UL*200,65536*65/100,{0,0x100000000ULL*20/7022}};
  //SQ_DIST dist2={256UL*70,{0,0x100000000ULL*30/7022}};

  SQ_DIST power_amp={256UL*1,65536*30/100,{0,0x100000000ULL*10/7022}};


  BIQUAD sand_remove={0,0,0,0,0,0,0,0,0};
  BIQUAD sand_remove2={0,0,0,0,0,0,0,0,0};

  BIQUAD out_eq[19];

  fout = CreateFile("out.wav",GENERIC_WRITE,0,NULL,CREATE_ALWAYS,FILE_ATTRIBUTE_NORMAL, 0);
  WriteFile(fout,hdr1,sizeof(hdr1),&br,NULL);
  fin = CreateFile("6.wav",GENERIC_READ,0,NULL,OPEN_EXISTING,FILE_ATTRIBUTE_NORMAL,0);
  SetFilePointer(fin,0x28,NULL,FILE_BEGIN);
  ReadFile(fin,&total_bytes,4,&br,NULL);
  //total_bytes=2*44100*25;
  inbuf=malloc(total_bytes);
  outbuf=malloc(total_bytes*2);
  ReadFile(fin,inbuf,total_bytes,&br,NULL);

  SetBiQuad(44100,1600,2000,14,&preeq_b);

  SetChebyshev(44100,4500,1,1,&sand_remove);
  SetChebyshev(44100,4500,1,1,&sand_remove2);

  memset(out_eq,0,sizeof(out_eq));

/*  SetBiQuad(44100,80,20*2,4,out_eq+0);
  SetBiQuad(44100,100,25*2,4,out_eq+1);
  SetBiQuad(44100,125,35*2,4,out_eq+2);
  SetBiQuad(44100,160,40*2,3,out_eq+3);
  SetBiQuad(44100,200,50*2,3,out_eq+4);
  SetBiQuad(44100,250,70*2,3,out_eq+5);
  SetBiQuad(44100,320,80*2,3,out_eq+6);
  SetBiQuad(44100,400,100*2,2,out_eq+7);
  SetBiQuad(44100,500,140*2,1,out_eq+8);
  SetBiQuad(44100,640,160*2,0,out_eq+9);
  SetBiQuad(44100,800,200*2,0,out_eq+10);
  SetBiQuad(44100,1000,250*2,0,out_eq+11);
  SetBiQuad(44100,1280,320*2,0,out_eq+12);
  SetBiQuad(44100,1600,400*2,1,out_eq+13);
  SetBiQuad(44100,2000,500*2,2,out_eq+14);
  SetBiQuad(44100,2560,640*2,3,out_eq+15);
  SetBiQuad(44100,3200,800*2,4,out_eq+16);
  SetBiQuad(44100,4000,1000*2,5,out_eq+17);
  SetBiQuad(44100,5000,1280*2,6,out_eq+18);*/

#define bw_scale (1.5)
/*  SetBiQuad(44100,80,20*bw_scale,7.3,out_eq+0);
  SetBiQuad(44100,100,25*bw_scale,7,out_eq+1);
  SetBiQuad(44100,125,35*bw_scale,6.6,out_eq+2);
  SetBiQuad(44100,160,40*bw_scale,5.9,out_eq+3);
  SetBiQuad(44100,200,50*bw_scale,4.8,out_eq+4);
  SetBiQuad(44100,250,70*bw_scale,3.3,out_eq+5);
  SetBiQuad(44100,320,80*bw_scale,2.1,out_eq+6);
  SetBiQuad(44100,400,100*bw_scale,1.5,out_eq+7);
  SetBiQuad(44100,500,125*bw_scale,0.9,out_eq+8);
  SetBiQuad(44100,640,160*bw_scale,-2.9,out_eq+9);
  SetBiQuad(44100,800,200*bw_scale,-0.9,out_eq+10);
  SetBiQuad(44100,1000,250*bw_scale,0.7,out_eq+11);
  SetBiQuad(44100,1280,320*bw_scale,2,out_eq+12);
  SetBiQuad(44100,1600,400*bw_scale,4.5,out_eq+13);
  SetBiQuad(44100,2000,500*bw_scale,4.4,out_eq+14);
  SetBiQuad(44100,2560,640*bw_scale,5.3,out_eq+15);
  SetBiQuad(44100,3200,800*bw_scale,5.2,out_eq+16);
  SetBiQuad(44100,4000,1000*bw_scale,6.2,out_eq+17);
  SetBiQuad(44100,5000,1280*bw_scale,5,out_eq+18);*/

  SetBiQuad(44100,80,20*bw_scale,8.3,out_eq+0);
  SetBiQuad(44100,100,25*bw_scale,8,out_eq+1);
  SetBiQuad(44100,125,35*bw_scale,7.6,out_eq+2);
  SetBiQuad(44100,160,40*bw_scale,6.9,out_eq+3);
  SetBiQuad(44100,200,50*bw_scale,6.8,out_eq+4);
  SetBiQuad(44100,250,70*bw_scale,5.3,out_eq+5);
  SetBiQuad(44100,320,80*bw_scale,4.1,out_eq+6);
  SetBiQuad(44100,400,100*bw_scale,3.5,out_eq+7);
  SetBiQuad(44100,500,125*bw_scale,0.9,out_eq+8);
  SetBiQuad(44100,640,160*bw_scale,-1.9,out_eq+9);
  SetBiQuad(44100,800,200*bw_scale,-1.9,out_eq+10);
  SetBiQuad(44100,1000,250*bw_scale,3.7,out_eq+11);
  SetBiQuad(44100,1280,320*bw_scale,2,out_eq+12);
  SetBiQuad(44100,1600,400*bw_scale,5.5,out_eq+13);
  SetBiQuad(44100,2000,500*bw_scale,7.4,out_eq+14);
  SetBiQuad(44100,2560,640*bw_scale,7.3,out_eq+15);
  SetBiQuad(44100,3200,800*bw_scale,7.2,out_eq+16);
  SetBiQuad(44100,4000,1000*bw_scale,7.2,out_eq+17);
  SetBiQuad(44100,5000,1280*bw_scale,8,out_eq+18);

  while(pos<(total_bytes))
  {
     int i;
     signed long sample=256UL*inbuf[rpos++];
     sample=HighPassFilter(sample,&dc_offset_remove);
     //Pre EQ
     sample=BiQuad(sample/6,&preeq_b);
     //Distortion
     sample=SquareDist(sample,&dist1)/8;
     sample=HighPassFilter(sample,&dc_offset_remove1);
     //Remove sand
     sample=BiQuad(sample,&sand_remove);
     //EQ
     for(i=0; i<sizeof(out_eq)/sizeof(BIQUAD); i++)
     {
          sample=BiQuad(sample,out_eq+i);
     }
     //Power Amp
     //DebugSample=sample;
     sample=SquareDist(sample,&power_amp);
     sample=HighPassFilter(sample,&dc_offset_remove2);
     //Remove sand
     sample=BiQuad(sample,&sand_remove2);
     //Write out
     outbuf[pos++]=sample/256;
     outbuf[pos++]=stereo_delay[stereo_delay_p]/256;
     //outbuf[pos++]=DebugSample/256;
     stereo_delay[stereo_delay_p]=sample;
     stereo_delay_p++;
     if (stereo_delay_p>=sizeof(stereo_delay)/sizeof(stereo_delay[0])) stereo_delay_p=0;
  }
  total_bytes*=2;
  WriteFile(fout,outbuf,total_bytes,&br,NULL);
  SetFilePointer(fout,0x28,NULL,FILE_BEGIN);
  WriteFile(fout,&total_bytes,4,&br,NULL);
  SetFilePointer(fout,0x04,NULL,FILE_BEGIN);
  total_bytes+=0x24;
  WriteFile(fout,&total_bytes,4,&br,NULL);
  CloseHandle(fout);
  CloseHandle(fin);
  free(inbuf);
  free(outbuf);
  return 0;
}
//---------------------------------------------------------------------------

Предэквализация тут в виде биквадратного фильтра (естественно, настраиваемого - центральная частота, добротность, затухание/усиление), затем собственно ограничитель по мотивам характеристик железного диффкаскада с регулировкой четные/нечетные, затем девятнадцатиполосный третьоктавный эквалайзер. Потом исключительно понты - еще один слабенький ограничитель, типа нелинейности выходного каскада. Еще в паре мест чебышевские фильтры НЧ для среза песка. Ах да, еще 10мс задержка во второй канал, ревер еще не прицепил сюда (но это уже мелочи).

Все считается в целых числах (кроме расчетов коэффициентов биквадратных фильтров, но это не реалтайм, а необходимо по мере кручения ручек) и требует порядка 120 операций MAC на отсчет, т.е. планируемый мной камень на тактовой частоте 48МГц будет обрабатывать отсчет примерно 6мкс, а время между отсчетами 22мкс (при дискретизации 44100 Гц). Так что вопрос о мощности проца особо не стоит.

В общем, имеем возможность забабахать первый в мире DIY цифровой гитарный процессор :good:

Делаем проект?

PS Если кто хочет попробовать покомпилить и поиграться с софтом, имейте в виду, что на вход надо кормить исключительно .wav-файлы 16 бит, 44100 Гц, моно и без всяких метаданных, ибо полноценный разбор входного файла я не делал

Наверх of the page up there ^

#73 ericdriven

Отправлено 26 мая 2010 - 15:47

проект был бы очень интересен, если будет проет то буду собирать

Наверх of the page up there ^

#74 Rst7

Отправлено 26 мая 2010 - 16:01

Проект будет. Кстати, надо бы обсудить кое-какие вопросы не плана обработки, а плана реализации управления и прочего. Щас я тему создам специально для этого, там и обсудим.

Наверх of the page up there ^

#75 stranger

Отправлено 04 июня 2010 - 08:32

Rst7 (26.5.2010, 16:01) писал:

Есть продолжение? Интересно ведь!
Кстати, если статья об эмуляции лампы на полевиках на сайте guitarjfet недоступна, то есть еще вот такая ссылка на блог гтлаб.

Наверх of the page up there ^

#76 Rst7

Отправлено 04 июня 2010 - 11:50

Цитата

Есть продолжение?

http://guitar-gear.r...?showtopic=1646

Наверх of the page up there ^

#77 Маклауд

Злой модер

Отправлено 06 октября 2010 - 12:42

Вот, кстати, собран по принципу, описанному в первой статье:

______ACDC.mp3 (671,39К)
Количество загрузок:858

Наверх of the page up there ^

#78 Кондратий

Вездесущий

Отправлено 06 октября 2010 - 15:12

Маклауд, звук хорош.
Звук похож очень, только все-таки чуток более резкий, чем лампа. А может мне кажется. : )

А принцип из первого поста нежизнеспособен?

Общаемся грамотно, уважительно и корректно.
Соблюдаем правила, не плодим бесполезных тем.
Хотим поблагодарить - нажимаем справа "+", наказать - нажимаем "-". Не стесняемся помогать модераторам кнопкой "Жалоба" слева.
Пусть наш форум будет лучше.

Наверх of the page up there ^

#79 Маклауд

Злой модер

Отправлено 06 октября 2010 - 15:20

Кондратий, ну дык в линию писал) Конечно, потом гора доработок напильником. Но честные! типо спикерсимуляция, панорама, ревер, нормализация и т.п.
Из первого поста лучше не использовать. Придумал ещё вариант, можно не только с операми использовать. Можно по идее любую диодную ограничалку на эту конструкцию сменить.

Наверх of the page up there ^

#80 Nickalaich

Отправлено 07 января 2015 - 00:28

Господа, а возможны такие апгрейды для всяких тда? Сымитировать перегруз оконечника. У меня есть пара TDA2822M для издевательств, но у неё внутри неизвестного номинала резистор в ООС.