ปัญหาของผมคือผมมีข้อมูล

ผมต้องการให้เป็นแบบนี้

ในคู่มือของ MATLAB บอกว่าให้ใช้คำสั่ง reshape ซึ่งมีวิธีการใช้ดังนี้

B = reshape(A,m,n)
โดยที่ A คือเมตริกซ์ที่ต้องการจะทำการเปลี่ยนแปลง
m และ n คือจำนวนแถวและหลัก

ถ้าผมใช้ B = reshape(y,[],1)  ผมจะได้

ซึ่งไม่ใช่สิ่งที่ต้องการ คิดอยู่ต้องนานว่าจะทำอย่างไร สรุปว่าก็แค่ B = reshape(y’,[],1) นั่นก็
คือทำการ transpose ก่อนนั่นเอง

Bila Anda ingin mendapatkan paper lengkap, Silahkan Submit permintaan Anda melalui menu "Comment" di bawah ini.

In an example in this post, we use Trapezoidal rule to solve an improper integral where the integrand becomes infinite.  The integral is $latex \int_{0}^{b} 1/sqrt{x} dx $.  The integrand becomes infinite at x=0.  Since x=0 would be one of the points where the integrand will be sought by the multiple-segment Trapezoidal rule, we choose the value of the integrand at x=0 to be zero (any other value would do too - a better assumption would be f(h), where h is the segment width in the multiple-segment Trapezoidal rule).

Here is a MATLAB program that shows you the exact value of the integral and then compares it with the multiple-segment Trapezoidal rule.  The convergence is slow but you can integrate improper integrals using Trapezoidal rule.

The MATLAB program that can be downloaded at http://numericalmethods.eng.usf.edu/blog/trapezoidal_improper_sqrt.m (better to download it as single quotes from the web-post do not translate correctly with the MATLAB editor).  The html file showing the mfile and the command window output is here: http://numericalmethods.eng.usf.edu/blog/trapezoidal_improper_sqrt.html

% Simulation : Can I use Trapezoidal rule for an improper integral?

% Language : Matlab 2007a

% Authors : Autar Kaw, http://numericalmethods.eng.usf.edu

% Mfile available at
% http://numericalmethods.eng.usf.edu/blog/trapezoidal_improper_sqrt.m

% Last Revised : October 8, 2008

% Abstract: This program shows use of multiple segment Trapezoidal
% rule to integrate 1/sqrt(x) from x=0 to b, b>0.

clc
clear all

disp('This program shows the convergence of getting the value of ')
disp('an improper integral using multiple segment Trapezoidal rule')
disp('Author: Autar K Kaw.  autarkaw.wordpress.com')

%INPUTS.  If you want to experiment, these are the only variables
% you should and can change.
% b  = Upper limit of integration
% m = Maximum number of segments is 2^m
b=9;
m=14;

% SIMULATION
fprintf('\nFinding the integral of 1/sqrt(x) with limits of integration as x=0 to x=%g',b)

% EXACT VALUE OF INTEGRAL
% integrand 1/sqrt(x)
syms x
f=1/sqrt(x);
a=0;
valexact=double(int(f,x,a,b));
fprintf('\n\nExact value of integral = %f',valexact)
disp( '  ')

f=inline('1/sqrt(x)');
%finding value of the integral using 16,...2^m segments
for k=4:1:m
n=2^k;
h=(b-a)/n;
sum=0;
for i=1:1:n-1
sum=sum+f(a+i*h);
end
% See below how f(a) is not added as f(a)=infinity.  Instead we
% use a value of f(a)=0.  How can we do that? Because as per integral calculus,
% using a different value of the function at one point or
% at finite number of points does not change the value of the
% integral.
sum=2*sum+0+f(b);
sum=(b-a)/(2*n)*sum;
fprintf('\nApproximate value of integral =%f with %g segments',sum,n)
end
disp('  ')

This post is brought to you by Holistic Numerical Methods: Numerical Methods for the STEM undergraduate at http://numericalmethods.eng.usf.edu.

