[python] 파이썬 적분 샘플 코드 (scipy.integrate.quad)

파이썬에서 적분(integral) 방법의 간단한 소개입니다. 아래와 같은 식을 사용자 정의 함수로 정의하고 적분값을 구해 보겠습니다.


아래는 코드 샘플입니다.


import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import quad


x = np.linspace(0, 2, 16)

def f(x): return x**2
y = np.zeros(x.size)
for i in range(x.size):
    y[i], error = quad(f, 0, x[i])
    
plt.plot(x, y)


실행결과는 아래와 같습니다.


끝.

댓글

댓글 쓰기

이 블로그의 인기 게시물

[PLC] PLC 아날로그 입출력 기본

이미지
 PLC 아날로그 입출력 개념 및 결선  아날로그 신호가 PLC로 들어오면 우선 A/D 컨버터(Analog-Digital Converter)를 거치게 됩니다. A/D 컨버터는 아날로그 형태의 신호를 PLC에서 처리 가능하도록 디지털 형태로 변환시켜 줍니다.  PLC에서 다루는 아날로그 신호는 아래와 같은 형태가 있습니다. - 전압(Voltage) - 전류(Current) - 저항(Resistance) 신호의 형태에 따라 신호가 취급되는 값의 범위가 있습니다. 가장 흔히 사용되는 형태는 전류 형태이며 4-20mA 값을 주로 사용합니다. A/D 컨버터가 아날로그 값을 디지털로 변환할 때 사용되는 bit수를 분해능이라 합니다. 일반적으로는 16bit 길이가 사용됩니다. 다시말해 4-20mA 사이의 16mA를 16bit로 쪼개어 표현한다는 의미입니다. 따라서 분해능에 사용된 bit수가 많을수록 더 세밀한 단위를 취급할 수 있습니다. 아래 그림은 4-20mA 출력을 가지는 온도 트랜스미터(temperature transmitter)의 예입니다. 온도 트랜스미터와 온도센서가 연결되어 있습니다. 온도 정보는 트랜스미터에서 섭씨 0~100도로 계량화(calibrated) 됩니다. 즉, 섭씨 0도는 4mA, 섭씨 100도는 20mA가 되는 것입니다. 위 그림에서 온도 정보는 PLC로 바로 들어오지 않고 트랜스미터를 거쳤습니다. 온도 정보를 PLC가 바로 해석하기 보다는 트랜스미터에서 해석하여 전압이나 전류 형태로 변환하여 PLC의 전압/전류 아날로그 카드로 연결됩니다. 만약 트랜스미터가 없다면, 온도센서의 정보를 해석하기 위해 PLC에는 온도입력 모듈이 별도로 필요하게 됩니다. 전압/전류 아날로그 입력카드를 사용하면, 온도센서가 아니더라도 다양한 형태의 아날로그 데이터를 트랜스미터 전압/전류 형태로 취급할 수 있으므로, 좀 더 유현한 구조가 됩니다. 아날로그 입력 모듈의 결선법은 입력 신호가 전압인지, 전류인지에 따라 다릅니다. 전압/전류은 저항과 연...
Read more »

전력(kW) 계산하기 (직류, 교류 단상, 교류 삼상)

전력(kW) 계산하기 (직류, 교류 단상, 교류 삼상) 전류(A)와 전압(V)를 알고 있는 경우, 전압원이 종류에 따라 전력(kW)을 어떻게 계산하는지 보겠습니다. DC amps to kilowatts calculation (직류) 직류는 단순히 전류와 전압을 곱하고 1000으로 나누면 kW 단위의 전력을 구할 수 있습니다. P (kW) = I (A) × V (V) / 1000 AC single phase amps to kilowatts calculation (교류 단상) 교류 단상의 경우 역률(Power Factor)를 곱한 뒤, 전류와 RMS 전압을 곱하고, 1000으로 나누면 kW 단위의 전력을 구할 수 있습니다. P (kW) = PF × I (A) × V (V) / 1000 AC three phase amps to kilowatts calculation (교류 삼상) Calculation with line to line voltage (상간 전압 이용시) 교류 삼상의 상간 전압을 이용하여 전력을 구하는 경우, √3과 역률(Power Factor)을 곱하고, 위상 전류(phase current)와 두 활선 사이의 RMS 전압을 곱합니다. 그리고 1000으로 나누어 kW 단위의 전력을 구합니다. P (kW) = √3 × PF × I (A) × V L-L(V) / 1000 Calculation with line to neutral voltage (활선-중립선 이용시) 교류 삼상의 활선-중립선 전압을 이용하여 전력을 구하는 경우, 3과 역률(Power Factor)을 곱하고, 위상 전류(phase current)와  두 활선 사이의 RMS 전압을 곱합니다. 그리고 1000으로 나누어 kW 단위의 전력을 구합니다. P (kW) = 3 × PF × I (A) × V L-N(V) / 1000 끝.
Read more »

