The Energy Branch, that I work in, was revealed via a short video article that displayed on a major Israeli news site called NRG (Ma’ariv).

http://www.nrg.co.il/online/1/ART2/244/584.html?hp=1&cat=459

See you there

During my work on one of my main projects, I found myself required to verifie voltage measurements and  to compare it with temperature measurements of an appliance.

Over my examination  I’ve used the Agilent 34980A with temperature multiplexer. I took samples over certain period of time and preserved the information on laptop in form of ascii file. Then I redo the same process, only this time I measured voltage instead of temperature.

Now I have two ascii files that contains my sampled temperature/voltage. Using matlab or Scilab I want to convert my voltage sampled information to temperature in order to confirm that they are approximately the same.

I have 8 scanned ports, but I’ll examine every port in a separate process.

Finding mathematical behavior for voltage – temperature slope

According to the literature, we know that the slope of thermocouple response curve changes over temperature. Using the fact that my thermocouples are from type T and examined between 0 – 80 Celsius degrees the slope acts much linear. Actually, Ihave checked it myself by giong to NIST appropriate page and plotting (both Matlab and Scilab) the curve in order to see that it’s actually linear.

Using “Basic Fitting” Matlab tool, I found both the slope and the equation of this linear behavior. We will get the same result by:

polyfit(x1,cels,1)

In Scilab it was bit difficult, but thanks to WolfData’s site I managed to receive the same result using ‘regres

s’ function.

So now I know for sure that I’m in the linear region, allowing me to write a simple equation that will convert my volatge samples to a temperature value, dumping the option of typing the full table of NIST.

In the next part I will examine (using the linear equation) my voltage and temperature samples  match suitability.

NOTE: You can enlarge any of the pictures by right clicking on it and ‘open image in a new tab’ option.

After thought I have decided to start writing my own blog in English, in order to participate others with my professional knowledge.

]]> http://vaxincity.wordpress.com/2011/04/30/moving-from-hebrew-to-english/feed/ 0 vaxincity http://vaxincity.wordpress.com/2011/04/18/mat_e_bug/ http://vaxincity.wordpress.com/2011/04/18/mat_e_bug/#comments Mon, 18 Apr 2011 07:25:45 +0000 vaxincity http://vaxincity.wordpress.com/?p=535 ]]> מספר פעמים קרה לי שביצעתי פעולה חשבונית כלשהי במטלאב וקיבלתי תוצאה מוצמדת בחזקה של e, למרות שלא רציתי בכך. החלטתי לעלות על הכתב את הפתרון הפשוט בכדי לא לשכוח אותו להבא.

בכדי להיפטר מההצמדה של ה e עלינו לבקש:

format long g

לאחר מכן לקרוא שוב לפעולה שחישבנו ולקבל את התוצאה במספר עשרוני מלא.

חג שמח לכולם

שבוע שעבר, כחלק מבדיקה כלשהי שבוצעה במקום עבודתי, נתבקשתי לאגור את הנתונים המתקבלים מאתר הניסוי ולבצע עליהם מספר מניפולציות מתמטיות כלשהן כאשר בסוף התהליך יבדקו מדדים כלשהם.

המידע נשלח מהשטח על ידי רכיבים דוגמים כלשהם, נשלח על גבי אמצעי תקשורת כלשהו למחשב מסויים המפענח את המידע (מאות בינארי לאות דצימלי) אשר נרשם בקובץ טקסט ומופרד על ידי טאבים. כמובן שכל הפרטים עד כה אינם מעניינים לצורך הרשומה ועל כן לא אפרט עליהם.

כעת, לאחר שברשותנו מספר קבצים שכאלה בעלי מאות ואלפי שורות של מספרים (זמן הדגימה ומספר פרמטרים שנמדדו באותה יחידת זמן בדידה) אנו מעוניינים לבצע עליהם מספר מניפולציות מתמטיות.

