
VBBV
PEC
Programming in MATLAB applied to interdisciplinary engineering fields
Statistics
RANK
82
of 257,819
REPUTATION
1,612
CONTRIBUTIONS
25 Questions
1,014 Answers
ANSWER ACCEPTANCE
84.0%
VOTES RECEIVED
151
RANK
5,074 of 17,766
REPUTATION
225
AVERAGE RATING
3.00
CONTRIBUTIONS
12 Files
DOWNLOADS
77
ALL TIME DOWNLOADS
2200
CONTRIBUTIONS
0 Posts
CONTRIBUTIONS
0 Public Channels
AVERAGE RATING
CONTRIBUTIONS
0 Highlights
AVERAGE NO. OF LIKES
Content Feed
phase shift using xcorr
t = 1 :0.1: 20; phaseshift= pi/2; s1 = sin(2*t); s2 = sin(2*t+ phaseshift); figure(1); clf; plot (t...
3 days ago | 0
I am trying to find the integral gain (ki) of the following transfer function.
syms ki ki = 10 % e.g numeric values , not symbolic sys = tf([30.672 175.694*ki],[1 38.672 175.694*ki]); S = stepinfo(sys); ...
3 days ago | 0
Hold on with subplot in a loop
A={[40;50;60] [70;80;90] [20;30;50]}; B={[45;55;65] [75;85;95] [25;35;55]}; T=[1;2;3]; figure(1) % use this outside of l...
4 days ago | 0
Question about the inability to display 16k octave bands(1/3 oct band)?
if fs==22050 fc=1000*2.^((-6:1:2)); fctxt={'16','32','63','125','250','500','1k','2k','4k'}; else fc=1000*...
6 days ago | 1
| accepted
Error using FDE12
s = 1:4; % define s vector passed as argument to function line below e.g. values sdot = @(t,s)[A-(e*s(1)*s(2))-(a*s(1)*s(3))+...
7 days ago | 0
An error in my code
clear syms x % syms F= atan(x); t5= taylor(F,x,'order',5); t7= taylor(F,x,'order',7); fplot(x^3/3+x,[-pi pi],'b') hold o...
8 days ago | 2
I do some revisions to my equation but I still got no answer. It says that there is no enough input argument. Can someone please help me?
x0=[150;150;150;150;150;150;150]; % x = 1:7; % give input vector to fucntion sol=fsolve(@rle,x0) % call the function using fu...
9 days ago | 0
| accepted
Pcolor Plot returning error Color data input must be a matrix.
h = pcolor(X1, Y1, z); % use X1, Y1 instead of x1 , y1, z must be 21x26 matrix
9 days ago | 0
how to plot solution of ODE eqution with out using fplot ?
% this code is used to show how to use ode in matlab, via a simple problem % for serios RC circuit with DC supply clc;clear a...
11 days ago | 1
| accepted
How can I store the whole diagonal and not only the first element (For loop)
% Part 3.1: Diagonal %A rng('default') r2 = randi(100,3,3); % Creating a 3x3 matrix, with random values from 1 to 100 %B ...
11 days ago | 0
| accepted
Why is ySol(t) = Dsolve(ode,cond) wrong? What should it be instead?
y=eval(ySol(t))% Use ySol(t) to evaluate
11 days ago | 0
Having trouble writing a multiple gravitational physics simulation code, as every planet effects one another.
a(q,:)=(sum(F(q:(q+o-1),:)))./mass(q); %Calculate acceleration of each planet v(q,:)=v(q,:)+a(q,:)*ts; %calula...
12 days ago | 0
Trouble graphing this non-linear function
x = 1:100; y = ((factorial(x)).^(1./x))./x; % use element wise division fplot(y,[0 100]) % no space between function and par...
12 days ago | 0
I am trying to a series generated from complex function but it gives to me an error.
clear; clc; syms k Delta=0; chi=0; theta=0; alpha=5; phi=0; bn=exp(-0.5*alpha^2)*(alpha^k)/sqrt(gamma(k+1)); bn1=exp(-0...
18 days ago | 0
using for loops to calculate compound interest with yearly contributions
N=17804; P=1000; C=1000; r=(10/100); I = 1; while P<=N y(I)=(P*(1+r)); % use an index for desired amount, P=y(I)...
23 days ago | 0
Semilogx function giving a linear scale
S = [1e-2 1e-1 1 10 100 1000 1e4]; % Remdesivir concentration, in uM Km = 0.5; % Equilibrium constant, in uM Vmax = 0.82; % ma...
24 days ago | 0
| accepted
Need help with syntax and how to use sigma in matlab
clear close all Z1=500 Z2=50 Z3=500 u=1 tau1=(2*Z2)/(Z1+Z2) tau2=(2*Z1)/(Z1+Z2) tau3=(2*Z3)/(Z3+Z2) p1=(Z2-Z1)/(Z2+Z1...
24 days ago | 1
| accepted
Unable to perform assignment because the size of the left side is 1-by-2 and the size of the right side is 2-by-2. (Line 32)
pp=@(t,x)[.23*x(1)-.0133*x(1)*x(2); -.4*x(2)+.0004*x(1)*x(2)] [t,y]=rk4sys(pp,[0,65],[550,20],0.3) figure(1) plot(t,y(:,1),...
25 days ago | 1
How to skip a data set in the legend of a plot?
clear; clc; clf; close all; x1 = [0.015 0.020 0.050 0.100 0.200 0.300 0.400 0.500]; y1 = [5.518E-11 7.215E-11 1.832E-10 ...
25 days ago | 1
| accepted
how to put for loop value into a vector
k = 0.168:0.168:1.68 ; for x=1:length(k) if k(x)>=0.03 && k(x)<=0.06 h(x)=1.437*k(x)^(-0.279); elseif k(x)>=0.07 && k(x)<...
25 days ago | 0
fplot invalid parameter error
m=100; k1=1500; k2=2500; k3=700; wb=4.4; Y=0.05; w1=sqrt(k1/m); w2=sqrt(k2/m); w3=sqrt(k3/m); r1=wb/w1; r2=wb/w2; r...
26 days ago | 1
Error with the command Geoplot
lat=cell2mat(A.Latitud);% if Latitud is the header name of column. Please check this. There Was extra parenthesis in my pr...
26 days ago | 0
Multiplying f(x) with f(-x)
cof = [1 15 7 2] syms x F(x) F = (cof(1)*x.^3+cof(2)*x.^2+cof(3)*x+cof(4)) FF = subs(F,x,-x) Fx = F.*FF % resulting produc...
26 days ago | 0
| accepted
Discontinuous display in Simulink Scope
Select the plot type as Auto in scope block. This plots the graph as continuous signals instead of discrete ones which you have...
27 days ago | 0
Why is my if counter not working
s= 0 x = linspace(-1,1, 50) n = linspace(1,10,50) L = 2 total = 0 s= zeros(size(n)) for i = 1:length(n) s(i) = (8/(p...
27 days ago | 0
I need to plot a graph of P, Rdia and Vo, please help me how to get graph. (for every 10 increment of pressure, there is 1 ohm increment of Rdia , and display the output vo)
Rdia = 3673; % dia resistance value in ohm ; Rtemp = 3618; % Temperature resistance value in ohm; R1 = 3600; % fixed resistanc...
27 days ago | 1
| accepted
Export results from several excel files into one excel file
Partikel_RL(xy,:)=Elementanzahl_RL(xy,:); % second for loop is not needed end
29 days ago | 0
Legend wont display lines properly in figure window
lgd = legend('Bundle','30% Opening','50% Opening', ... '65% Opening','85% Opening','100% Opening'); Use without a { } ...
1 month ago | 0
gamma function error in calculation
syms x a Y=sym(zeros(1)); Y(1)=0; a=1/2 for i=1:4 if i==5 A=1; B = 0; Y(i+1)=simplify((gamm...
1 month ago | 0
Function works with one input but not the other
F(2,:) = Thrust ./ MFR ; F(1,:) = (y .* A_choke) ./ MFR; Check the size of output matrix
1 month ago | 0