ลองเริ่มต้นง่ายๆนะครับ

สมมติว่ามีฟังก์ชัน y=f(x)=3+5x-x^2 และต้องการหาค่าสูงสุด

เราจะเขียน code หาค่าสูงสุดได้อย่างไร

ผมแนะนำดังนี้ครับ

- ต้องเริ่มจาก ตรวจสอบว่าเรามีพื้นฐานการเขียนโปรแกรมมาแล้วอย่างไรบ้าง

- หากไม่มีพื้นฐานเลย อันนี้ก็ต้องเริ่มต้นนับหนึ่งกันเลยทีดียว

ที่สำคัญ เราต้องเริ่มต้นเองก่อน หากมีปัญหา น่าจะมีคนช่วยดู ช่วยปรับ code ให้ได้

MathGuy
ร่วมแบ่งปันความรู้และความเห็นแล้ว 1988 ครั้ง - ได้รับดาวแล้ว 323 ดวง - โหวตเพิ่มดาว

ดาวน์โหลดหนังสือ SCILAB ฟรี ได้ที่
http://home.npru.ac.th/~t3058/webscilab/index.php

gte723s
ร่วมแบ่งปันความรู้และความเห็นแล้ว 18 ครั้ง - ได้รับดาวแล้ว 150 ดวง - โหวตเพิ่มดาว

เพิ่งเล่น scilab วันนี้เอง พอดีเห็นว่ามันเหมือนกับ matlab มากๆ เลยมาทำการบ้านน้า MathGuy ดู

ปรากฏว่าเขียนฟังก์ชันแล้ว มันไม่รู้จะเซฟไว้ที่ใหนให้ เพราะผมเขียนใน scipad เลย เขียนในรูป input แทน

ลองดูนะครับว่าพอใช้ได้มั้ย

N=input('Enter three numbers you want to find GCD in form [n1 n2 n3]:')
n1=N(1);n2=N(2);n3=N(3);
i=1;
while(modulo(n2,(n1/i))~=0)
  i=i+1;
end
GCD1=n1/i;

j=1;
while(modulo(n3,(GCD1/j))~=0)
 j=j+1;
end
GCD=GCD1/j

ขอบคุณอาจารย์ปิยะ มากครับที่ให้ Link ที่มีประโยชน์มาก พร้อมทั้งให้ ดาวน์โหลดหนังสือ SCILAB ของอาจารย์ได้ฟรีอีกด้วย

ผมเห็นมีอบรมของอาจารย์สำหรับแนะนำการใช้โปรแกรม SCILAB ช่วงวันที่ 23-24 สิงหาคมนี้ ขออนุญาตประชาสัมพันธ์ตรงนี้ด้วยนะครับ

...............................................................................


แบบฟอร์มตอบรับเข้าร่วมการอบรมเชิงปฏิบัติการ
“การแก้ไขโจทย์คณิตศาสตร์และฟิสิกส์ด้วยโปรแกรม SCILAB  ครั้งที่ ๓/๒๕๕๑”
วันที่ ๒๓ – ๒๔  สิงหาคม ๒๕๕๑
ณ อาคารวิศวกรรมและเทคโนโลยี ห้อง EB ๖๐๘
ค่าลงทะเบียนท่านละ ๑๘๐ บาท

ชื่อ ……………………………………… นามสกุล ………………………………………………
ตำแหน่ง ………………………………… สังกัดหน่วยงาน ……………………………………...
โทรศัพท์ ……………………….. โทรสาร ……………………….. มือถือ ……………………...
E-mail: …………………………………………………………………………………………….

ลงชื่อ …………………………………….
( ………………………………………. )
*********************************************

