C語言 例題5-1利用 Euler尤拉 修正近似法 解一階常微分方程ODE y' = -y+x+1 0<= x <= 1 , y(0)=1 , h=0.2
/***********************************************
****EULER METHOD FOR DIFFERENTIAL EQUATIONS*****
***********************************************/
#include<stdio.h>
#include<math.h>
/*Define the RHS of the first order differential equation here(Ex: dy/dx=f(x,y)) */
double f(double x, double y){
return -y+x+1;
}
double predict(double x, double y, double h)
{
// value of next y(predicted) is returned
double y1p = y + h * f(x, y);
return y1p;
}
// using Modified Euler method
double correct(double x, double y, double x1, double y1,double h)
{
// (x, y) are of previous step
// and x1 is the increased x for next step
// and y1 is predicted y for next step
double e = 0.00001;
double y1c = y1;
do {
y1 = y1c;
y1c = y + 0.5 * h * (f(x, y) + f(x1, y1));
} while (fabs(y1c - y1) > e);
// every iteration is correcting the value
// of y using average slope
return y1c;
}
int main()
{
// here x and y are the initial
// given condition, so x=0 and y=15
double x = 0, y = 1;
// final value of x for which y is needed
double xn = 1;
// step size
double h = 0.2;
int i=1;
while (x < xn) {
double x1 = x + h;
double y1p = predict(x, y, h);
double y1c = correct(x, y, x1, y1p, h);
printf("Euler Modified method----The value of y , at x = %2.2lf is :%2.4lf\n" ,x1 ,y1c);
x = x1;
y = y1c;
i++;
}
return 0;
}
輸出畫面
Euler Modified method----The value of y , at x = 0.20 is :1.0182
Euler Modified method----The value of y , at x = 0.40 is :1.0694
Euler Modified method----The value of y , at x = 0.60 is :1.1477
Euler Modified method----The value of y , at x = 0.80 is :1.2481
Euler Modified method----The value of y , at x = 1.00 is :1.3666
2019年5月12日 星期日
2016年12月18日 星期日
2016年1月9日 星期六
2015年11月18日 星期三
2015年11月14日 星期六
實作練習作業4 (電路板1)
請下載MyDesign CH4 FIG4-12電路圖檔案
下載點 MyDesign CH4 FIG4-12 圖中有3個元件外型圖 錯誤
並參考以下網頁
http://alex9ufoexploer.blogspot.tw/2015/10/protel99se.html
從第7步驟 開始 (步驟1-步驟6 已經處理好了--畫電路圖步驟)
7) 新增一個電路板PCB (File ->New)
到
完成後 mail alex9ufo@gmail.com
下載點 MyDesign CH4 FIG4-12 圖中有3個元件外型圖 錯誤
並參考以下網頁
http://alex9ufoexploer.blogspot.tw/2015/10/protel99se.html
從第7步驟 開始 (步驟1-步驟6 已經處理好了--畫電路圖步驟)
7) 新增一個電路板PCB (File ->New)
到
16) 完成的電路板如下 列印輸出 或是產生鑽孔檔
=============================
利用 View -> Board in 3D 檢視電路板
完成後 mail alex9ufo@gmail.com
2015年10月31日 星期六
電腦繪圖
週
次
|
起 訖 日 期
|
預定進度
|
||||
章節課次
|
作業
次數
|
日常
考查
|
起訖頁數
|
|||
自08 月24日
至08 月28日
|
教學預備週
|
|||||
1
|
自08 月31日
至09 月04日
|
8/31(一)註冊開學、9/1(二)正式上課
|
||||
2
|
自09 月07日
至09 月11日
|
CH1 Protel系統軟體介紹
|
1-1~1-18
|
|||
3
|
自09 月14日
至09 月18日
|
CH2 Protel軟體的檔案說明
CH3 快速畫一個簡單電路圖
|
2-1~2-24
3-1~3-34
|
|||
4
|
自09 月21日
至09 月25日
|
CH4編輯功能的詳細說明
|
1
|
1
|
4-1~4-25
|
|
5
|
自09 月28日
至10月02日
|
9/28(一)中秋節補假
|
||||
6
|
自10月05日
至10 月09日
|
CH5畫一個具有匯流排的電路圖
10/6(二)~10/8(四)第一次段考
|
5-1-~5-16
|
|||
7
|
自10 月12日
至10 月16日
|
檢討段考試卷
CH6畫電路圖的重要功能
|
6-1~6-21
|
|||
8
|
自10 月19日
至10 月23日
|
CH7畫一個階層電路
|
7-1~7-19
|
|||
9
|
自10 月26日
至10 月30日
|
CH8 編輯元件庫的電路圖元件
|
8-1~8-24
|
|||
10
|
自11月02日
至11 月06日
|
CH9 PCB設計項目說明
|
9-1~9-24
|
|||
11
|
自11 月09日
至11 月13日
|
CH10 PCB電路板的設計流程
|
2
|
2
|
10-1~10-26
|
|
12
|
自11 月16日
至11 月20日
|
CH11 PCB佈局層說明
|
11-1~11-26
|
|||
13
|
自11 月23日
至11 月27日
|
CH12 放置和佈線功能說明
11/25(三)~11/27(五)第二次段考
|
12-1~12-12
|
|||
14
|
自11 月30日
至12月04日
|
CH12 放置和佈線功能說明
檢討段考試卷
|
12-12~12-28
|
|||
15
|
自12 月07日
至12 月11日
|
CH13 元件外形圖的詳細介紹
|
13-1~13-19
|
|||
16
|
自12 月14日
至12 月18日
|
CH14 輸出列印和設計規則說明
|
14-1~14-36
|
|||
17
|
自12 月21日
至12 月25日
|
CH15 畫一個時脈產生器的印刷電路板
|
15-1~15-25
|
|||
18
|
自12 月28日
至01月01日
|
1/1(五)開國紀念日
CH16 畫一個介面卡電路
|
3
|
3
|
16-1~16-38
|
|
19
|
自01月04日
至01 月08日
|
CH17 建立元件外形圖的元件庫
|
17-1~17-21
|
|||
20
|
自01 月11日
至01 月15日
|
期末考總複習
1/14(四)、1/15(五)期末考
|
||||
21
|
自01 月18日
至01 月22日
|
1/18(一)~1/20(三)期末考
|
||||
2015年10月27日 星期二
新建電路圖元件 實驗8-1 8-2
新建電路圖元件
在什麼情況下?需要畫電路圖元件,有兩種情況需要畫電路圖:
1. 找不到適用的元件,需要畫一個新的電路圖元件。
2. 找到類似的元件,可以修改這個電路圖元件。
要新建電路圖元件時,你需要準備哪些資料?說明如下:
1. 只要畫電路圖元件,必須知道元件接腳的資料,所以在編輯元件之前,
要準備好下列表格的內容,如下表所示(以電容元件為例):
接腳名稱(Name) 接腳編號(Number) 電氣特性(Electrical Type) 接腳設定
1 1 Passive
2 2 Passive
接腳設定欄位可以設定下列資料,如下:
(1) Dot Symbol:接腳是反相接腳?
(2) Clk Symbol:接腳是時脈接腳?
(3) Hidden:接腳是隱藏接腳?
上面的接腳資料要填入Pin對話盒中
2. 如果要完成PCB電路板,則必須知道電路圖元件的元件外形圖(Footprint),這個資料可以填入Component Text Fields對話盒中 。
常用的元件都應該找得到電路圖元件,所以要畫的電路圖元件通常是比較
少見的元件,例如:特殊IC、新設計的IC、感測器元件...等,以元件組
成單元分別說明,如下:
1. 元件圖形:只是一個元件符號,所以可以隨便畫,只要能代表這個元
件,最簡單的元件圖形可以是一個矩形。
2. 元件接腳:在畫電路圖元件之前,一定要先準備好接腳資料的表格內
容,如上表所示。
3. 元件特性值:輸入這個元件的特性值。
接腳的圓點必須在外側,不可以和矩形相接,如果要旋轉接腳,只要按SPACE鍵,就可以逆時針旋轉90度。
放好接腳後,要再移動或旋轉接腳,只要游標在接腳上,按住Mouse左鍵,可以移動接腳,按SPACE鍵,可以旋轉接腳(此時仍要按住Mouse左鍵),再放開Mouse左鍵。
File -> New Schematic library Documents
完成後 並且使用新建電路圖元件的電路圖 下載點
在什麼情況下?需要畫電路圖元件,有兩種情況需要畫電路圖:
1. 找不到適用的元件,需要畫一個新的電路圖元件。
2. 找到類似的元件,可以修改這個電路圖元件。
要新建電路圖元件時,你需要準備哪些資料?說明如下:
1. 只要畫電路圖元件,必須知道元件接腳的資料,所以在編輯元件之前,
要準備好下列表格的內容,如下表所示(以電容元件為例):
接腳名稱(Name) 接腳編號(Number) 電氣特性(Electrical Type) 接腳設定
1 1 Passive
2 2 Passive
接腳設定欄位可以設定下列資料,如下:
(1) Dot Symbol:接腳是反相接腳?
(2) Clk Symbol:接腳是時脈接腳?
(3) Hidden:接腳是隱藏接腳?
上面的接腳資料要填入Pin對話盒中
2. 如果要完成PCB電路板,則必須知道電路圖元件的元件外形圖(Footprint),這個資料可以填入Component Text Fields對話盒中 。
常用的元件都應該找得到電路圖元件,所以要畫的電路圖元件通常是比較
少見的元件,例如:特殊IC、新設計的IC、感測器元件...等,以元件組
成單元分別說明,如下:
1. 元件圖形:只是一個元件符號,所以可以隨便畫,只要能代表這個元
件,最簡單的元件圖形可以是一個矩形。
2. 元件接腳:在畫電路圖元件之前,一定要先準備好接腳資料的表格內
容,如上表所示。
3. 元件特性值:輸入這個元件的特性值。
![]() |
| 新增說明文字 |
***** 畫圖過程中 需留意*******
放好接腳後,要再移動或旋轉接腳,只要游標在接腳上,按住Mouse左鍵,可以移動接腳,按SPACE鍵,可以旋轉接腳(此時仍要按住Mouse左鍵),再放開Mouse左鍵。
File -> New Schematic library Documents
SHOW HIDDEN PINS
完成後 並且使用新建電路圖元件的電路圖 下載點
訂閱:
文章 (Atom)
MQTT 協定與 Modbus 通訊的遠端監控與控制系統架構
MQTT 協定與 Modbus 通訊的遠端監控與控制系統架構 這張圖片展示了一個 結合 MQTT 協定與 Modbus 通訊的遠端監控與控制系統架構 (主要透過 Node-RED 進行資料整合)。 系統包含三個核心部分,其運作功能說明如下: 1. ESP32 終端設備(硬體控制層...
-
數位IC設計入門-Verilog combinational logic 8 to 1 Multiplexer 多工器 Behavioral Modeling (& Test Bench) //數位IC設計入門-Verilog combinationa...
-
Line 發報機 Python TKinter (CONFIG_LINE Message API_2.py) import serial import serial.tools.list_ports import tkinter as tk from tkinter import...
-
python pip 不是内部或外部命令 -- 解決方法 要安裝 Pyqt5 1. 首先,開啟命令提示字元。 2. 輸入 pip3 install pyqt5 好像不能執行 ! ! 錯誤顯示 : ‘ pip3 ’ 不是內部或外部命令、可執行的程式或批...













