Trying to get function to iterate matrix












0















I am writing code to model a compression test. There are 211 xdata points and to make the formula work I need to calculate F at for each xdata point and use the value I1 for each iteration of F to put in the formula.



Equation for F



L is a value that varies from 1 to 0.8 depending on the xdata.
What I would like is to get a different value for I1 at each xdata point and then use that in the function.
This is working fine for L in the function, but I think because I1 is defined before the function starts there is an issue.



Any help would be greatly appreciated!



syms W1
syms W2
format long
format compact

data = xlsread('sample8.xlsx'); % two columns of data
xdata = data(:,1);
ydata = data(:,2);

Area = 133; % area of surface mm^2

Force = ydata; %mN

L = 1 - (xdata ./ 100); %Lambda
Lambda_inverse = L.^-1;

F = blkdiag(L.^-1, L.^-1, L.^2);

FT = F.';
B = F*FT;

I1 = trace(B);

fun = @(W,xdata) -Area.*( (W(1).* exp( (W(2)) .* (I1 - 3)) ) .* (L.^2 .* L.^-1) );

%starting guess
Wguess = [5,5];

[W,fminres] = lsqcurvefit(fun,Wguess,xdata,ydata)









share|improve this question

























  • The issue is not clear. What formula are you trying to use? Does your code crash or produce incorrect output? If it crashes, at which line, with what message and what is the type & size of input variables to the faulty line? If the output is incorrect, up to what point can you track the various variables to be correct and from what point are they bad?

    – Brice
    Nov 22 '18 at 8:48
















0















I am writing code to model a compression test. There are 211 xdata points and to make the formula work I need to calculate F at for each xdata point and use the value I1 for each iteration of F to put in the formula.



Equation for F



L is a value that varies from 1 to 0.8 depending on the xdata.
What I would like is to get a different value for I1 at each xdata point and then use that in the function.
This is working fine for L in the function, but I think because I1 is defined before the function starts there is an issue.



Any help would be greatly appreciated!



syms W1
syms W2
format long
format compact

data = xlsread('sample8.xlsx'); % two columns of data
xdata = data(:,1);
ydata = data(:,2);

Area = 133; % area of surface mm^2

Force = ydata; %mN

L = 1 - (xdata ./ 100); %Lambda
Lambda_inverse = L.^-1;

F = blkdiag(L.^-1, L.^-1, L.^2);

FT = F.';
B = F*FT;

I1 = trace(B);

fun = @(W,xdata) -Area.*( (W(1).* exp( (W(2)) .* (I1 - 3)) ) .* (L.^2 .* L.^-1) );

%starting guess
Wguess = [5,5];

[W,fminres] = lsqcurvefit(fun,Wguess,xdata,ydata)









share|improve this question

























  • The issue is not clear. What formula are you trying to use? Does your code crash or produce incorrect output? If it crashes, at which line, with what message and what is the type & size of input variables to the faulty line? If the output is incorrect, up to what point can you track the various variables to be correct and from what point are they bad?

    – Brice
    Nov 22 '18 at 8:48














0












0








0








I am writing code to model a compression test. There are 211 xdata points and to make the formula work I need to calculate F at for each xdata point and use the value I1 for each iteration of F to put in the formula.



Equation for F



L is a value that varies from 1 to 0.8 depending on the xdata.
What I would like is to get a different value for I1 at each xdata point and then use that in the function.
This is working fine for L in the function, but I think because I1 is defined before the function starts there is an issue.



Any help would be greatly appreciated!



syms W1
syms W2
format long
format compact

data = xlsread('sample8.xlsx'); % two columns of data
xdata = data(:,1);
ydata = data(:,2);

Area = 133; % area of surface mm^2

Force = ydata; %mN

L = 1 - (xdata ./ 100); %Lambda
Lambda_inverse = L.^-1;

F = blkdiag(L.^-1, L.^-1, L.^2);

FT = F.';
B = F*FT;

I1 = trace(B);

fun = @(W,xdata) -Area.*( (W(1).* exp( (W(2)) .* (I1 - 3)) ) .* (L.^2 .* L.^-1) );