หมายเหตุ
ผู้เข้าอบรมจะได้รับเอกสารและซีดีประกอบการฝึกอบรม ท่านละ ๑ ชุด
พร้อมทั้งอาหารกลางวัน ๒ มื้อ และอาหารว่าง ๔ มื้อ
ผู้เข้าอบรมสามารถชำระค่าลงทะเบียนได้ ณ จุดลงทะเบียน หน้าห้องอบรม EB ๕๐๘
กรุณาส่งใบตอบรับมาที่เบอร์โทรสาร ๐๓๔ – ๒๖๑ – ๐๖๕  ก่อนวันที่จะเข้ารับการอบรม ๗ วัน
ตรวจสอบรายชื่อผู้เข้ารับการอบรม หรือสอบถามข้อมูลเพิ่มเติมได้ที่ ผศ.ดร.ปิยะ โควินท์ทวีวัฒน์  โทรศัพท์ ๐๘๙ – ๔๕๖ – ๕๐๕๐

กำหนดการอบรม
เรื่อง การแก้ไขโจทย์คณิตศาสตร์และฟิสิกส์ด้วยโปรแกรม SCILAB  ครั้งที่ ๓/๒๕๕๑
ณ อาคารวิศวกรรมและเทคโนโลยี ห้อง EB ๖๐๘

วันที่ ๒๓ สิงหาคม ๒๕๕๑
 
เวลา หัวข้อของการอบรม   
08:00 – 08:20 น. ลงทะเบียน   
08:20 – 08:30 น. พิธีเปิดการอบรม   
08:30 – 10:30 น. แนะนำโปรแกรมภาษา SCILAB
ประเภทของข้อมูล การสร้างตัวแปรสเกล่าร์ เวกเตอร์ และเมทริกซ์
การแก้โจทย์สมการหลายตัวแปร (ระบบสมการเชิงเส้น)
การแก้โจทย์สมการโพลิโนเมียล   
10:30 – 10:45 น. อาหารว่าง   
10:45 – 12:00 น. การใช้งานโอเปอเรเตอร์ทางคณิตศาสตร์
–0 โอเปอเรเตอร์ Arithmetic
–1 โอเปอเรเตอร์ Relational & Logical   
12:00 – 13:00 น. พักรับประทานอาหารกลางวัน   
13:00 – 14:30 น. การใช้งานฟังก์ชันพื้นฐานทางคณิตศาสตร์
–2 ฟัง–3 ก์ชันพื้นฐานทาง–4 คณิตศาสตร์ทั่วไป
–5 ฟัง–6 ก์ชันตรีโกณมิติ
–7 ฟัง–8 ก์ชันพื้นฐานทาง–9 สถิติ
การแก้โจทย์ทางด้านตรีโกณมิติและทางด้านสถิติ   
14:30 – 14:45 น. อาหารว่าง   
14:45 – 16:30 น. ตัวอย่างการประยุกต์ใช้งาน
–10 การหาตัวคูณร่วมน้อย (ค.ร.น.) และตัวหารร่วมมาก (ห.ร.ม.)
–11 การแยกตัวประกอบของ–12 สมการพหุนาม
–13 การหาปริพันธ์จำกัดเขต
–14 การแก้สมการอนุพันธ์อันดับหนึ่ง–15  

วันที่ ๒๔ สิงหาคม ๒๕๕๑
 
เวลา หัวข้อของการอบรม   
08:30 – 10:30 น. การเขียนโปรแกรมด้วย SCILAB
–16 การใช้คำสั่ง–17 การวนลูป
–18 การใช้คำสั่ง–19 ทดสอบเง–20 ื่อนไข   
10:30 – 10:45 น. อาหารว่าง   
10:45 – 12:00 น. การพัฒนาฟังก์ชันย่อยขึ้นมาใช้งานเฉพาะทาง   
12:00 – 13:00 น. พักรับประทานอาหารกลางวัน   
13:00 – 14:30 น. พื้นฐานการวาดกราฟ
–21 การวาดกราฟสอง–22 มิติ
–23 การวาดกราฟสามมิติ   
14:30 – 14:45 น. อาหารว่าง   
14:45 – 16:30 น. การแก้โจทย์สมการการเคลื่อนที่
การแก้สมการภาคตัดกรวยแบบต่างๆ   
16:30 – 16:45 น. พิธีปิดการอบรม 

MathGuy
ร่วมแบ่งปันความรู้และความเห็นแล้ว 1988 ครั้ง - ได้รับดาวแล้ว 323 ดวง - โหวตเพิ่มดาว