공압 속도 제어: 미터인 vs 미터아웃

이미지
Understanding Flow Controls: Meter-In vs Meter-Out 에어실린더에 작용하는 힘은 압력조절기(pressure regulator)에 의해 결정됩니다. 그리고 피스톤이 움직이는 속도는 속도제어(flow control) 밸브를 통해 이루어집니다. 실린더의 각 포트마다 속도제어 밸브가 설치되어 인출과 인입 속도를 제어하게 됩니다. 아래는 속도제어 밸브의 구조입니다. 위 그림에 대한 기호는 아래와 같이 표기됩니다. 미터아웃(Meter-Out) 방식의 실린더 속도제어 아래 그림은 미터아웃 방식에 대한 도식입니다. 미터인(Meter-In) 방식의 실린더 속도제어 아래 그림은 미터인 방식에 대한 도식입니다. 끝.
Read more »

제너 다이오드에 저항을 연결하는 이유

이미지
Why do we connect a resistor before a Zener diode? 아래 회로를 통해 저항의 크기를 변화시키면 어떠한 패턴이 발생하는지 살펴보겠습니다. 아래 그림에서 3개의 병렬회로에는 각각 12V가 인가됩니다. 그리고 동일한 제너 다이오드가 연결되어 있고 각각의 제너 다이오드 앞단에 60R, 6R, 0.6R 저항이 연결되어 있습니다. 그림에서 알수 있듯이 저항의 값을 줄이면 제너 다이오드에 걸리는 전류값이 증가합니다. 제너 다이오드에 흐르는 전류가 증가하면 덩달아 제너 다이오드에서 소비전력이 높아질 것이고, 과열로 인하여 다이오드가 손상될 것입니다. 따라서 제너 다이오드 앞단에 연결된 저항은 다이오드에 과전류가 흐르는 것을 막기 위해 사용합니다. 아래 데이터시트를 통해 좀더 구체적인 예를 들어 보겠습니다. 첫번째 항목인 P D 는 손실전력을 의미합니다. 위 시트상 P D =500mW입니다. 제너 전압 V Z 가 결정되어 있다면 아래 조건을 만족하도록 저항을 선택해야 합니다. 예를 들어 제너 전압이 V Z =5.1V 인 경우, 아래와 같이 최대 전류를 구할 수 있습니다. 이제 98mA를 초과하지 않도록 적절한 저항을 선택하여 제너 다이오드 앞단에 연결해 주면 됩니다. 끝.
Read more »

회로 차단기 용량 선정하는 방법

이미지
Sizing the Circuit Breakers 회로 차단기의 용량을 선정하는 방법 몇가지를 소개합니다. Using the Specification Sticker 회로 차단기에 연결할 제품에 제원 스티커가 있고 소비전류가 표기되어 있다면, 가장 쉽게 계산할 수 있습니다. 예를 들어, 스티커에 0.1 A @ 115 VAC, 60 Hz로 표기된 제품이 115 VAC용 회로에서 12개 사용되고 하나의 회로 차단기에 연결되는 경우라면 아래와 같이 계산할 수 있습니다. 0.1 A x 12 = 1.2 A Using the Watts/Volts Method (유효전력) 제품의 소비전력(W)값을 전압(V)로 나누어 소비전류를 구하는 방법이 있습니다. 예를 들어, 정격전압 120 VAC에서 소비전력 240 W인 제품의 소비전류는 아래와 같이 계산될 수 있습니다. 240 W / 120 V = 2 A Considering the Power Factor (역률) 위에서 사용된, 소비전력을 전압으로 나누는 방법은 역률이 1에 가까운 경우에만 유효합니다.역률이 낮은 경우 실제 소비전력은 위 방법으로 구한 값보다 높은 전류값을 가집니다. 그러므로, 역률이 0.9 이하인 제품이 회로 차단기에 연결되는 경우 역률를 고려하여 계산해야 합니다. 예를 들어, 정격전압 120 VAC에서 소비전력 240 W인 제품의 역률이 0.7이라면 아래와 같이 계산될 수 있습니다. 2 A / 0.7 = 2.8571 A 이 값은 역률를 고려하지 않는 경우에 비해 약 50% 높은 값입니다. Using the Volt Amps Method (피상전력) 만약 제품의 제원 스티커에 소비전력이 VA 단위로 표기되어 있다면 이 값을 전압(V)로 나누어 소비전류를 구하면 됩니다. 예를 들어, 소비전력이 360 VA로 표기된 120 V 제품의 경우 아래와 같이 계산하면 됩니다. 360 VA / 120 V = 3 A 이 경우에는 역률을 고려하지 않아도 됩니다. Selecting the Circuit Breaker Nationa...
Read more »