An abridged (for low cost) book on Numerical Methods with Applications will be in print (includes problem sets, TOC, index) on December 10, 2008 and available at lulu storefront.

Subscribe to the blog via a reader or email to stay updated with this blog. Let the information follow you.

Bila Anda ingin tulisan lengkap, Silahkan Submit permintaan Anda melalui menu "Comment" di bawah ini.

I have been playing with JMP a bit, but honestly, it's a bit too "high level" and while it's neat for data exploration, it really made me nuts that two weeks after I got a license, they were *offering* me a special deal on the impending upgrade. Nothing like dropping a bundle on instantly outdated software. Way to go guys.

So for me, it's been a pleasure to use R, Python (with iPython of course), and ferret on my Mac Pro. If I ever take a break from playing on the computer I'll post up what I installed on the new Mac in terms of scientific software.

We would also like you to be famaliar with IDL, MatLab, Mathematica and other scientific software langauges.  If you have a degree in physics or astronomy this is an advanatge.  You must be a US citizen.

Contact me at stevep at brainscouts dot com if you are interested.

Thanks
Steve Parker

MATLAB guide, 2nd edition--by Higham, Desmond J. & Higham, Nicholas J.

I give most of the credit for this information to my MATLAB professeur Dr. Sutherland from the University of Utah and Mcgraw hills in their text book.

We are seeking a Quality Assurance engineer to join our Budapest team in development of the next generation of equity analytics products. The candidate must have strong experience in test-planning and execution for enterprise software products. Practical knowledge of test strategies, test planning, test execution, and test automation is required. Candidate should be experienced in use of scripting tools like HP/Mercury QTP and VB, as well as with test management tools like TestDirector and Perforce. Experience with MATLAB is highly desirable.

The candidate will work closely with our global development team-members in Berkeley to develop highly quantitative applications products. The candidate must be fluent in English and have strong oral and written communication skills, interpersonal skills, and software development skills. He or she needs to be flexible and adaptable, and should have demonstrated success working in a team environment.

Mandatory Requirements

• Bachelors Degree in Computer Science or related discipline.
• Strong Quality Assurance engineering experience on the Windows XP/Vista platform.
• Experience with test automation, including use of HP/Mercury QTP and VB.
• 5-7 years of relevant experience in professional software development environment.
• High degree of motivation, with strong communications skills and ability to learn quickly and work independently.

Desirable Skills

• Knowledge of Equities and investment concepts.
• Experience with MATLAB.

laszlo_spiller@kellyservices.hu

BMI Categories:

BMI <19 (underweight)
19≤BMI≤25 (healthy)
25<BMI≤30 (overweight)
BMI is >30 (obese)

Here is a MATLAB program that calculates the BMI of a person, classifies him/her in the BMI category, and then suggests a healthy weight.

It is better to download the program as single quotes in the pasted version do not translate correctly when pasted into a mfile editor of MATLAB.  See the html version for clarity and sample output.

% So you want my phone number and my BMI?  How shallow can you get?

% Worksheet : Finding the BMI of  a person, classifying the BMI and
% suggesting a healthy weight
% Language : Matlab 2007a
% Authors : Autar Kaw
% Last Revised : September 27, 2008

%  Abstract: Finding the body mass index of a person, classifying their health
% and recommending a target weight.

%
%  In 1998, the federal government developed the body mass index (BMI)
%  to determine ideal weights.  Body mass index is calculated as
%  703 times the weight in pounds divided by the square of the height in inches,
%  the obtained number is then rounded to the nearest whole number
%  (23.5 will be rounded to 24; 23.1 will be rounded to 23; 23.52 will be rounded to 24).
%
% Write a MATLAB program to do the following:
%
% Assign a value to weight in lbs, and height in inches and
% then calculate BMI as a rounded integer.
% Output a variable called health_id as
%     0 if the person's BMI <19 (underweight)
%     1 if the person’s BMI is 19≤BMI≤25 (healthy)
%     2 if the person's BMI is 25<BMI≤30 (overweight)
%     3 if the person’s BMI is >30 (obese)
% Output also a variable hw for healthy weight in rounded integer lbs for all the conditions.
% Use fprintf command with explanation for the inputs and outputs
clc
clear all
% INPUTS
% Weight in lbs
weight=180;
% Height in inches
height=69;

