هفتِ بیجار

[coronaa]

آخه این که مختصات یک خط سه بعدیه نه یک حجم سه بعدی! باز هم دوباره میگم اگه شما با داده های اصلی که دارید کار می کنید از try استفاده کنید دیگه ارور نمیده

داده های من همین چیزی هست که گذاشتم. با try خطا نمیده ولی چیزی هم پلات نمیکنه!

Kp =linspace(0,60,50)';

Ki =linspace(0,200,50)';

Kd =linspace(0,1000,50)';

% -----------------



NODES = [Kp Ki Kd];

try

plotSlices(NODES,(-1:1))

catch

end

شما داده ها رو یکدفعه تولید کردین ولی من جدا جدا داشتم و پیش هم گذاشتم. چه فرقی میکنه؟
من تصویر اصلی رو هم از اینجا گرفتم.

[rahnema1]

مگه نگفتید یکی از پارامترها به این صورته؟
https://barnamenevis.org/attachment.p...2&d=1430328678
یه کم دارم گیج می شم

[coronaa]

مگه نگفتید یکی از پارامترها به این صورته؟
https://barnamenevis.org/attachment.p...2&d=1430328678
یه کم دارم گیج می شم

درسته ولی من فکر کنم مقاله اومد فقط رنج پارامترها رو بصورت نموداری نشون داد. دیگه نیومد نقطه به نقطه رو پلات کنه!
همون قسمتی که های لایت کردم. منم میخوام بگم رنج این سه تا پارامتر در این محدوده هست و پلات بشه. مگه اون نمودار اصلی اینو نشون نمیده؟

[rahnema1]

شما بالاخره تونستید پلات را به صورت سه بعدی بچرخونید؟
چرا این که میگم انجام نمیدید؟ لطفا با try اون سه تا پارامتر را با تابع plotSlices که سطوح اون 10:10:90 باشه پلات کنید. خواهشا
با drawMesh هم پلات کنید و بچرخونید و تصویر سه بعدی را اینجا بذارید

[coronaa]

منظورت کدوم پارامترها بود؟
من با این داده ها انجام دادم :

clc; clear; close all;



Kp =linspace(0,60,50)';

Ki =linspace(0,200,50)';

Kd =linspace(0,1000,50)';



NODES = [Kp Ki Kd];



try

    plotSlices(NODES, (10:10:90))

catch

end

% -------------------------------

figure

FACES = minConvexHull(NODES);

try

    drawMesh(NODES,FACES)

catch

end

خطای زیر رو داد :

??? Error using ==> qhullmx

qhull precision error: initial facet 3 is coplanar with the interior point



ERRONEOUS FACET:





While executing:  | qhull Qt



Options selected for Qhull 2003.1 2003/12/30:

  Qtriangulate  _pre-merge  _zero-centrum  _max-width 1e+003

  Error-roundoff 1.1e-012  _one-merge 7.5e-012  _near-inside 3.7e-011

  Visible-distance 2.1e-012  U-coplanar-distance 2.1e-012

  Width-outside 4.3e-012  _wide-facet 1.3e-011





The input to qhull appears to be less than 3 dimensional, or a

computation has overflowed.



Qhull could not construct a clearly convex simplex from points:





The center point is coplanar with a facet, or a vertex is coplanar

with a neighboring facet.  The maximum round off error for

computing distances is 1.1e-012.  The center point, facets and distances

to the center point are as follows:









facet

 p42

 p49

 p0

 distance= -7.1e-015



facet

 p9

 p49

 p0

 distance= -7.1e-015



facet

 p9

 p42

 p0

 distance=    0



facet

 p9

 p42

 p49

 distance= -1.4e-014





These points either have a maximum or minimum x-coordinate, or

they maximize the determinant for k coordinates.  Trial points

are first selected from points that maximize a coordinate.





The min and max coordinates for each dimension are:



  0:         0        60  difference=   60



  1:         0       200  difference=  200



  2:         0      1000  difference= 1000





If the input should be full dimensional, you have several options that

may determine an initial simplex:

  - use 'QJ'  to joggle the input and make it full dimensional

  - use 'QbB' to scale the points to the unit cube

  - use 'QR0' to randomly rotate the input for different maximum points

  - use 'Qs'  to search all points for the initial simplex

  - use 'En'  to specify a maximum roundoff error less than 1.1e-012.

  - trace execution with 'T3' to see the determinant for each point.