PLC 출력 형태

PLC Output Units PLC 출력 유닛은 보통 3가지 형태입니다: - 릴레이(Relay) - 트랜지스터(Transistor) - 트라이액(Triac) PLC의 출력부에 부하(Load)를 연결하기 전에 부하의 사양을 확인해야 합니다. 부하의 최대 소비전류값이 PLC 출력부의 허용전류 이내인지 살펴봅니다. 릴레이 출력(Relay Outputs) 릴레이 출력은 가장 흔하게 사용됩니다. 릴레이는 정해진 극성이 없고 AC와 DC 모두 사용할 수 있습니다. 트랜지스터 출력(Transistor Outputs) 트랜지스터 출력은 DC전류만 취급할 수 있습니다. PLC는 트랜지스터의 베이스(BASE)에 약한 전류를 흘리므로써 회로가 폐로(Closed)되어 출력 장치로 전류가 흐르게 됩니다. 트랜지스터 출력은 릴레이 출력보다 상대적으로 작은 부하만 사용할 수 있습니다. 만약 사용하고자 하는 부하가 트랜지스터 출력부의 허용범위를 초과한다면 외부 릴레이를 사용하여 간접적으로 부하를 구동시키면 됩니다. 트랜지스터 출력 타입에 따라 NPN, PNP 구분이 있습니다. 또 다른 구분으로는 BJT와 MOSFET이 있습니다. BJT 타입은 MOSFET 타입보다 작은 스위칭 용량을 가지지만 MOSFET보다 빠른 스위칭이 가능합니다. 트랜지스터 출력은 적은 전류를 취급하지만 릴레이보다 속도가 빠르고 수명이 깁니다. 트라이액 출력(Triac Output) 트라이액 출력은 AC 부하만 사용할 수 있습니다. 릴레이 출력보다 속도가 빠르고 수명이 긴 장점이 있습니다. 유도성 부하(Inductive loads)는 동작할 때 역전류(back current)를 가지는 경향이 있습니다. 역전류는 릴레이 출력부에 영향을 미칠 수 있으므로, 다이오드(diode), 베리스터(varistor) 등의 보호회로를 사용하여 PLC 출력부를 보호합니다.
Read more »

[PLC] 채터링 현상과 입력 필터

이미지
전자 회로 내의 스위치나 릴레이의 접점이 붙거나 떨어질 때 기계적인 진동에 의해 실제로는 매우 짧은 시간 안에 접점이 붙었다가 떨어지는 것을 반복하는 현상입니다. PLC 측면에서 본다면 이는 특정 입력단자의 전기신호가 짧은 시간 안에 ON/OFF가 수 많이 반복되는 형태로 인식됩니다. 디지털 입력 필터란 채터링에 의한 노이즈를 제거하기 위한 소프트웨어적 방법이며, 대부분의 PLC 프로그래밍 소프트웨어에서 설정이 가능합니다. 예를 들어 입력 필터 시간을 3ms로 설정한다면, 위 그림에서 최소한 3ms 동안 신호의 ON 상태를 유지해야 PLC에서 해당 입력을 ON으로 인정하다는 의미입니다. 반대로 최소한 3ms 동안 신호의 OFF 상태를 유지해야 PLC에서 해당 입력을 OFF로 인정하게 됩니다. LS ELECTRIC PLC의 경우 XG5000 소프트웨어를 통해 기본 파라미터 항목에서 "표준 입력필터" 항목을 통해 설정할 수 있으며, 통상 3ms가 사용됩니다. 끝.
Read more »

NPN, PNP 트랜지스터 차이점