המהלכים המתמטיים התחלקו לשניים: מהלכים שהם בעצם נוסחא (כלומר לא תלוי תוכנה ספציפית) ומהלכים של קבלת המידע לתוכנה המתמטית ועיבודו הראשוני. המהלך הראשוני לא מעניין אותנו מכיוון שיכול להתבצע על כל תוכנה מתמטית, בחלק השני לא הייתי בטוח האם מטלאב ואוקטב יתפקדו בצורה זהה וזה בדיוק מה שעניין אותי לבדוק.

כמובן שהיה ניתן לבצע את כל המהלכים על גבי אקסל, אבל חיפשתי פתרון מכני או אוטומטי לכל קבוצה של דגימות וזאת מתוך נקודת הנחה שאין לי מטלאב על גבי המחשב בעבודה. כלומר אם אני מצליח לבצע את הפעולות על גבי אוקטב, אני חופשי להתשמש באוקטב בעבודה ולחסוך לעצמי זמן יקר.

שלב ראשון: קריאת הנתונים

הנתונים מגיעים בקובץ טקסט או csv. קריאת הקובץ לתוך מטלאב ואוקטב בהתאמה היא:

textread(‘c:\Documents and Settings\1.txt’  , ‘ %f %f’)

load(“-ascii”, “c:\Documents and Settings\1.txt”)

נקבל על המסך רישום של המטריצה שנקראה מתוך הקובץ.

שלב שני: מציאת הערך הקרוב ביותר במטריצה שלנו לערך נקוב שהוחלט עליו מראש

כעת אנחנו רוצים לקחת את כל שלל הנתונים שקיבלנו (מתעלמים מזמן הדגימה, רק תוצאת הדגימה בזמן נתון מעניינת אותנו). ישנו ערך כלשהו שקבענו ואנו רוצים לדעת מהי הדגימה שקיבלה את הערך הקרוב ביותר לערך שקבענו מראש.

לשם כך ביצעתי במטלאב את הפעולות:

x=מטריצת הדגימות שהכנסנו מתוך קובץ

y=2   (נניח שהערך שקבענו הנו 2)

z=abs(x-y)  (מחסרים הדגימות מהערך הנקוב, מפעילים ערך מוחלט מכיוון שרק המרחק היחסי מעניין אותנו)

min(z)   (מחפשים את הערך הקטן ביותר, שהוא בעצם הדגימה הקרובה ביותר לערך שקבענו)

לשמחתי הקוד המקביל באוקטב זהה לחלוטין לקוד של מטלאב (התמונות הן משלבי הנסיון שלי ולכן “היגרלתי” מטריצה בתמונות הללו):

שארית המהלכים היו מתמטיים לחלוטין ולכן לא טרחתי לעלותם על הכתב. הנקודה החשובה היא שבמקום העבודה לא היה לי על גבי המחשב מטלאב לביצוע פעולה מתמטית קבועה כל שהיא וכעת אוקטב פתרה את הבעיה ותפסה את מקומה.

אני חייב לציין שמנסיוני לא תמיד כך הדבר וישנם מצבים שפונקציות מסויימות אינן קיימות באוקטב. יש לבדוק כל עניין לגופו.

המשך סופ”ש נעים.

וקסמן

הכנות מקדימות

1. אני יוצא מתוך נקודת הנחה שאנחנו מכירים את פקודת wget והורדנו כבר ממקור אמין את קובץ ה tar.bz2 שמכיל בתוכו את הפריסה שלנו.

2. תחת משתמש על נוריד את הקובץ לספריה

cd /usr/src

rm linux

tar xvfj <kernel-source-file>.tar.bz2

כעת לאחר שמחקנו את קישור ה linux הקודם ופרסנו את קובץ הקרנל נוכל להתחיל בעבודה ממש.

הידור הקרנל

נגדיר קישור חדש למשתנה ה linux שלנו

ln -s <linux-releas> linux

נגדיר הרשאות, וגדיר שינויים רצויים בקרנל (menuconfig) ונתחיל להדר