If the input is lower dimensional:

  - use 'QJ' to joggle the input and make it full dimensional

  - use 'Qbk:0Bk:0' to delete coordinate k from the input.  You should

    pick the coordinate with the least range.  The hull will have the

    correct topology.

  - determine the flat containing the points, rotate the points

    into a coordinate plane, and delete the other coordinates.

  - add one or more points to make the input full dimensional.







Error in ==> convhulln at 71

[k,vv] = qhullmx(x', opt);



Error in ==> minConvexHull at 44

hull = convhulln(nodes);



Error in ==> plot4 at 15

FACES = minConvexHull(NODES);

[coronaa]

به اینصورت اجرا کردم :

clc; clear; close all;



NODES = randn(50,3);



try

    plotSlices(NODES, (10:10:90))

catch

end

% -------------------------------

figure

FACES = minConvexHull(NODES);

try

    drawMesh(NODES,FACES)

catch

end

که دستور اولی رو پلات نکرد (خطایی هم نداد) و دستور دومی رو بصورت زیر پلات کرد :

plot5.jpg

[coronaa]

به این صورت که انجام دادم :

clc; clear; close all;



NODES = randn(50,3);



plotSlices(NODES, (-1:1))

% -------------------------------

figure

FACES = minConvexHull(NODES);

try

    drawMesh(NODES,FACES)

catch

end

هر دو رو پلات کرد ولی level اون چیزی که خواستین نیست. با اون level خطا میده که توی پست قبلی انجام دادم.

1.jpg

2.jpg

[coronaa]

من میخوام با این رنج های زیر نموداری مثل مقاله رسم بشه :

Kp =linspace(0,60,50)';

Ki =linspace(0,200,50)';

Kd =linspace(0,1000,50)';

میشه؟؟؟

level رو میشه هر طوری انتخاب کرد؟

[rahnema1]

همه چیز را ترسیم کردید به جز اون چیزی که باید ترسیم کنید :)
از linspace استفاده نکنید
سه تا بردار دارید که یکی از اونها به این شکل هست:
https://barnamenevis.org/attachment.p...2&d=1430328678
اون سه تا را به عنوان سه ستون NODES در نظر بگیرید و با plotSlices و drawMesh ترسیم کنید و بچرخونید و حتما از try هم استفاده کنید

[coronaa]

مگه توی اون تصویر اصلی نیومد فقط رنج پارامترها رو نشون بده؟ فکر نکنم اومده باشه که نقاط بین این دو رنج رو هم پلات کرده باشه!
به هر حال من همون چیزی که گفتید رو قرار دادم.

clc; clear; close all



load Data_Conventional



Kp =kd1co;

Ki =kp1co;

Kd =ki1co;

% -----------------



NODES = [Kp Ki Kd];

try

    plotSlices(NODES,(-1:1))

catch

end



figure

FACES = minConvexHull(NODES);

try

    drawMesh(NODES,FACES)

catch

end

قسمت اولی رو پلات نمیکنه ولی قسمت دومی بصورت زیر شد :

rrrr.jpg

[coronaa]

مهندس ببخشید من هنوز نتونستم کامپایلر رو نصب کنم. یعنی sdk که گفته بودی رو هم دانلود کردم ولی وسط راه یه پیغامی میده به این صورت :
ضاهراً یه فایل وجود نداره. درسته؟

problem.jpg


problem2.jpg

[rahnema1]

اون -1:1 برای اون مثال دیگه بود واسه داده های خودتون اینجور بذارید

    plotSlices(NODES,(5:5:50))

واسه نصب هم گزینه دوم را انتخاب کنیدو مسیر فولدر sdk را معرفی کنید

[coronaa]

اون -1:1 برای اون مثال دیگه بود واسه داده های خودتون اینجور بذارید

2	plotSlices(NODES,(5:5:50))

انجام دادم. اتفاق خاصی نیوفتاد!
من داده ها رو توی اکسل ذخیره کنم و اینجا آپ کنم میتونید خودتون با اکتاو بررسی کنید؟

واسه نصب هم گزینه دوم را انتخاب کنیدو مسیر فولدر sdk را معرفی کنید

تصویر دومی که گذاشتم واسه همین بود. وقتی از مسیر کامپیوتر فراخونی میکنم خطای تصویر دوم رو میده.

[rahnema1]

ظاهرا دانلود نکردید لینکش داخل همون فایل
Win8.1SDKFiles.txt هست
http://download.microsoft.com/downlo....29750/Windows Software Development Kit-x86_en-us.msp