ขอบคุณอาจารย์ปิยะค่ะ ดาวน์โหลดหนังสือแล้ว จะหาเวลาศึกษาสักตั้ง

ครูไผ่
ร่วมแบ่งปันความรู้และความเห็นแล้ว 2312 ครั้ง - ได้รับดาวแล้ว 216 ดวง - โหวตเพิ่มดาว

จะเขียนโปรแกรมที่หาโดเมนของปัญหา ที่ถูกจำกัดด้วยสมการเชิงเส้น ตั้งแต่ 2 เส้นขึ้นไปค่ะ เช่น max F(x,y)=x-4y
                     max f(x,y)=x+y
          โดยมีเงื่อนไข
                -x -y   <= -3
            -3x + 2y  >= -4
             -2x + y   <= 0
              2x + y   <= 12
ไม่รู้ว่าจะต้องเขียนโปรแกรมอย่างไรดีนะค่ะ ถามหาค่าเริ่มต้นของ  x ทุกที อ.ช่วยหน่อยนะค่ะ

การแก้ปัญหา linear programming
มีวิธีการแก้ที่สำคัญคือ Simplex method
หากเข้าใจขั้นตอนวิธีการนี้ ก็สามารถเขียนโปแกรมตามนั้นได้เลย


..................................................................................

การหาโดเมน หรือ feasible region (บริเวณที่สอดคล้องเงื่อนไข)

ในกรณี 2 ตัวแปร ให้วาดรูปในระนาบ XY ดูก่อน

MathGuy
ร่วมแบ่งปันความรู้และความเห็นแล้ว 1988 ครั้ง - ได้รับดาวแล้ว 323 ดวง - โหวตเพิ่มดาว

พอดีเรียนอยู่ปี 3 นะคะ ได้ทำโปรเจคเรื่องpoints orientation & scale in image อยากทราบรายละเอียดเกี่ยวกับเรื่องนี้นะคะพอช่วยได้มั้ยคะ แล้วจะใช้โปรแกรมอะไรเขียนดีคะ ขอบคุณล่วงหน้าคะที่ตอบปัญหาให้คะ

ต้องรู้ว่า ปัญหาของเราคืออะไร
มีเงื่อนไข มีอะไรกำหนดมาให้บ้าง แล้วต้องการทำให้ได้อะไร
มีขั้นตอน มีการคำนวณ การแสดงผลอะไรเกี่ยวข้องบ้าง

เมื่อเข้าใจปัญหา ก็จะพอตอบได้ว่า ควรจะใช้เครื่องมือ หรือโปรแกรมอะไร

สำหรับปัญหาการคำนวณ การแสดงผลทั่วๆไป Scilab ก็น่าจะใช้ได้

MathGuy
ร่วมแบ่งปันความรู้และความเห็นแล้ว 1988 ครั้ง - ได้รับดาวแล้ว 323 ดวง - โหวตเพิ่มดาว

โจทย์มีอยู่ว่าc=88 ,m=52 ,t=1ถึง12sec ,dt=1sec,g=9.81,ให้หาvที่tต่างๆ จากสูตร
v=(g*m/c*(1-exp(-(c/m*t)))); %exat soln
 v(t+dt) = ((g-(c*v(t))/m)*dt)+v(t)%numerical
ได้ลองเขียนเป็นmatlabมาแล้วค่ะ คือ
g=9.81; %ค่าจี
m=52; %มวล
c=88; %ค่าซี
t=0:0.1:12; %เวลา
vex=(g*m/c*(1-exp(-(c/m*t)))); %exat soln
%plot(t,v)

%numerical
v = zeros(1,12)
dt = 0.1;
for i = 2:13,
    v(i) = (g-c*(v(i-1))/m)*dt+(v(i-1));
end

plot(t,v,'r')
hold on
plot(t,vex)
อยากทราบว่าถ้าจะเขียนเป็นscilabจะต้องเขียนยังไงบ้างค่ะ

อีกข้อคือจากโจทย์ข้างบนให้หาความเร็วของนักกระโดดร่มก่อนถึงพื้นดินที่1km
ที่t=0,v=0 เขียนเป็นscilabยังไงค่ะ

