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

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

بدست آوردن پاسخ فرکانسی تابع همراه با کد متلب 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

نظرات 6 + ارسال نظر
مارال یکشنبه 2 دی 1397 ساعت 23:58

با عرض سلام وادب
با تشکر از سایت بسیار خوبتون بنده برنامه ای با متلب برای سازه ای سه طبقه برای بدست آواردن میزان جابجایی سازه و فرکانس ها و شکل مودی سازه نیاز دارم.
خیلی از لطفتون ممنون میشم اگر به بنده در این رابطه کمک نمایید

salam
linke zir moraje konid

http://mkhosraviani.blogsky.com/1392/10/03/post-32/

آرش سه‌شنبه 13 آذر 1397 ساعت 23:40

سلام
شما Sxy رو در اسکریپت قید نکردید و برنامه ایراد میگیره

gheyd shode

کوروش چهارشنبه 27 دی 1396 ساعت 02:17

با سلام و احترام
میخواستم برای طبقه اول و دهم سازه تابع پاسخ فرکانس رسم کنم و برای 2 حالت بدون میراگر جرمی و با tmd مقایسه کنم (قله های کوچکتر نمودار در حالت با tmd نسبت به بدون tmd نشانگر تاثیر بهینه میراگر در کاهش جابه جایی در اثر زلزله) ممنون میشم راهنمایی کنید کدوم قسمت کد باید تغییر بدم
با سپاس

salam

khodeton code ro motavajeh beshid jaye input va output haton avaz mishe ke X input hast hamon niro va Y output haast ke pasokhe shomast shetab masalan

امیر سه‌شنبه 11 مهر 1396 ساعت 12:20

با سلام ...
اگر بخاهیم با متلب برنامه ای بنویسیم که مثل برنامه سایزمو سیگنال محتوای فرکانسی رو ترسیم کنه ،،یعنی محتوای فرکانسی در مقابل دامنه فوریه ،،آیا می تونید راهنمائی کنید؟

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

سلام وقت بخیر
آقای دکتر میخواستم برای یک سازه سه طبقه فولادی تابع پاسخ فرکانس بکشم و اشکال مودی و فرکانس های سازه رو بکشم ممنون میشم بگین کجای کد رو باید تغییر بدهم

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

abedin سه‌شنبه 12 بهمن 1395 ساعت 17:57

با سلام خدمت جناب دکتر
من میخواستم برای یه خرپا دو بعدی تابع پاسخ فرکانسیشو بدست بیارم، ممنون میشم منو راهنمایی کنید.
(البته من شکل مود ها و فرکانس های خرپامو بدست آوردم)

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 ثبت نام کنید. (راهنما)
ایمیل شما بعد از ثبت نمایش داده نخواهد شد