[rahnema1]

داده ها را بفرستید

[coronaa]

داده ها را بفرستید

http://www.sharefile.ir/uploads/1431184856.xlsx

تشکر

[rahnema1]

http://www.sharefile.ir/uploads/1431184856.xlsx

تشکر

شما level را به صورت 22:3:47 وارد کنید نشون میده

[coronaa]

درست شد.
ولی خب دلیلش چیه که level رو این اعداد انتخاب کردین و بر چه اساسی؟
دوم اینکه از روی نمودار حالا چطوری میشه فهمید که Kp به سمت چه عددی میل کرد؟

ppp.jpg

[rahnema1]

گفتم که از min و max استفاده کنید تا محدوده متغیر ها را متوجه بشید از روی نمودار mesh یا نمودار نقطه ای هم میشد فهمید محدوده اونها چه قدره
اینکه kp به سمت چه عددی میل کرده نمیدونم از روی نمودار که مشخصه

[coronaa]

گفتم که از min و max استفاده کنید تا محدوده متغیر ها را متوجه بشید از روی نمودار mesh یا نمودار نقطه ای هم میشد فهمید محدوده اونها چه قدره
اینکه kp به سمت چه عددی میل کرده نمیدونم از روی نمودار که مشخصه

ممنونم.
من منظورم چیز دیگه ای بود. فرص کنید کدی در دسترس نیست و خواننده میواد از روی نمدار متوجه بشه که رنجشون چیه!
مثل همون چیزی که توی مقاله هست. پس اصلن واسه چی این نمودار رو رسم کرده بود؟
بازم دستتون درد نکنه. خیلی زحمت کشیدین.

[coronaa]

سلام
مهندس ببخشید. اون لینکی که واسه sdk گذاشته بودین رو دانلود کردم و توی فولدر مربوطه قرار دادم ولی همچنان همون خطا رو میده!
من فایلهایی که دانلود کردم به صورت زیره :

1.jpg

فایلی که شما گفتید، شماره 1 هست و چیزی که ظاهرن موقع نصب میخواد فراخونی کنه شماره 2 هست. فایله وجود داره ولی خطا میده.
نمیشه کامپایلر دیگه ای نصب کرد؟

3.jpg

[rahnema1]

نمیدونم مشکل کجاست. اگه بتونید توی تنظیمات folder option در ویندوز گزینه نشون دادن extension ها رافعال کنید تا مشخص بشه پسوند فایل exe هست یا msu
من خودم تا به حال این را نصب نکردم توی اون سایت لینکش را دیدم گفتم شاید به درد شما بخوره حالا اگه تو نصبش مشکل دارید فکر کنم بهتره یه عدد dvd ویژوال استادیو کامل البته همون نسخه ای که می خواهید فکر کنم 2012 بود را یه جا پیدا کنید که این سی دی ها همه جا ریخته و فراوان هست

[coronaa]

نمیدونم مشکل کجاست. اگه بتونید توی تنظیمات folder option در ویندوز گزینه نشون دادن extension ها رافعال کنید تا مشخص بشه پسوند فایل exe هست یا msu
من خودم تا به حال این را نصب نکردم توی اون سایت لینکش را دیدم گفتم شاید به درد شما بخوره حالا اگه تو نصبش مشکل دارید فکر کنم بهتره یه عدد dvd ویژوال استادیو کامل البته همون نسخه ای که می خواهید فکر کنم 2012 بود را یه جا پیدا کنید که این سی دی ها همه جا ریخته و فراوان هست

سلام مهندس عزیز
بالاخره پس از روزها تلاش و کوشش تونستم بر طرف کنم. یه ++c دیگه پیدا کردم و نصب کردم که جواب داد خدا رو شکر.
دست شما هم درد نکنه
:)

[coronaa]

سلام آقای مهندس. وقت شما بخیر.
ببخشید یه مشکل دیگه داشتم. من فرمول زیر رو میخوام کد بزنم و پلات کنم :

soc.jpg

بصورت زیر نوشتم ولی توی لحظه ی صفر خطا میگیره :

clc; close all;



tic

Time=linspace(1,10,250000);

P=rand(1,250000);

soc=zeros(1,250000);

soc(0)=100;

for i=1:length(Time)

soc(i)=(soc(i-1)+(P(i)*(Time(i)-Time(i-1))))/100;

end

socc=soc';

plot(Time,socc);

toc