chown -R root:root <linux-release>

make menuconfig

make all

make module_install

make install

יצירת תמונת initrd

update-initramfs -c -k ‘<release number>’

הוספת קישור לגראב

כעת לפני שנאתחל את המערכת עלינו להוסיף קישור לגראב על ידי

cd /boot/grub

nano menu.lst

כל מה שנותר זה לאתחל את המערכת לקרנל החדש ולהתחיל להשתעשע…

]]> http://vaxincity.wordpress.com/2010/06/30/kernel_source_makeing/feed/ 11 vaxincity http://vaxincity.wordpress.com/2010/06/30/made-in-dizengoff/ http://vaxincity.wordpress.com/2010/06/30/made-in-dizengoff/#comments Wed, 30 Jun 2010 16:51:03 +0000 vaxincity http://vaxincity.wordpress.com/?p=475 ]]>

לפני כמה ימים בשעת צהריים של יום שישי עם כוס בירה צוננת ביד יצא לי לראות את אחד ההרכבים הקבארטיים הבודדים (אם לא היחיד) שעלו על במות תל אביב בשנים האחרונות.

מדובר על הצמד made in dizengoff שהחל לפעול בתל אביב בחודשים האחרונים. שי הררי (שירה וגיטרה) מספר כי ההרכב צמח מתוך צורך אמיתי למוזיקה אחרת ושונה מהקו השולט כיום בארץ. זמן רב חיפש צלע נוספת לתפקיד השירה וכאשר מצא את שרית גרינברג (שירה וקאבארטיזציה) היה זה מעל ומעבר למה שציפה. שרית תרמה חומרים כתובים בעלי הווי יונה וולכי שכזה, ללא פחד לעסוק בנושאים שרבים לא רוצים לגעת בהם וגם סגנון שירה מאוד ייחודי.

במהלך העבודה המשותפת שלהם נבנה סגנון במה ושירה שמזכירים הרכבי קונספט משנות השמונים. שרית (שהיא שחקנית במקצועה) ושי (שהוא מנהל שיווק)  נעטפים בתחפושת הן בפן ההתנהגותי והן בפן הלבושי. הפרוייקט המשותף נלקח צעד קדימה כאשר הצליחו לגייס צלמת מקצועית (הדר דוהר – קרדיטים על כל התמונות הנפלאות) שתתעד את ההופעות ועוד משהו שהוטלה עלי צנזורה לספר עליו, ממש ברגע האחרון. שי מציין שכעת הם עובדים עם מתופף חדש שיצטרף בקרוב להופעות. בנוסף הם מחפשים אחר נגן גיטרה מעניין שיוכל להשתלב ברוח ההרכב המתגבש.

“ככל שיצטרפו יותר כלים, הקונספט ישתנה ויתאים עצמו להווי אשר יצמח” אומר שי ומוסיף “שאין ספק שגם המופע יראה שונה ומושקע יותר ככל שהזמן יעבור”.

בנתיים שיתוף הפעולה ביניהם זורם גם לכיוונים של אחרים. לאחרונה שי הצטרף בתור נגן ליווי להצגה מצליחה ששרית מעלה מזה מספר שנים והיד עוד נטויה. עושה רושם שיש חיבור טוב ביניהם על הבמה.

במהלך ההופעה צילמתי כמה קטעי וידאו לשם התרשמות. אני מתנצל מראש על האיכות וזאת מכיוון שההקלטה נעשתה על ידי מצלמה שלא מיועדת לעניין.

למעוניינים להצטרף להרכב או להתעדכן בנוגע להופעות קרובות יש אפשרות לדבר עם האנשים הנפלאים הללו דרך הפייסבוק (בצעו חיפושי שמות באנגלית).