%REST OF THE PROGRAM
bmi=weight/height^2*703.0;
bmi=round(bmi);

% Assigning the proper health_id and finding the suggested healthy weight
% health_id= BMI category
% hw = healthy weight
% You can use 4 separate if-end statements or switch case statement also.
if (bmi<19)
health_id=0;
hw=19*height^2/703;
elseif (bmi>=19 & bmi<=25);
health_id=1;
elseif (bmi>25 & bmi<=30);
health_id=2;
hw=25*height^2/703;
else
health_id=3;
hw=25*height^2/703;
end
hw=round(hw);

% Printing the outputs
disp('The BMI MATLAB program')
disp('Author: Autar K Kaw')
disp('http://numericalmethods.eng.usf.edu')
disp('September 28, 2008')
disp(' ')
fprintf('\nWeight of person =%6.0f lbs',weight)
fprintf('\nHeight of person =%6.0f inches',height)
fprintf('\nYour BMI is  =%g ',bmi)
if health_id==0
fprintf('\nYou are underweight.  \nYour target weight is =%6.0f lbs',hw)
elseif health_id==1
fprintf('\nYou are a healthy weight.  \nYour target weight is =%6.0f lbs',hw)
elseif health_id==2
fprintf('\nYou are overweight.  \nYour target weight is =%6.0f lbs',hw)
elseif health_id==3
fprintf('\nYou are obese.  \nYour target weight is =%6.0f lbs',hw)
end

________________________________________________________________________________________________

This post is brought to you by Holistic Numerical Methods: Numerical Methods for the STEM undergraduate at http://numericalmethods.eng.usf.edu.

An abridged (for low cost) book on Numerical Methods with Applications will be in print (includes problem sets, TOC, index) on December 10, 2008 and available at lulu storefront.

Subscribe to the blog via a reader or email to stay updated with this blog. Let the information follow you.

Then I just use that time value for drawing the head map function.

There are two different things to look at, and a couple of different issues to consider.
The first is the small task we were considering initially considering, i.e. looking at a smallish sqd file, translating it to matlab format, and then running some reduction on that.
There is also the broader (and ultimately more interesting) question, of how to take a complete sqd file and process it. In order to do that I will want to denoise a sqd file and create another sqd file which will contain the now-reduced data.
This allows us to run the noise reduction algorithm on the complete file before we extract any kind of trigger information or anything, reducing the possibility of removing signal even with pretty aggressive noise removal.

Nayef had a file that did this second thing, using his adaptive filter as the denoising technique, although in principle it should be possible to swap out which NR (noise reduction) algorithm I want to use while keeping his framework, thus allowing me to have a sqd file as both input and output.

Jonathan has given me the routines for doing TSPCA and SNS noise reduction, which he prefers. They do different things:

TSPCA uses the three reference channels and tries to remove the noise component from the signal channels by comparing them to some combination of the three reference channels. In its more complete form you?re allowed to include not only a linear combination of the three reference channels, but actually a phase delayed or advanced version of those channels. This is to take into account of the fact that the cables have different lengths, etc. Thus by including these various leads and lags we get a much better noise suppression than we would otherwise.

SNS does not use the reference channels at all. It simply looks for aberrations in a given signal channel caused by glitches in the system. If there is a major aberration in one channel without there being any corresponding deviation in any of the neighboring channels, it is considered a glitch and corrected.

There is also the low-pass filter (from MEG160) which we will make using the ellip function. That has a number of parameters which we will allow Matlab to tune in order to give us 20dB suppression at 30Hz. 4 and 0.5 are reasonable values for the first two parameters (a 4th order polynomial, 0.5dB of peak to peak ripple, whatever that may mean).