이미지
Difference Between an NPN and a PNP Transistor 위 그림에서 기호는 아래 의미입니다. B: BASE C: COLLECTOR E: EMITTER NPN, PNP 둘다 BJT(bipolar junction transistor: 접합형  트랜지스터)입니다. 전류를 이용하여 제어가 되며, 전류를 증폭하는 기능이 있습니다. base에 전류가 흐르면, emitter와 collector 사이에는 더 많은 전류가 흐르게 됩니다. 즉, base는 emitter-collector 간 전류를 흐르게 하는 스위치 역할과 더불어 전류 증폭기 역할을 하게 됩니다. 이는 NPN, PNP 구분에 상관없이 동일한 특성입니다. NPN, PNP 트랜지스터의 차이는, 트랜지스터 핀에 전원이 인가되는 방법의 차이에 기인합니다. NPN 트랜지스터는 collector와 base에 양전압을 연결합니다.  NPN 트랜지스터의 base에 양전압이 인가되면 emitter 방향으로 전류가 공급(sourcing, 소스)가 되면서 collector-emitter간 연결이 만들어지고, 전류도 흐르게 됩니다. 다시 정리하자면, base에 임계전압인 0.7V이상의 전압이 인가되면 전류는 base에서 emitter로 흐르게 되고, 덩달아 collector-emitter간에 흐르는 전류도 점차 증가하게 되는 형태입니다. PNP 트랜지스터는 emitter에 양전압을 연결하고, base에 음전압(또는, emitter에 인가된 전압보다 낮은 전압)을 연결합니다. PNP 트랜지스터는 base에 음전압이 인가되면 base에서 ground 방향으로 전류가 빠지면서(sink, 싱크) emitter-collector간 연결이 만들어지고 전류가 흐르게 됩니다. NPN, PNP 트랜지스터는 emitter-collector간 전류의 방향이 다르기 때문에, 위 그림에서 확인할 수 있듯이 부하(load)가 연결되는 위치 또한 다르게 됩니다. 끝.
Read more »

[전기 기초] 저항의 정격전력(Watt) 표기의 의미

이미지
저항을 보면 1/2W, 1/4W 등이 표기되어 있습니다. 이게 무엇을 의미하는지 조사한 내용을 공유합니다.  1/2W, 1/4W 등의 표시는 저항이 최대로 버틸 수 있는 정격전력(power rating)을 의미합니다. 저항을 통해 전류가 흐를 때 저항에는 열이 발생됩니다. 만약 전류량이 커서 많은 열이 발생하게 되면 전자소자들이 녹기 시작하고 쇼트를 발생시킬 수도 있습니다. 우선 전력 P에 관한 다음 식들을 살펴보겠습니다. \(P=VI\) \(P=\frac{V^{2}}{R}\) \(p=I^{2}R\) 예를 들어, 800Ω 저항이 있고 12V에서 동작하는 LED를 구동시키는 회로가 있다고 가정해 보겠습니다. 이때 전력은 아래와 같이 구할 수 있습니다. \(P=\frac{V^{2}}{R}=\frac{(12^{2})V}{800Ω}=0.18W\) 위 회로에서는 1/4W 저항을 선정하면 충분하겠습니다. 다른 예로, 150Ω 저항에 15V로 동작하는 모터회로를 가정해 봅시다. \(P=\frac{V^{2}}{R}=\frac{(15^{2})V}{150Ω}=1.5W\) 위 모터회로에는 2W의 저항을 선정하면 될 것 같습니다. 주로 사용되는 정격전력에는 1/4W(0.25W), 1/2W(0.5W), 1W, 2W, 5W, 25W, ... 등이 있습니다. 끝.
Read more »

[python] 파이썬 pyplot 2차원 그래프 샘플 코드

이미지
예제 코드를 실행하기 위한 환경입니다. 라이브러리: matplotlib, numpy IDE : Spyder 데이터를 포함하는 2개의 배열로부터 2차원 그래프를 생성하겠습니다. 아래는 샘플 데이터입니다. 파이썬 코드는 아래와 같습니다. import matplotlib.pyplot as plt import numpy as np m = np.array([1.0, 2.0, 4.0, 6.0, 9.0, 11.0]) V = np.array([0.13, 0.26, 0.50, 0.77, 1.15, 1.36]) plt.plot(m, V, 'o') plt.xlabel('m (kg)') plt.ylabel('V (l)') plt.show() 실행하게 되면 아래와 같은 그래프를 출력할 수 있습니다. 만약 아래와 같이 함수를 포함하는 데이터를 이용하여 그래프를 생성하는 경우를 다루어 보겠습니다. 함수의 입력값에 대한 수열과 각 입력값에 대한 출력값 수열을 생성하고 그래프를 생성하면 됩니다. 아래는 샘플 코드입니다. import matplotlib.pyplot as plt import numpy as np x0 = 0.0 x1 = 10.0 dx = 0.1 n = 100 # x축 데이터용 n x 1 행렬을 만들고 0으로 초기화 x = np.zeros((n, 1), float) # y축 데이터용 n x 1 행렬을 만들고 0으로 초기화 y = np.zeros((n, 1), float) for i in range(n): x[i] = x0 + i*dx y[i] = np.sin(x[i]) plt.plot(x, y) plt.show() 아래는 실행 결과입니다. 위 코드는 linspace 함수를 사용하면 보다 간결하게 작성할 수 있습니다. import matplotlib.pyplot as plt import numpy as np x = np.linspace(0, 10, 100) y = np.sin(x) plt...
Read more »