Community Profile

# Torsten

Last seen: Today 2013 年からアクティブ

バッジを表示

#### Content Feed

Eigenvalues in symbolic matrix
Replace %% 2.stability analysis J=jacobian([gp(1)*vp(1)./(kp(1)+N0)*N*P1 - mp(1)*P1^2==0, ... gp(2)*vp(2)./(kp(2)...

| 採用済み

sum of Bessel function
m = 0:250; x = linspace(-5*pi,5*pi,1000).'; A = 2./((2*m+1)*pi).*besselj(0,pi/2*m).*sin(1000*(2*m+1).*x-(2*m+1)*pi/2); f = su...

Finite difference temperature distribution with TDMA
You solve the system iteratively. But you were told to solve it with the Thomas-Algorithm. So apply this algorithm to the matrix...

| 採用済み

How to solve a system of 8 equations having integrals!
d0 = ones(1,8); d = fsolve(@fun,d0) norm(fun(d)) function res = fun(D) d=1.18; E=27000000; L=121.53; I=(pi*d^4)...

| 採用済み

In fmincon I set the linear constraint x1+x2+x3=12, but the sum of decision variables in the iteration result does not meet the condition
Tighten your ConstraintTolerance, and you'll see that fmincon converges to an infeasible point for n=12 and n=13. 1st： (main.m)...

1日 前 | 2

| 採用済み

unable to plot RungeKutta trajectory graph, potentially initial conditions?
Your code produces NaN values for z(1,:) and z(3,:), most probably caused by thetad and thetav. alpha = pi/4; tend = 2; v0 = ...

1日 前 | 0

bvp4c/bvp5c error
Did you create a function with name "eps" ? Then you should rename it because "eps" is an internal MATLAB function.

1日 前 | 0

How to solve f'(x)=0?
You want to find x symbolically or numerically ? Numerically use "fzero". Symbolically use "solve".

1日 前 | 0

Trying to find area under my plot
https://de.mathworks.com/matlabcentral/answers/514107-how-to-find-area-enclosed-by-a-closed-curve To get the positive and negat...

2日 前 | 0

I want to use a matrix and a function to then create another function
theta1 = -50; theta2 = 30; theta3 = 15; theta4 = -75; theta5 = -45; theta6 = -60; DH = [theta1, 486.5, 150, -90;... ...

2日 前 | 0

Solving system of 9 nonlinear equaitons in 16 variables
x0 = -10*ones(16,1); AB = [13 15 -1;9 24 0;7 2 35]; options = optimset('TolFun',1e-16,'TolX',1e-16); x = fmincon(@(x)fun(x,AB...

2日 前 | 2

the same value resulted from the linear optimization problem
You define w(k)=beta(k)*(-k_(k) *w_k(k)-k_(k)*log2(lamda(k)*(1-tau^2))*Pt + 2*k_(k)*log2(sqrt(segma_squared))-v_k(k)); and w...

2日 前 | 0

Hello guys, I need to solve this system of nonlinear symbolic equation; the variables known are X,Y,Z and I already know the solution for q1; so I just need to find q2 and q3.
syms d1 d5 a2 a3 a4 q1 q2 q3 X Y Z % Square eq1 and eq2, add both sides of the equations and take sqrt % Check for other solu...

2日 前 | 0

| 採用済み

What's the simplest way to determine the mass of a sphere from a plot of its density varying in r?
If r and rho are column vectors where r is ordered as r(1) < r(2) < ... < r(n) and rho is the density at position r, an approxim...

2日 前 | 2

| 採用済み

Index exceeds the number of array elements. Index must not exceed 1.
Pb is a scalar. So Pb(ii) does not exist. You can use Fbump(ii)=Pb/VA(ii); %this line of code is givig me the error but I don...

2日 前 | 0

Erro in calling function using horner()
I don't understand what you intend by using the lines [pz,b] = horner(a,x); [dpz,b]= horner(b,x); [pz ,a] = horner(a,x) The...

2日 前 | 0

Solving an integral when the parameters is a matrix - Controllability Gramian
You mean syms t t0 tf A = [0 1 0 0;0 0 -4.91 0;0 0 0 1;0 0 73.55 0]; B = [0;2;0;-10]; [M,J] = jordan(A); phi_s=M * expm(J...

2日 前 | 0

| 採用済み

sum of Bessel function
syms x syms m integer A = 2/((2*m+1)*pi)*besselj(0,pi/2*m)*sin((2*m+1)*1000*x-(2*m+1)*pi/2) g = symsum(A,m,0,Inf) fplot(g...

2日 前 | 0

Solving Laplace equation using Finite difference
If L = 3 and H = 1, you cannot choose the same number of points in x and y direction if you update U as if you use a grid with e...

2日 前 | 0

| 採用済み

solve matrix differential equations with ode45
syms t y1(t) y2(t) A = [0 1;8 -2]; Y = dsolve(diff([y1;y2],t)==A*[y1;y2]) y1 = Y.y1; y2 = Y.y2; diff(y1,t)-y2 diff(y2,t)-8...

3日 前 | 0

How do I use meshgrid and quiver3 in a system of differential equations?
https://web.sha1.bfh.science/Labs/F2/ODE/VectorFields.pdf

3日 前 | 0

How can I solve this matrix question using for loop?
Not the most efficient solution, but maybe best to understand what's happening. A=[0 0 0 ;6 5 7;12 14 30;35 40 42;75 55 50;85 6...

3日 前 | 0

| 採用済み

myfunction not working when applied on a meshgrid
x= 0:0.01:2*pi; %629 element. y= 0:0.001:1; %1001 element. [X, Y]= meshgrid(x,y); z=zeros(1001,629); %% for X and Y are bot...

3日 前 | 1

| 採用済み

How to find the actual value of d,b,t and D in the maximum value of sac?
close all clear all clc % Constants c = 340; % Speed of sound rho = 1.204; % ...

3日 前 | 0

Need to solve this with Runge-Kutta 4th order method in MatLab.
fy=@(x,y,z) z; fz=@(x,y,z) (1-y^2)*z-y; x(1)=0; z(1)=0; y(1)=1; h=0.1; xfinal=10; N=ceil((xfinal-x(1))/h); for j=1:N ...

3日 前 | 0

| 採用済み

If statement should continue for some values even if "else" is true
You must have a flag which indicates whether the engine is on or off. function [w,T] = EnginePower(B_SoC) persistent on ...

4日 前 | 0

Make a loop with a function to vary one value
tspan = [0 :0.01:20]; y0 = 0; opt = odeset('RelTol',1e-12,'AbsTol',1e-16,'MaxStep',0.001); Eta = [10, 50, 100, 500, 1000]; f...

4日 前 | 0

| 採用済み

Verify the Divergence theorem in MATLAB
Nice exercise. syms X Y Z x y z syms a phi positive F = [X^2*Y^2,Y^2*Z^2,Z^2*X^2]; divF = diff(F(1),X)+diff(F(2),Y)+diff(F(...

4日 前 | 0

How to get the derivate using bvp4c
I assumed Alpha = A in your code. clc; p=0.01; Betaf= 207; Betas = 17; Beta = 0.5; kof = 0.613; kos = 400; m = 1; b2 = ...

4日 前 | 0

| 採用済み

How can I more fully vectorize to eliminate for loops?
Wans2 = (sum(G.*(W.'-V),2)).'

5日 前 | 0

| 採用済み