Ես փորձում եմ հասկանալ, թե հնարավո՞ր է Octave-ն ավելի արդյունավետ օգտագործել՝ հեռացնելով «for» օղակը, որը ես օգտագործում եմ X մատրիցայի յուրաքանչյուր տողի բանաձևը հաշվարկելու համար.
myscalar = 0
for i = 1:size(X, 1),
myscalar += X(i, :) * y(i) % y is a vector of dimension size(X, 1)
...
Բանաձևն ավելի բարդ է, քան սկալարի ավելացումը: Հարցն այստեղ իրականում այն է, թե ինչպես կարելի է կրկնել X տողերի միջով առանց ինդեքսի, որպեսզի ես կարողանամ վերացնել for հանգույցը:
X(i,:) վեկտորը y(i) սկալյարով: Ձեր myscalar-ն իրականում վեկտոր է: Եթե X-ը միայն մեկ սյունակ չունի, բայց ինչո՞ւ են այն դնում: Իսկ ինչո՞ւ եք փոխադրում բազմապատկման համար: Բացի այդ, կարիք չկա օգտագործել size (X, 1)-ը, Octave-ն արդեն ունի rows և columns ֆունկցիաները հենց այս նպատակների համար: 09.09.2012 myscalar-ը պետք է սկալյար լինի: Բազմապատկումների գումարո՞ւմ։ Մատրիցա՞ Ինչ չափի: Ես չեմ հասկանում, թե ինչ եք ուզում։ 09.09.2012 myscalar-ի հետ, քանի որ դուք պարզապես վերագրում էիք դրա վրա նախկինում, ոչ թե ավելացնում դրա վրա (=՝ +=-ի փոխարեն): 09.09.2012 Այո, դուք կարող եք օգտագործել հեռարձակումը դրա համար (դուք անհրաժեշտ է 3.6.0 կամ ավելի ուշ): Եթե դուք գիտեք python, սա նույնն է (բացատրություն python-ից) . Պարզապես բազմապատկեք մատրիցը սյունակով: Ի վերջո, cumsum-ն ավելացնում է, բայց մենք ուզում ենք միայն վերջին շարքը:
newx = X .* y;
myscalars = cumsum (newx, 1) (end,:);
կամ մեկ տողում՝ առանց ջերմաստիճանի փոփոխականների
myscalars = cumsum (X .* y, 1) (end,:);
Եթե չափերը ճիշտ են, հեռարձակումն իրականացվում է ավտոմատ կերպով: Օրինակ:
octave> a = [ 1 2 3
1 2 3
1 2 3];
octave> b = [ 1 0 2];
octave> a .* b'
warning: product: automatic broadcasting operation applied
ans =
1 0 6
1 0 6
1 0 6
octave> a .* b
warning: product: automatic broadcasting operation applied
ans =
1 2 3
0 0 0
2 4 6
Զգուշացման պատճառն այն է, որ դա նոր գործառույթ է, որը կարող է շփոթեցնել օգտատերերին և գոյություն չունի Matlab-ում: Դուք կարող եք ընդմիշտ անջատել այն՝ ավելացնելով warning ("off", "Octave:broadcast") ձեր .octaverc ֆայլին
Բոլորի համար, ովքեր օգտագործում են Octave-ի հին տարբերակը, նույնը կարող է իրականացվել՝ անմիջապես զանգահարելով bsxfunին:
myscalars = cumsum (bsxfun (@times, X, y), 1) (end,:);
bsxfun-ն իրականացվում է: Ես խմբագրել եմ իմ պատասխանը՝ ավելի հին տարբերակների համար այլընտրանք ցույց տալու համար: 06.12.2013