%starting guess
Wguess = [5,5];

[W,fminres] = lsqcurvefit(fun,Wguess,xdata,ydata)









share|improve this question
















I am writing code to model a compression test. There are 211 xdata points and to make the formula work I need to calculate F at for each xdata point and use the value I1 for each iteration of F to put in the formula.



Equation for F



L is a value that varies from 1 to 0.8 depending on the xdata.
What I would like is to get a different value for I1 at each xdata point and then use that in the function.
This is working fine for L in the function, but I think because I1 is defined before the function starts there is an issue.



Any help would be greatly appreciated!



syms W1
syms W2
format long
format compact

data = xlsread('sample8.xlsx'); % two columns of data
xdata = data(:,1);
ydata = data(:,2);

Area = 133; % area of surface mm^2

Force = ydata; %mN

L = 1 - (xdata ./ 100); %Lambda
Lambda_inverse = L.^-1;

F = blkdiag(L.^-1, L.^-1, L.^2);

FT = F.';
B = F*FT;

I1 = trace(B);

fun = @(W,xdata) -Area.*( (W(1).* exp( (W(2)) .* (I1 - 3)) ) .* (L.^2 .* L.^-1) );

%starting guess
Wguess = [5,5];

[W,fminres] = lsqcurvefit(fun,Wguess,xdata,ydata)






matlab function matrix






share|improve this question















share|improve this question













share|improve this question




share|improve this question








edited Nov 21 '18 at 19:31









Wolfie

16.1k51744




16.1k51744










asked Nov 21 '18 at 18:56









James BlackwellJames Blackwell

96




96













  • The issue is not clear. What formula are you trying to use? Does your code crash or produce incorrect output? If it crashes, at which line, with what message and what is the type & size of input variables to the faulty line? If the output is incorrect, up to what point can you track the various variables to be correct and from what point are they bad?

    – Brice
    Nov 22 '18 at 8:48



















  • The issue is not clear. What formula are you trying to use? Does your code crash or produce incorrect output? If it crashes, at which line, with what message and what is the type & size of input variables to the faulty line? If the output is incorrect, up to what point can you track the various variables to be correct and from what point are they bad?

    – Brice
    Nov 22 '18 at 8:48

















The issue is not clear. What formula are you trying to use? Does your code crash or produce incorrect output? If it crashes, at which line, with what message and what is the type & size of input variables to the faulty line? If the output is incorrect, up to what point can you track the various variables to be correct and from what point are they bad?

– Brice
Nov 22 '18 at 8:48





The issue is not clear. What formula are you trying to use? Does your code crash or produce incorrect output? If it crashes, at which line, with what message and what is the type & size of input variables to the faulty line? If the output is incorrect, up to what point can you track the various variables to be correct and from what point are they bad?

– Brice
Nov 22 '18 at 8:48












0






active

oldest

votes











Your Answer






StackExchange.ifUsing("editor", function () {
StackExchange.using("externalEditor", function () {
StackExchange.using("snippets", function () {
StackExchange.snippets.init();
});
});
}, "code-snippets");

StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "1"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});

function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});


}
});














draft saved

draft discarded


















StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53418845%2ftrying-to-get-function-to-iterate-matrix%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown

























0






active

oldest

votes








0






active

oldest

votes









active

oldest

votes






active

oldest

votes
















draft saved

draft discarded




















































Thanks for contributing an answer to Stack Overflow!


  • Please be sure to answer the question. Provide details and share your research!

But avoid



  • Asking for help, clarification, or responding to other answers.

  • Making statements based on opinion; back them up with references or personal experience.


To learn more, see our tips on writing great answers.




draft saved


draft discarded














StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53418845%2ftrying-to-get-function-to-iterate-matrix%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown





















































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown

































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown







Popular posts from this blog

MongoDB - Not Authorized To Execute Command

in spring boot 2.1 many test slices are not allowed anymore due to multiple @BootstrapWith

How to fix TextFormField cause rebuild widget in Flutter