خطا :

??? Attempted to access soc(0); index must be a positive integer or

logical.



Error in ==> Soc at 7

soc(0)=100;

[rahnema1]

سلام
این را امتحان کنید

soc(2:end) = (soc(1:(end -1)) + P(2:end) .*deta(2:end)) * 0.01;

soc = cumsum(soc);

[coronaa]

deta الآن چیه؟
البته من از لحظه ی اول میخوام. وگرنه از لحظه ی دوم رو بدست میارم. وقتی که بصورت زیر باشه :

soc(1)=100;

for i=2:length(Time)

soc(i)=(soc(i-1)+(P(i)*(Time(i)-Time(i-1))))/100;

end

[rahnema1]

deta الآن چیه؟

تو همون فرمول دلتا داریم دیگه بر اساس همون عکس فرمول که گذاشتید من نوشتم

[coronaa]

تو همون فرمول دلتا داریم دیگه بر اساس همون عکس فرمول که گذاشتید من نوشتم

درسته ولی دلتا رو که شما تعریف نکردین!
مگه متلب این عبارت رو میشناسه؟

[coronaa]

سلام
این را امتحان کنید

soc(2:end) = (soc(1:(end -1)) + P(2:end) .*deta(2:end)) * 0.01;

soc = cumsum(soc);

من اصلن متوجه نشدم که شما سمت چپ معادله رو چرا برداری نوشتین؟؟ ما که از همون اول مقادیر soc رو نداریم!
این رابطه بازگشتیه و باید با مقدار فبل خودش حساب بشه.

[rahnema1]

این را امتحان کنید

P = rand(1,250000);

soc=zeros(1,250000);

soc(1) = 100;

delta = diff(1:250000);

soc(2:end) = (soc(1:(end -1)) + P(2:end) .* delta) * 0.01;

soc = cumsum(soc);

[coronaa]

این را امتحان کنید

P = rand(1,250000);

soc=zeros(1,250000);

soc(1) = 100;

delta = diff(1:250000);

soc(2:end) = (soc(1:(end -1)) + P(2:end) .* delta) * 0.01;

soc = cumsum(soc);

مرسی. همونی شد که توی پست 146 گذاشتم.

[coronaa]

سلام مهندس
وقتت بخیر.
ببخشید من یه حلقه دارم که هفت بار اجرا میشه و در هر مرحله (سطر) تعداد ستون های متفاوتی داره. بصورت تصویر زیر.
من چطوری باید سلولش رو تعریف کنم؟ باید یه سلول 7 در 3 بشه. چون نهایتن تعداد ستونها به 3 میرسه.
البته جای اون حروفها اعداد قرار میگیره.

aaa.jpg

[rahnema1]

سلام
درست متوجه نشدم ولی مثلا در این مثال x چیزی هست که قراره در هر سطر جا داده بشه

jadval = {};

for i = 1 : 7

	jadval{i} = x;

end

[coronaa]

برنامه ی من 7 بار با حلقه for اجرا میشه.
در حلقه اول، یه مقدار وجود داره که در سطر اول و ستون اول قرار میگیره. یعنی در سلول a
در حلقه دوم و سوم هم همینطور جای b و c قرار میگیرن.
در حلقه چهارم، دو مقدار بدست میاد که در سطر چهارم و ستونهای اول و سوم قرار میگیرن. یعنی جای d , e
...
درحلقه هفتم، برای x سه مقدار داریم. یعنی یک بردار سه 1 در سه هست، که در سطر 7 و ستون یک تا سه قرار میگیره.
اینو چطوری باید توی یک سلول و یا آرایه ذخیره کنم؟
آرایه باشه بهتره. منظورم بصورت ماتریسی هست.

[rahnema1]

خب یک آرایه خالی به همین ابعاد 7 در 3 ایجاد کنید و در حلقه مقدار دهی کنید . الان من درست متوجه نمیشم

[coronaa]

میخوام مقدار offer توی یک آرایه قرار بگیره. مثل شکلی که قرار دادم :

alpha0=[.4,.52,.6];

beta0=[1.2 1.2 1.2];

lambda=1;

D0=[5.3,3.6,1.3];

m0=[1,1,1];

k=.38;



% % -----------------------------

n=[1 2 3];

for i=1:numel(n)

    S{i}=nchoosek(n,i);

end

% % -----------------------------

Coalition=1;

profit=[];

 jadval = zeros(7,3);