בנתיים שיהיה לכולנו המשך שבוע טוב.

]]> http://vaxincity.wordpress.com/2010/06/30/made-in-dizengoff/feed/ 2 vaxincity dizengoff in action נותנים את הנשמה http://vaxincity.wordpress.com/2010/06/20/proj_ker/ http://vaxincity.wordpress.com/2010/06/20/proj_ker/#comments Sun, 20 Jun 2010 20:43:59 +0000 vaxincity http://vaxincity.wordpress.com/?p=471 ]]> כבר זמן מה אני מנסה ללמוד ,לסירוגין, כתיבת מודולים ודרייברים בלינוקס (בעיקר על גבי דביאן). כל פעם אני נתקל בבעיה אחרת שגורמת לי לבזבז זמן רב ומדי פעם לנטוש את הלמידה לכמה שבועות או חודשים.

אני מנצל את הפלטפורמה הזו בכדי לבקש ממי שמעוניין ללמוד ,בשילוב כוחות עימי, את התחום הזה ליצור איתי קשר בכתובת הדואר האלקטרוני שבדף הפרטים עלי.

המשך שבוע נעים

רציתי ליצור משווה כחלק מפרוייקט גדול יותר. המשווה שלי בודק 32 סיביות, מה שאילץ אותי לחשוב על פתרון מהיר ופשוט לכתיבת הקוד על ידי וקטורים. במרבית שפות התכנות ניתן ליצור וקטור מסויים במקום לכתוב סיבית אחר סיבית.

פירוט הקוד

מדובר על קוד מאוד פשוט ואלגנטי, אשר למעשה מתבסס על חלק ה assign שבו מתבצעת פעולת ההשוואה ממש וחלק ה parameter שבו אנו מגדירים את הוקטור.

module comp_32

parameter word_size = 32;

input [word_size-1:0]

output A_greater_B, A_little_B, A_equal_B;

assign A_greater_B = (A>B),

A_little_B = (A<B),

A_equal_B = (A==B);

endmodule

עיבוד הרעיון לשימוש שונה

כעת נניח שאנו רוצים משווה שיבדוק לנו רגישות או שינויים כלשהם בכל סיבית וסיבית (על פני שתי סיביות). נוכל לקחת את הקוד הקודם ובעזרת שינויים מינוריים להלביש אותו על הדרישה

module comp_32_edg_sense

input A1,A0,B1,B0;

output A_greater_B, A_little_B, A_equal_B;

reg A_greater_B, A_little_B, A_equal_B;

always @(A0 or A1 or B0 or B1)

begin

A_little_B = ({A1,A0}<{B1,B0});

A_greater_B = ({A1,A0}>{B1,B0});

A_equal_B = ({A1,A0}=={B1,B0});

endmodule

כמה נקודות חשובות

כעת נעלה ברמת החשיבה. נניח והרכיב הוא חלק ממערך גדול יותר של רכיבים. מה יקרה בעת הפעלת רכיב הספציפי הזה? אות התחלתי (ובלתי רצוי) עלול להיפלט במוצא !

תיקון הקוד יעשה על ידי איפוס של כל המוצאים בגוף המודול. פרט לכך נלך צעד קדימה ונשכלל את קוד ההפעלה על פי עקרון לולאת “אם”

module comp_32_edg_sense

input [1:0] A,B;

output A_greater_B, A_little_B, A_equal_B;

reg A_greater_B, A_little_B, A_equal_B;

always @(A or B)

begin

A_little_B = 0;

A_greater_B = 0;

A_equal_B = 0;

if (A==B)

A_equal_B=1;

else if (A>B)

A_greater_B=1;

else

A_little_B=1;

end

endmodule

הערה מקצועית

בדוגמא האחרונה נשים לב לכך שהלולאה בנויה בצורה כזו שאם לא מתקיימים א וב נגיע למצב האחרון כברירת מחדל (A קטן מ B). במערכות קריטיות שאין שום מרווח לטעות ולו הקטנה ביותר, הייתי ממליץ לסגור היטב גם את האופציה הזו רק לשם הסדר הטוב.