และอยากทราบว่าmatlabและscilabต่างกันยังไง
ช่วยหน่อยนะคะ ขอบคุณค่ะ

Matlab เป็นโปรแกรมที่เราต้องซื้อ

ส่วน Scilab เป็น Freeware ใช้ได้ฟรี และมีขนาดเล็กกว่ามาก

การคำนวณทั่วไปๆ สามารถใช้ได้ดีพอๆกัน คำสั่งต่างๆ ก็เหมือนกัน
(อาจจะมีต่างกันบ้าง ที่คำสั่งการ plot กราฟบางคำสั่ง)

code บางที copy ของ Matlab เอาไปใช้ได้เลย 

ทดลองดูครับ ดีที่สุด

(เข้าไปที่ website ของ Scilab แล้ว Dowload โปรแกรม)



ลักษณะ code จะเหมือนภาษา C แต่ถูกปรับให้ง่ายและสะดวกกว่า

MathGuy
ร่วมแบ่งปันความรู้และความเห็นแล้ว 1988 ครั้ง - ได้รับดาวแล้ว 323 ดวง - โหวตเพิ่มดาว

ใน คหพต. 193  มีอยู่ว่า  "...อยากทราบว่าถ้าจะเขียนเป็นscilabจะต้องเขียนยังไงบ้างค่ะ"

คำตอบใน คหพต.194  เป็น  "Matlab เป็นโปรแกรมที่เราต้องซื้อ
ส่วน Scilab เป็น Freeware ใช้ได้ฟรี และมีขนาดเล็กกว่ามาก....ฯลฯ..."



น่าจะตอบสั้นๆแค่นี้
(เข้าไปที่ website ของ Scilab แล้ว Dowload โปรแกรม)

NpEducate
ร่วมแบ่งปันความรู้และความเห็นแล้ว 632 ครั้ง - ได้รับดาวแล้ว 186 ดวง - โหวตเพิ่มดาว

ขอตอบคุณ sunshine จากความคิดเห็น 193 เพิ่มเติม

(หวังว่าคงไม่ช้าเกินไปนะครับ)
.....................................................................................

ที่ต้องปรับเปลี่ยน มีดังนี้
(1) การเขียน comment จาก % ใน Matlab ให้แก้เป็น // สำหรับ Scilab (เหมือนภาษา C)

(2) การจะสั่ง plot( x, y) ทั้ง x และ y ต้องมีขนาดเท่ากัน

แก้ไขได้ด้วยการใช้ linspace( , , )

t=linspace(0,12,120); 

จะมีค่า t จำนวน 120 ค่า เริ่มจาก 0 ถึง 12  (เพิ่มทีละ 0.1 ตรงกับค่า dt)

ซึ่งจะมีขนาดตรงกับ

vex=zeros(1,120);

และ

v = zeros(1,120);

การ plot ใช้คำสั่ง subplot( , , ) ช่วย

.........................................................................................

โค๊ดที่ปรับและแก้ไข ให้ใช้กับ scilab ได้ เป็นดังนี้ครับ

g=9.81; //ค่าจี
m=52; //มวล
c=88; //ค่าซี

t=linspace(0,12,120);
vex=zeros(1,120);

for i=1:120, //เวลา
  vex(i)=(g*m/c*(1-exp(-(c/m*t(i))))); //exat soln
end

//numerical
v = zeros(1,120);
dt = 0.1;
for i = 2:120,
    v(i) = (g-c*(v(i-1))/m)*dt+(v(i-1));
end

subplot(1,2,1)
plot(t,v,'r')

subplot(1,2,2)
plot(t,vex,'b')

MathGuy
ร่วมแบ่งปันความรู้และความเห็นแล้ว 1988 ครั้ง - ได้รับดาวแล้ว 323 ดวง - โหวตเพิ่มดาว

ได้ผลการรัน ดังนี้ครับ

MathGuy
ร่วมแบ่งปันความรู้และความเห็นแล้ว 1988 ครั้ง - ได้รับดาวแล้ว 323 ดวง - โหวตเพิ่มดาว