Finally there is de-meaning, which takes out the DC component.

The file that does this is called sqd2mat.m and is stored in /User/dbhertz/Documents/Matlab/dbhertz/

That will be where I plan on storing my code, at least for the present.

I have actually changed that file so that it now takes in one file and turns it into one .mat file. That?s more efficient and useful, at least as long as I?m dealing with small files (which is what I?m doing for the present. I?m not sure why Nayef was splitting things up into 160 files before, but this way seems better for what I want to do, in any case.

So what do I use in order to open a sqd file and translate its contents into a matlab array?

As far as I can tell, some basic i/o functions are handled by the four files written by Nayef:
mat2sqd.m
sqd2mat.m
adapt_noise_supp_fix.m
adapt_noise_supp_mref_fix.m

The last two seem to be variants on each other, although the specifics of that are a bit beyond me.

The first two seem to be (templates, at least) for translating a file from the MEG160 to the matlab format and the other way.

So sqd2mat.m seems to address point 1 from the last post.

כמו כל סטודנט טוב הלכתי ל"אופן אופיס" הטוב והתחלתי לכתוב. הבעיה נתעוררה בכמה אופנים:

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

• מבחינה ויזואלית התוצר הסופי לא נראה רציני, העבודה נראית כמו עבודה של תיכוניסט. לא טוב.

חלופות?

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

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

אז מה זה לאטך בדיוק?

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

לשם פשטות אתן כמה דוגמאות לכוחה של תוכנה זו:

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

מאיפה משיגים את זה??

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

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

איך זה עובד

אז בלאטך לא ממש כותבים כמו במעבד תמלילים. בלאטך נותנים פקודות. לדוגמא פתיחת מסמך בד"כ תגיע עם:

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

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

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

כמה דוגמאות על מה מדובר

קיבצתי כמה דוגמאות שממחישות את חוזקה של לאטך.

בעיות בלאטך

ויש כמה:

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

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

1. General Math Cheat Sheet (iPaper and other formats)
2. Elementary Algebra Cheat Sheet (PDF)
3. Trigonometry Cheat Sheet (PDF)
4. Calculus Cheat Sheet (PDF)
5. Derivatives and Integrals Cheat Sheet (PDF)
6. Laplace Transforms Cheat Sheet (PDF)
7. Abstract Algebra Cheat Sheet (PDF)
8. Probability Theory Cheat Sheet (PDF)
9. Matlab Cheat Sheet (PDF)
10. Mathematica Cheat Sheet (PDF)
11. Maple Cheat Sheet (PDF)
12. Maxima Cheat Sheet (HTML web page)
13. LaTeX Cheat Sheet (several formats)

I hope you find them helpful. I have personally used the LaTeX sheet often.

Check them out.

Update: Another source for math review is Pauls Online Math Notes.

So what did we do in MEG160?

1. Opened a file. This was just a small file, as I recall, although for MATLAB we?re going to be opening a larger file. I have not been able to find a larger file to try opening with MEG160 for the moment, so I?m just messing around with one small file with MEG160 for the moment. I?ll have to try opening larger files with MEG160 and doing this process with that at some point, since I?ll then have to do the following steps and include a sub-step there.
   • Include a trigger-finding procedure to extract the signals? I haven?t done this for the small file, as far as I can tell, unless it?s done automatically.
2. Remove the noise. This is just called ?Noise Reduction?. In MATLAB I?m going to have to call a couple of functions to do this.
3. Low Pass Filter (LPF) of 30 Hz.
4. Baseline Correction (using pre-test time of -100ms to 0ms).

So, this is what I need to do with my MATLAB script as well.

NOTE: MEG numbers the channels from 0-159. MATLAB numbers them as 1-160. MATLAB starts all arrays and everything at 1, not 0.

There are two components of the Noise Reduction.

1. TSPCA: Compares data channels to reference channels (158-160)