X
تبلیغات
رایتل

مهندسی عمران-زلزله(محمدجواد خسرویانی)

بدست آوردن پاسخ فرکانسی تابع همراه با کد متلب Frequency response function

Systems respond differently to inputs of different frequencies. Some systems may amplify components of certain frequencies, and attenuate components of other frequencies. The way that the system output is related to the system input for different frequencies is called the frequency response of the system.

The frequency response is the relationship between the system input and output in the Fourier Domain.

Fourier Transfer Block.svg

In this system, X(jω) is the system input, Y(jω) is the system output, and H(jω) is the frequency response. We can define the relationship between these functions as:

Y(j\omega) = H(j\omega)X(j\omega)


\frac{Y(j\omega)}{X(j\omega)} = H(j\omega)



at the follow the codes for Matlab would show the FRF in an example


clear all
A1=20; A2=30; f1=5; f2=15; wn1=2*pi*f1; wn2=2*pi*f2;
zeta1=0.02; zeta2=0.01;
wd1=sqrt(1-zeta1^2)*wn1; wd2=sqrt(1-zeta2^2)*wn2;
fs=100; T1=10; t1=[0:1/fs:T1-1/fs]; h=(A1/wd1)*exp(-zeta1*wn1*t1).*sin(wd1*t1)+(A2/wd2)*exp(-zeta2*wn2*t1).*sin(wd2*t1);

T= 2000; % T=10000;

randn('state',0);
x=randn(1,T*fs);
y=filter(h,1,x)/fs; % scaled appropriately
Tr=4; N=Tr*fs; % Tr=20;

[Sxx, ~]=cpsd(x,x, hanning(N),N/2, N, fs, 'twosided');
[Syy, f]=cpsd(y,y, hanning(N),N/2, N, fs, 'twosided');

H=fft(h)/fs; % scaled appropriately.

f1=fs*(0:length(H)-1)/length(H);

figure (1)
plot(f,10*log10(fs*Sxx), f, zeros(size(f)), 'r:')

xlabel('Frequency (Hz)')

ylabel('Estimate of \itS x x\rm(\itf\rm) (dB)')

axis([0 30 -10 10])

figure(2)

plot(f,10*log10(Syy),f1,10*log10(abs(H).^2/fs), 'r:')

xlabel('Frequency (Hz)')

ylabel('Estimate of \itS y y\rm(\itf\rm) (dB)')

axis([0 30 -100 -20])

figure(3)

plot(f,10*log10(abs(Sxy)),f1,10*log10(abs(H)/fs), 'r:')

xlabel('Frequency (Hz)')

ylabel('Estimate of |\itS x y\rm(\itf\rm)| (dB)')

axis([0 30 -60 -20])

figure(4)

plot(f,unwrap(angle(Sxy)),f1,unwrap(angle(H)), 'r:')

xlabel('Frequency (Hz)')

ylabel('Estimate of arg\itS x y\rm(\itf\rm))(rad)')

axis([0 30 -3.5 0]) 

http://s6.picofile.com/file/8190471700/12222.jpg



http://s6.picofile.com/file/8190471734/212112.jpg

تاریخ ارسال: سه‌شنبه 5 خرداد 1394 ساعت 20:19 | نویسنده: M.khosraviani | چاپ مطلب
نظرات (3)
سه‌شنبه 11 مهر 1396 12:20
امیر
امتیاز: 0 0
لینک نظر
با سلام ...
اگر بخاهیم با متلب برنامه ای بنویسیم که مثل برنامه سایزمو سیگنال محتوای فرکانسی رو ترسیم کنه ،،یعنی محتوای فرکانسی در مقابل دامنه فوریه ،،آیا می تونید راهنمائی کنید؟
پاسخ:
salam
be post zir dar blogam berid kamelan codi neveshte shode ke hamin kar ro mikone
http://mkhosraviani.blogsky.com/1395/04/25/post-79/1
سه‌شنبه 3 مرداد 1396 14:08
نوید
امتیاز: 0 0
لینک نظر
سلام وقت بخیر
آقای دکتر میخواستم برای یک سازه سه طبقه فولادی تابع پاسخ فرکانس بکشم و اشکال مودی و فرکانس های سازه رو بکشم ممنون میشم بگین کجای کد رو باید تغییر بدهم
پاسخ:
salam
tebghe tozihat Y9w) khoroji sazas ke shoma masalan pasokhe freuency tabaghe 3 ro midin va X(W) ham vorodiye sazas ke darvaghe hamon tahrik ya bri hast ke be saze vared mikonid
shoma bayad dar code jahye h ke bar hast bare khodeton ro bezarid hamchenin be mesal tavajoh konid motavajeh mishid ke in mesal omade bare randomi tolid karde va pasokhesh ro gerefte va FRF ro rasm karde
سه‌شنبه 12 بهمن 1395 17:57
abedin
امتیاز: 0 0
لینک نظر
با سلام خدمت جناب دکتر
من میخواستم برای یه خرپا دو بعدی تابع پاسخ فرکانسیشو بدست بیارم، ممنون میشم منو راهنمایی کنید.
(البته من شکل مود ها و فرکانس های خرپامو بدست آوردم)
پاسخ:
salam berahati bar asase tozihate blog mitonid mohasebe konid

tabe frequency ke H bashe be insorate
FFT( pasokh) taghsim bar meghdare bar
ke agar dar code blog ham bebinid be hamin sorate
H=fft(h)/fs
برای نمایش آواتار خود در این وبلاگ در سایت Gravatar.com ثبت نام کنید. (راهنما)
نام :
پست الکترونیک :
وب/وبلاگ :
ایمیل شما بعد از ثبت نمایش داده نخواهد شد