Metamath Proof Explorer


Theorem fourierswlem

Description: The Fourier series for the square wave F converges to Y , a simpler expression for this special case. (Contributed by Glauco Siliprandi, 11-Dec-2019)

Ref Expression
Hypotheses fourierswlem.t T = 2 π
fourierswlem.f F = x if x mod T < π 1 1
fourierswlem.x X
fourierswlem.y Y = if X mod π = 0 0 F X
Assertion fourierswlem Y = if X mod T 0 π 1 1 + F X 2

Proof

Step Hyp Ref Expression
1 fourierswlem.t T = 2 π
2 fourierswlem.f F = x if x mod T < π 1 1
3 fourierswlem.x X
4 fourierswlem.y Y = if X mod π = 0 0 F X
5 simpr X mod π = 0 2 X π 2 X π
6 2z 2
7 6 a1i X mod π = 0 2 X π 2
8 pirp π +
9 mod0 X π + X mod π = 0 X π
10 3 8 9 mp2an X mod π = 0 X π
11 10 birani X mod π = 0 2 X π X π
12 divides 2 X π 2 X π k k 2 = X π
13 7 11 12 syl2anc X mod π = 0 2 X π 2 X π k k 2 = X π
14 5 13 mpbid X mod π = 0 2 X π k k 2 = X π
15 2cnd k 2
16 picn π
17 16 a1i k π
18 zcn k k
19 15 17 18 mulassd k 2 π k = 2 π k
20 17 18 mulcld k π k
21 15 20 mulcomd k 2 π k = π k 2
22 19 21 eqtrd k 2 π k = π k 2
23 22 adantr k k 2 = X π 2 π k = π k 2
24 17 18 15 mulassd k π k 2 = π k 2
25 24 adantr k k 2 = X π π k 2 = π k 2
26 id k 2 = X π k 2 = X π
27 26 eqcomd k 2 = X π X π = k 2
28 27 adantl k k 2 = X π X π = k 2
29 3 recni X
30 29 a1i k k 2 = X π X
31 16 a1i k k 2 = X π π
32 18 adantr k k 2 = X π k
33 2cnd k k 2 = X π 2
34 32 33 mulcld k k 2 = X π k 2
35 pire π
36 pipos 0 < π
37 35 36 gt0ne0ii π 0
38 37 a1i k k 2 = X π π 0
39 30 31 34 38 divmuld k k 2 = X π X π = k 2 π k 2 = X
40 28 39 mpbid k k 2 = X π π k 2 = X
41 23 25 40 3eqtrrd k k 2 = X π X = 2 π k
42 1 a1i k k 2 = X π T = 2 π
43 41 42 oveq12d k k 2 = X π X T = 2 π k 2 π
44 15 17 mulcld k 2 π
45 2ne0 2 0
46 45 a1i k 2 0
47 37 a1i k π 0
48 15 17 46 47 mulne0d k 2 π 0
49 18 44 48 divcan3d k 2 π k 2 π = k
50 49 adantr k k 2 = X π 2 π k 2 π = k
51 43 50 eqtrd k k 2 = X π X T = k
52 simpl k k 2 = X π k
53 51 52 eqeltrd k k 2 = X π X T
54 53 ex k k 2 = X π X T
55 54 a1i X mod π = 0 2 X π k k 2 = X π X T
56 55 rexlimdv X mod π = 0 2 X π k k 2 = X π X T
57 14 56 mpd X mod π = 0 2 X π X T
58 2re 2
59 58 35 remulcli 2 π
60 1 59 eqeltri T
61 2pos 0 < 2
62 58 35 61 36 mulgt0ii 0 < 2 π
63 62 1 breqtrri 0 < T
64 60 63 elrpii T +
65 mod0 X T + X mod T = 0 X T
66 3 64 65 mp2an X mod T = 0 X T
67 57 66 sylibr X mod π = 0 2 X π X mod T = 0
68 67 orcd X mod π = 0 2 X π X mod T = 0 X mod T = π
69 odd2np1 X π ¬ 2 X π k 2 k + 1 = X π
70 10 69 sylbi X mod π = 0 ¬ 2 X π k 2 k + 1 = X π
71 70 biimpa X mod π = 0 ¬ 2 X π k 2 k + 1 = X π
72 15 18 mulcld k 2 k
73 72 adantr k 2 k + 1 = X π 2 k
74 1cnd k 2 k + 1 = X π 1
75 16 a1i k 2 k + 1 = X π π
76 73 74 75 adddird k 2 k + 1 = X π 2 k + 1 π = 2 k π + 1 π
77 15 18 mulcomd k 2 k = k 2
78 77 oveq1d k 2 k π = k 2 π
79 18 15 17 mulassd k k 2 π = k 2 π
80 1 eqcomi 2 π = T
81 80 a1i k 2 π = T
82 81 oveq2d k k 2 π = k T
83 78 79 82 3eqtrd k 2 k π = k T
84 16 mullidi 1 π = π
85 84 a1i k 1 π = π
86 83 85 oveq12d k 2 k π + 1 π = k T + π
87 86 adantr k 2 k + 1 = X π 2 k π + 1 π = k T + π
88 1 44 eqeltrid k T
89 18 88 mulcld k k T
90 89 17 addcomd k k T + π = π + k T
91 90 adantr k 2 k + 1 = X π k T + π = π + k T
92 76 87 91 3eqtrrd k 2 k + 1 = X π π + k T = 2 k + 1 π
93 peano2cn 2 k 2 k + 1
94 72 93 syl k 2 k + 1
95 94 17 mulcomd k 2 k + 1 π = π 2 k + 1
96 95 adantr k 2 k + 1 = X π 2 k + 1 π = π 2 k + 1
97 id 2 k + 1 = X π 2 k + 1 = X π
98 97 eqcomd 2 k + 1 = X π X π = 2 k + 1
99 98 adantl k 2 k + 1 = X π X π = 2 k + 1
100 29 a1i k 2 k + 1 = X π X
101 94 adantr k 2 k + 1 = X π 2 k + 1
102 37 a1i k 2 k + 1 = X π π 0
103 100 75 101 102 divmuld k 2 k + 1 = X π X π = 2 k + 1 π 2 k + 1 = X
104 99 103 mpbid k 2 k + 1 = X π π 2 k + 1 = X
105 92 96 104 3eqtrrd k 2 k + 1 = X π X = π + k T
106 105 oveq1d k 2 k + 1 = X π X mod T = π + k T mod T
107 modcyc π T + k π + k T mod T = π mod T
108 35 64 107 mp3an12 k π + k T mod T = π mod T
109 108 adantr k 2 k + 1 = X π π + k T mod T = π mod T
110 35 a1i k 2 k + 1 = X π π
111 64 a1i k 2 k + 1 = X π T +
112 0re 0
113 112 35 36 ltleii 0 π
114 113 a1i k 2 k + 1 = X π 0 π
115 2timesgt π + π < 2 π
116 8 115 ax-mp π < 2 π
117 116 1 breqtrri π < T
118 117 a1i k 2 k + 1 = X π π < T
119 modid π T + 0 π π < T π mod T = π
120 110 111 114 118 119 syl22anc k 2 k + 1 = X π π mod T = π
121 106 109 120 3eqtrd k 2 k + 1 = X π X mod T = π
122 121 ex k 2 k + 1 = X π X mod T = π
123 122 a1i X mod π = 0 ¬ 2 X π k 2 k + 1 = X π X mod T = π
124 123 rexlimdv X mod π = 0 ¬ 2 X π k 2 k + 1 = X π X mod T = π
125 71 124 mpd X mod π = 0 ¬ 2 X π X mod T = π
126 125 olcd X mod π = 0 ¬ 2 X π X mod T = 0 X mod T = π
127 68 126 pm2.61dan X mod π = 0 X mod T = 0 X mod T = π
128 0xr 0 *
129 35 rexri π *
130 iocgtlb 0 * π * X mod T 0 π 0 < X mod T
131 128 129 130 mp3an12 X mod T 0 π 0 < X mod T
132 131 gt0ne0d X mod T 0 π X mod T 0
133 132 neneqd X mod T 0 π ¬ X mod T = 0
134 pm2.53 X mod T = 0 X mod T = π ¬ X mod T = 0 X mod T = π
135 134 imp X mod T = 0 X mod T = π ¬ X mod T = 0 X mod T = π
136 127 133 135 syl2anr X mod T 0 π X mod π = 0 X mod T = π
137 128 a1i X mod T = π 0 *
138 129 a1i X mod T = π π *
139 modcl X T + X mod T
140 3 64 139 mp2an X mod T
141 140 rexri X mod T *
142 141 a1i X mod T = π X mod T *
143 id X mod T = π X mod T = π
144 36 143 breqtrrid X mod T = π 0 < X mod T
145 35 eqlei2 X mod T = π X mod T π
146 137 138 142 144 145 eliocd X mod T = π X mod T 0 π
147 146 iftrued X mod T = π if X mod T 0 π 1 1 = 1
148 147 adantl X mod π = 0 X mod T = π if X mod T 0 π 1 1 = 1
149 oveq1 x = X x mod T = X mod T
150 149 breq1d x = X x mod T < π X mod T < π
151 150 ifbid x = X if x mod T < π 1 1 = if X mod T < π 1 1
152 1ex 1 V
153 negex 1 V
154 152 153 ifex if X mod T < π 1 1 V
155 151 2 154 fvmpt X F X = if X mod T < π 1 1
156 3 155 ax-mp F X = if X mod T < π 1 1
157 140 a1i X mod T < π X mod T
158 id X mod T < π X mod T < π
159 157 158 ltned X mod T < π X mod T π
160 159 necon2bi X mod T = π ¬ X mod T < π
161 160 iffalsed X mod T = π if X mod T < π 1 1 = 1
162 156 161 eqtrid X mod T = π F X = 1
163 162 adantl X mod π = 0 X mod T = π F X = 1
164 148 163 oveq12d X mod π = 0 X mod T = π if X mod T 0 π 1 1 + F X = 1 + -1
165 1pneg1e0 1 + -1 = 0
166 164 165 eqtrdi X mod π = 0 X mod T = π if X mod T 0 π 1 1 + F X = 0
167 166 oveq1d X mod π = 0 X mod T = π if X mod T 0 π 1 1 + F X 2 = 0 2
168 167 adantll X mod T 0 π X mod π = 0 X mod T = π if X mod T 0 π 1 1 + F X 2 = 0 2
169 2cn 2
170 169 45 div0i 0 2 = 0
171 170 a1i X mod T 0 π X mod π = 0 X mod T = π 0 2 = 0
172 iftrue X mod π = 0 if X mod π = 0 0 F X = 0
173 4 172 eqtr2id X mod π = 0 0 = Y
174 173 ad2antlr X mod T 0 π X mod π = 0 X mod T = π 0 = Y
175 168 171 174 3eqtrrd X mod T 0 π X mod π = 0 X mod T = π Y = if X mod T 0 π 1 1 + F X 2
176 136 175 mpdan X mod T 0 π X mod π = 0 Y = if X mod T 0 π 1 1 + F X 2
177 iftrue X mod T 0 π if X mod T 0 π 1 1 = 1
178 177 adantr X mod T 0 π ¬ X mod π = 0 if X mod T 0 π 1 1 = 1
179 140 a1i X mod T 0 π ¬ X mod π = 0 X mod T
180 35 a1i X mod T 0 π ¬ X mod π = 0 π
181 iocleub 0 * π * X mod T 0 π X mod T π
182 128 129 181 mp3an12 X mod T 0 π X mod T π
183 182 adantr X mod T 0 π ¬ X mod π = 0 X mod T π
184 ax-1cn 1
185 184 16 mulcomi 1 π = π 1
186 84 185 eqtr3i π = π 1
187 186 oveq1i π + π 2 X T = π 1 + π 2 X T
188 169 16 mulcomi 2 π = π 2
189 1 188 eqtri T = π 2
190 189 oveq1i T X T = π 2 X T
191 112 63 gtneii T 0
192 3 60 191 redivcli X T
193 flcl X T X T
194 192 193 ax-mp X T
195 zcn X T X T
196 194 195 ax-mp X T
197 16 169 196 mulassi π 2 X T = π 2 X T
198 190 197 eqtri T X T = π 2 X T
199 198 oveq2i π + T X T = π + π 2 X T
200 169 196 mulcli 2 X T
201 16 184 200 adddii π 1 + 2 X T = π 1 + π 2 X T
202 187 199 201 3eqtr4ri π 1 + 2 X T = π + T X T
203 202 a1i π = X mod T π 1 + 2 X T = π + T X T
204 id π = X mod T π = X mod T
205 modval X T + X mod T = X T X T
206 3 64 205 mp2an X mod T = X T X T
207 204 206 eqtrdi π = X mod T π = X T X T
208 207 oveq1d π = X mod T π + T X T = X - T X T + T X T
209 29 a1i π = X mod T X
210 60 recni T
211 210 196 mulcli T X T
212 211 a1i π = X mod T T X T
213 209 212 npcand π = X mod T X - T X T + T X T = X
214 203 208 213 3eqtrrd π = X mod T X = π 1 + 2 X T
215 214 oveq1d π = X mod T X π = π 1 + 2 X T π
216 184 200 addcli 1 + 2 X T
217 216 16 37 divcan3i π 1 + 2 X T π = 1 + 2 X T
218 215 217 eqtrdi π = X mod T X π = 1 + 2 X T
219 1z 1
220 zmulcl 2 X T 2 X T
221 6 194 220 mp2an 2 X T
222 zaddcl 1 2 X T 1 + 2 X T
223 219 221 222 mp2an 1 + 2 X T
224 223 a1i π = X mod T 1 + 2 X T
225 218 224 eqeltrd π = X mod T X π
226 225 10 sylibr π = X mod T X mod π = 0
227 226 necon3bi ¬ X mod π = 0 π X mod T
228 227 adantl X mod T 0 π ¬ X mod π = 0 π X mod T
229 179 180 183 228 leneltd X mod T 0 π ¬ X mod π = 0 X mod T < π
230 iftrue X mod T < π if X mod T < π 1 1 = 1
231 156 230 eqtrid X mod T < π F X = 1
232 229 231 syl X mod T 0 π ¬ X mod π = 0 F X = 1
233 178 232 oveq12d X mod T 0 π ¬ X mod π = 0 if X mod T 0 π 1 1 + F X = 1 + 1
234 233 oveq1d X mod T 0 π ¬ X mod π = 0 if X mod T 0 π 1 1 + F X 2 = 1 + 1 2
235 1p1e2 1 + 1 = 2
236 235 oveq1i 1 + 1 2 = 2 2
237 2div2e1 2 2 = 1
238 236 237 eqtr2i 1 = 1 + 1 2
239 232 238 eqtr2di X mod T 0 π ¬ X mod π = 0 1 + 1 2 = F X
240 iffalse ¬ X mod π = 0 if X mod π = 0 0 F X = F X
241 4 240 eqtr2id ¬ X mod π = 0 F X = Y
242 241 adantl X mod T 0 π ¬ X mod π = 0 F X = Y
243 234 239 242 3eqtrrd X mod T 0 π ¬ X mod π = 0 Y = if X mod T 0 π 1 1 + F X 2
244 176 243 pm2.61dan X mod T 0 π Y = if X mod T 0 π 1 1 + F X 2
245 132 necon2bi X mod T = 0 ¬ X mod T 0 π
246 245 iffalsed X mod T = 0 if X mod T 0 π 1 1 = 1
247 id X mod T = 0 X mod T = 0
248 247 36 eqbrtrdi X mod T = 0 X mod T < π
249 248 iftrued X mod T = 0 if X mod T < π 1 1 = 1
250 156 249 eqtrid X mod T = 0 F X = 1
251 246 250 oveq12d X mod T = 0 if X mod T 0 π 1 1 + F X = - 1 + 1
252 251 oveq1d X mod T = 0 if X mod T 0 π 1 1 + F X 2 = - 1 + 1 2
253 neg1cn 1
254 184 253 165 addcomli - 1 + 1 = 0
255 254 oveq1i - 1 + 1 2 = 0 2
256 255 170 eqtri - 1 + 1 2 = 0
257 256 a1i X mod T = 0 - 1 + 1 2 = 0
258 1 oveq2i X T = X 2 π
259 2cnne0 2 2 0
260 16 37 pm3.2i π π 0
261 divdiv1 X 2 2 0 π π 0 X 2 π = X 2 π
262 29 259 260 261 mp3an X 2 π = X 2 π
263 29 169 16 45 37 divdiv32i X 2 π = X π 2
264 258 262 263 3eqtr2i X T = X π 2
265 264 oveq2i 2 X T = 2 X π 2
266 29 16 37 divcli X π
267 266 169 45 divcan2i 2 X π 2 = X π
268 265 267 eqtr2i X π = 2 X T
269 6 a1i X T 2
270 id X T X T
271 269 270 zmulcld X T 2 X T
272 268 271 eqeltrid X T X π
273 66 272 sylbi X mod T = 0 X π
274 273 10 sylibr X mod T = 0 X mod π = 0
275 274 iftrued X mod T = 0 if X mod π = 0 0 F X = 0
276 4 275 eqtr2id X mod T = 0 0 = Y
277 252 257 276 3eqtrrd X mod T = 0 Y = if X mod T 0 π 1 1 + F X 2
278 277 adantl ¬ X mod T 0 π X mod T = 0 Y = if X mod T 0 π 1 1 + F X 2
279 129 a1i ¬ X mod T 0 π ¬ X mod T = 0 π *
280 60 rexri T *
281 280 a1i ¬ X mod T 0 π ¬ X mod T = 0 T *
282 140 a1i ¬ X mod T 0 π ¬ X mod T = 0 X mod T
283 pm4.56 ¬ X mod T 0 π ¬ X mod T = 0 ¬ X mod T 0 π X mod T = 0
284 283 biimpi ¬ X mod T 0 π ¬ X mod T = 0 ¬ X mod T 0 π X mod T = 0
285 olc X mod T = 0 X mod T 0 π X mod T = 0
286 285 adantl X mod T π X mod T = 0 X mod T 0 π X mod T = 0
287 128 a1i X mod T π ¬ X mod T = 0 0 *
288 129 a1i X mod T π ¬ X mod T = 0 π *
289 141 a1i X mod T π ¬ X mod T = 0 X mod T *
290 0red ¬ X mod T = 0 0
291 140 a1i ¬ X mod T = 0 X mod T
292 modge0 X T + 0 X mod T
293 3 64 292 mp2an 0 X mod T
294 293 a1i ¬ X mod T = 0 0 X mod T
295 neqne ¬ X mod T = 0 X mod T 0
296 290 291 294 295 leneltd ¬ X mod T = 0 0 < X mod T
297 296 adantl X mod T π ¬ X mod T = 0 0 < X mod T
298 simpl X mod T π ¬ X mod T = 0 X mod T π
299 287 288 289 297 298 eliocd X mod T π ¬ X mod T = 0 X mod T 0 π
300 299 orcd X mod T π ¬ X mod T = 0 X mod T 0 π X mod T = 0
301 286 300 pm2.61dan X mod T π X mod T 0 π X mod T = 0
302 284 301 nsyl ¬ X mod T 0 π ¬ X mod T = 0 ¬ X mod T π
303 35 a1i ¬ X mod T 0 π ¬ X mod T = 0 π
304 303 282 ltnled ¬ X mod T 0 π ¬ X mod T = 0 π < X mod T ¬ X mod T π
305 302 304 mpbird ¬ X mod T 0 π ¬ X mod T = 0 π < X mod T
306 modlt X T + X mod T < T
307 3 64 306 mp2an X mod T < T
308 307 a1i ¬ X mod T 0 π ¬ X mod T = 0 X mod T < T
309 279 281 282 305 308 eliood ¬ X mod T 0 π ¬ X mod T = 0 X mod T π T
310 128 a1i X mod T π T 0 *
311 35 a1i X mod T π T π
312 141 a1i X mod T π T X mod T *
313 ioogtlb π * T * X mod T π T π < X mod T
314 129 280 313 mp3an12 X mod T π T π < X mod T
315 310 311 312 314 gtnelioc X mod T π T ¬ X mod T 0 π
316 315 iffalsed X mod T π T if X mod T 0 π 1 1 = 1
317 140 a1i X mod T π T X mod T
318 311 317 314 ltnsymd X mod T π T ¬ X mod T < π
319 318 iffalsed X mod T π T if X mod T < π 1 1 = 1
320 156 319 eqtrid X mod T π T F X = 1
321 316 320 oveq12d X mod T π T if X mod T 0 π 1 1 + F X = - 1 + -1
322 321 oveq1d X mod T π T if X mod T 0 π 1 1 + F X 2 = - 1 + -1 2
323 df-2 2 = 1 + 1
324 323 negeqi 2 = 1 + 1
325 184 184 negdii 1 + 1 = - 1 + -1
326 324 325 eqtr2i - 1 + -1 = 2
327 326 oveq1i - 1 + -1 2 = 2 2
328 divneg 2 2 2 0 2 2 = 2 2
329 169 169 45 328 mp3an 2 2 = 2 2
330 237 negeqi 2 2 = 1
331 327 329 330 3eqtr2i - 1 + -1 2 = 1
332 331 a1i X mod T π T - 1 + -1 2 = 1
333 4 a1i X mod T π T Y = if X mod π = 0 0 F X
334 311 317 ltnled X mod T π T π < X mod T ¬ X mod T π
335 314 334 mpbid X mod T π T ¬ X mod T π
336 247 113 eqbrtrdi X mod T = 0 X mod T π
337 336 adantl X mod π = 0 X mod T = 0 X mod T π
338 127 orcanai X mod π = 0 ¬ X mod T = 0 X mod T = π
339 338 145 syl X mod π = 0 ¬ X mod T = 0 X mod T π
340 337 339 pm2.61dan X mod π = 0 X mod T π
341 335 340 nsyl X mod T π T ¬ X mod π = 0
342 341 iffalsed X mod T π T if X mod π = 0 0 F X = F X
343 333 342 320 3eqtrrd X mod T π T 1 = Y
344 322 332 343 3eqtrrd X mod T π T Y = if X mod T 0 π 1 1 + F X 2
345 309 344 syl ¬ X mod T 0 π ¬ X mod T = 0 Y = if X mod T 0 π 1 1 + F X 2
346 278 345 pm2.61dan ¬ X mod T 0 π Y = if X mod T 0 π 1 1 + F X 2
347 244 346 pm2.61i Y = if X mod T 0 π 1 1 + F X 2