for i=1:numel(S)

    for j=1:size(S{i},1)

        alpha=[];

        beta=[];

        

        for kk=1:size(S{i},2)

            alpha=[alpha ones(1,m0(S{i}(j,kk))).*alpha0(S{i}(j,kk))];

        end

        for kk=1:size(S{i},2)

            beta=[beta ones(1,m0(S{i}(j,kk))).*beta0(S{i}(j,kk))];

        end

        m=sum(m0(S{i}(j,:)));

        D=sum(D0(S{i}(j,:)));

        fun_MG0 = @(q,alpha,beta,k,D) -sum((2*alpha.*q+beta + k .*(D-sum(q))).*q -alpha.*(q).^2-beta.*q);

        fun_MG = @(q) fun_MG0(q,alpha,beta,k,D) ;

        lb=zeros(1,m); % LowerBound

        ub=ones(1,m);  % UpperBound

        [x1, fx1]=fmincon(fun_MG,ones(1,m)*.99,[],[],[],[],lb,ub);

        profit(Coalition)=(-fx1);

        Coalition = Coalition + 1;

        fprintf('\n Emount of Energy :\n');

        disp(x1);

        fprintf('\n---------------------------------------\n');

        fprintf('\nOffer Without Saticfiction Function :\n');

        disp(2*alpha.*x1+ beta);

        fprintf('\n---------------------------------------\n');

        fprintf('\nOffer With Saticfiction Function :\n');

        offerr=2*alpha.*x1+ beta + k .*(D- sum(x1));

        disp(2*alpha.*x1+ beta + k .*(D- sum(x1)));

        fprintf('\n---------------------------------------\n');

         

    end 

end

[rahnema1]

فکر کنم اینجور بذارید درست می شه

offer{i}(j) = 2*alpha.*x1+ beta + k .*(D- sum(x1));

[coronaa]

??? Subscripted assignment dimension mismatch.



 Error in ==> bbbbbbbb at 46

        offer{i}(j) = 2*alpha.*x1+ beta + k .*(D- sum(x1));

[rahnema1]

این یکی

alpha0=[.4,.52,.6];

beta0=[1.2 1.2 1.2];

lambda=1;

D0=[5.3,3.6,1.3];

m0=[1,1,1];

k=.38;



% % -----------------------------

n=[1 2 3];

for i=1:numel(n)

    S{i}=nchoosek(n,i);

end

% % -----------------------------

Coalition=1;

profit=[];

 jadval = zeros(7,3);

 offer={};

for i=1:numel(S)

offerr_array=[];

    for j=1:size(S{i},1)

        alpha=[];

        beta=[];

        

        for kk=1:size(S{i},2)

            alpha=[alpha ones(1,m0(S{i}(j,kk))).*alpha0(S{i}(j,kk))];

        end

        for kk=1:size(S{i},2)

            beta=[beta ones(1,m0(S{i}(j,kk))).*beta0(S{i}(j,kk))];

        end

        m=sum(m0(S{i}(j,:)));

        D=sum(D0(S{i}(j,:)));

        fun_MG0 = @(q,alpha,beta,k,D) -sum((2*alpha.*q+beta + k .*(D-sum(q))).*q -alpha.*(q).^2-beta.*q);

        fun_MG = @(q) fun_MG0(q,alpha,beta,k,D) ;

        lb=zeros(1,m); % LowerBound

        ub=ones(1,m);  % UpperBound

        [x1, fx1]=fmincon(fun_MG,ones(1,m)*.99,[],[],[],[],lb,ub);

        profit(Coalition)=(-fx1);

        Coalition = Coalition + 1;

        fprintf('\n Emount of Energy :\n');

        disp(x1);

        fprintf('\n---------------------------------------\n');

        fprintf('\nOffer Without Saticfiction Function :\n');

        disp(2*alpha.*x1+ beta);

        fprintf('\n---------------------------------------\n');

        fprintf('\nOffer With Saticfiction Function :\n');

        offerr_array(j)=2*alpha.*x1+ beta + k .*(D- sum(x1));

        disp(2*alpha.*x1+ beta + k .*(D- sum(x1)));

        fprintf('\n---------------------------------------\n');

         

    end 

    offer{i} = offerr_array;

end

[coronaa]

???  In an assignment  A(I) = B, the number of elements in B and

 I must be the same.



Error in ==> aaaaaaaaaaaaa at 46

        offerr_array(j)=2*alpha.*x1+ beta + k .*(D- sum(x1));